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NOMENCLATURE 


2 
A         surface  area  of  wire  (in.  ) 


A         vertical  distance  from  zeroline  to  n   current  peak  as 


n 


m 


measured  on  photographs  (in.) 


a  defined  by  a  =  -(1/C),   (farads   ) 

B  defined  by  B  =  R/2L,   (see"  ) 

C.  capacitance  of  RC  integrator  used  in  current  measuring 
circuit  (farads) 

C  heat  capacity  (joules/gram.  K) 

c  velocity  of  light  (meters/sec) 

E  energy  deposited  in  exploding  wire  (joules) 

E  initial  energy  stored  in  capacitor  bank  (joules) 

E  increase  in  energy  density  above  standard  conditions  needed 


to  melt  wire  material  (eV/atom) 


E         estimated  integration  error  by  Simpson's  rule 

E         increase  in  energy  density  above  standard  conditions  needed 
to  vaporize  wire  material  (eV/atom) 

I1        energy  density  required  to  ionize  electrons  from  wire 
material  (eV/atom) 

I         current  calibration  (kA/photo  division) 

i         current  (kA) 

J  (x)     Bessel  function  of  the  first  kind  of  order  zero 
o 

J1 (x)     Bessel  function  of  the  first  kind  of  order  one 

1/2 
K        defined  by  K  =  V  (C/L)     ,   (ampere) 

K  defined  by  equation  (F-ll) 

k  thermal  conductivity  (joules/meter«sec •  K) 

L  inductance  of  discharge  circuit  (henries) 

P  power  (rate  of  energy  deposited  in  wire) ,   (megawatts) 
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NOMENCLATURE   (continued) 

2 
p         pressure  (dyne/cm  ) 

Q         rate  of  conductive  and  convective  heat  losses  from  wire 
(joules/sec) 

Q'        energy  deposited  per  unit  length  of  wire  (joules/in.) 

R         resistance  of  discharge  circuit  (ohms) 

R  ,      calibrating  resistance  (ohms) 
cal 

R         extraneous  circuit  resistance  (ohms) 
e 

R.        resistance  of  RC  integrator  used  in  current  measuring  circuit 
(ohms) 

R         resistance  of  wire  (ohms) 
w 

r         radial  distance  from  center  of  wire  (in.) 

r         radius  of  wire  (in.) 
o 

T         absolute  temperature  (  K) 


t.        time  at  end  of  j    time  interval  after  initiation  of  discharge 


t  time  (sec) 

J  (sec) 

t  maximum  time  interval  of  interest  during  a  wire  explosion  (sec) 
m 

U  internal  energy  of  exploding  wire  (joules) 

u(x-x  )    unit  step  function 
o 

V  voltage  on  capacitor  bank  (volts) 

V  volt/cm  setting  on  amplifier  L   (see  Figure  6) 
S  P 

V  volt/cm  setting  on  amplifier  L   (see  Figure  6) 

L  S 

S 

V  initial  charging  voltage  on  capacitor  (volts) 
W  hydrodynamic  energy  of  the  wire  (joules) 


NOMENCLATURE  (continued) 

w         weight  of  n   value  (used  for  averaging) 

X         defined  by  equation  A-10  (dimensionless  ratio) 

Y         defined  by  equation  A-24  (dimensionless  ratio) 

y.        value  of  function  at  end  of  i   time  interval 

Z         atomic  number 

z         defined  by  z  =  In  V(t) 

2 
a         thermal  diffusivity  (in. /sec) 

At        time  interval  between  measured  points  in  measuring  voltage 
across  calibrated  resistor  (sec) 

AU.       internal  energy  change  at  j    transition  (joules) 

e  emissivity 

0(r,t)  temperature  variable  defined  by  equation  (F-2) 

u  permeability  (rnhcm   ) 

tt  3.141592654 

p  density 

o  electrical  conductivity  (ohm  «m  ) 

a  Stefan-Bolt zmann  constant  =  5.7  x  10   erg/sec»cm  •  K 

a(x)  estimated  standard  deviation  of  single  observation  for  variable  x 

2 

a  (x)     estimated  variance  of  single  observation  for  variable  x 

a(x)      estimated  standard  deviation  of  mean  for  variable  x 

2  — 
a  (x)     estimated  variance  of  mean  for  variable  x 

t        period  of  oscillation  of  discharge  circuit  (sec) 

$        ratio  of  current  at  wire  center  to  steady  state  current  (defined 
by  equation  G-l) 

1/2         -1 
w         frequency,  u  =  (1/LC)      ,   (sec   ) 


1.0   Introduction 

Exploding  wires  have  been  used  for  many  practical  purposes  as 
well  as  being  studied  because  of  basic  scientific  interest  in  the 
phenomena  itself.   Examples  of  practical  applications  are  high  tem- 
perature sources  for  chemical  reactions,  metal  forming  explosives, 
and  as  initiators  of  (a)  explosive  devices,  (b)  shock  waves,  and 
(c)  short  duration  intense  light  sources.   In  addition,  some   re- 
search has  been  aimed  at  the  possibility  of  initiating  a  fusion  pro- 
cess with  an  exploding  wire. 

The  high  energy  densities  associated  with  an  exploding  wire 
indicate  that  they  could  be  used  to  simulate  a  nuclear  explosion. 
Peaceful  applications  of  nuclear  explosions  include  mining,  exca- 
vation, energy  sources,  isotope  production,  and  in  contribution  to 
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basic  scientific  knowledge.      However,  the  amount  and  the  chemi- 
cal form  of  the  fallout  resulting  from  a  nuclear  explosion,  even  a 
typical  underground  shot,  prohibits  their  application  as  noted  above. 
It  is  desirable  for  all  applications  of  nuclear  explosions  to  control 
the  chemical  form  and  the  amount  of  radioactive  fallout. 

In  a  typical  underground  shot  experimental  evidence  indicates 
that  water  vapor  and  radioactive  gasses  are  vented  from  the  cavity 
between  thirty  seconds  and  two  minutes  after  the  shot;  this  is  the 
time  interval  of  interest  in  simulating  the  chemical  effects  so  that 
the  radioactive  products  can  be  tied  up  in  nonvolatile  and  insoluable 
form.   Thus,  the  temperature  for  the  simulation  of  the  chemical  en- 
vironment that  determines  the  volatile  products  of  a  nuclear  explo- 
sion is  a  few  thousand  degrees  and  the  pressure  is  about  that  of  the 


static  soil  load  above  the  shot.   It  is  under  these  conditions  in 
the  environment  of  the  explosion  that  the  chemical  nature  of  the 
fallout  materials  is  determined.   Meyer  and  Svensson  of  the  Lawrence 
Radiation  Laboratory  have  shown  through  experimental  work  that  ex- 
ploding wires  may  be  used  to  simulate  the  pressure  and  temperature 
conditions  that  are  important  in  determining  the  chemistry  of  the 

o 

nuclear  explosion.    However,  suitable  apparatus  needs  to  be  devel- 
oped which  will  trap  the  pressure  and  heat  energy  released  in  the 
wire  explosions  to  do  meaningful  work. 

The  objectives  of  the  research  reported  here  were  to: 

1)  determine  the  characteristics  of  a  capacitor  bank  discharge  into 
an  exploding  wire  including 

(a)  the  effect  of  wire  size  and  composition  on  energy  deliv- 
ered and  the  rate  of  energy  delivery. 

(b)  the  effect  of  circuit  parameters  -  capacitance,  induc- 
tance, and  charging  voltage. 

(c)  designing  of  adequate  fast  switching  circuitry. 

2)  determine  the  amount  of  energy  and  rate  of  deposition  of  energy 
needed  to  produce  temperatures  and  pressures  analogous  to  those  of 
a  nuclear  explosion. 

3)  design  and  determine  the  effect  of  a  sample  cell  (including  mag- 
netic pinch)  on  the  energy  density. 

4)  determine  the  feasibility  of  simulating  the  chemical  environment 
of  a  nuclear  explosion  by  exploding  wire  techniques. 

5)  compare  alternative  possibilities  for  simulating  a  nuclear  ex- 
plosion. 


2.0  Literature  Survey 
2.1  Mechanism  of  Wire  Explosion  -  Physical  Description 

The  first  paper  on  exploding  wires  was  read  before  the  Royal 
Society  of  London  in  1773  by  Edward  Nairne.    At  that  time  equip- 
ment was  not  available  for  detailed  study  of  the  microsecond  phe- 
nomena.  Significant  quantitative  research  on  exploding  wires  was 
not  reported  until  after  1950,  when  microsecond  measuring  techniques 
were  sufficiently  developed.   The  development  of  these  techniques 
made  possible  accurate  and  quantitative  measurements  of  current, 
voltage,  and  other  factors  associated  with  this  high  speed  phenom- 
enon. 

Since  1950  a  large  amount  of  interest  has  developed  in  exploding 
wires,  sometimes  as  a  means  to  an  end  and  sometimes  for  the  sake  of 
the  phenomena  itself.   Exploding  wires  are  used  in  many  different 
fields,  some  of  which  are  commercial.   Thus,  although  exploding  wire 
phenomena  is  a  relatively  new  area  of  quantitative  scientific  study, 
there  is  much  information  available  in  the  literature  on  the  charac- 
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teristics  of  exploding  wires. 

A  complete  explanation  is  not  available  to  describe  the  exploding 
wire  mechanism.   An  equation  of  state  is  needed  but  the  necessary 
data  are  not  available  to  formulate  one.   It  is  generally  agreed  that 
the  explosion  of  a  wire  is  a  complex  phenomenon  resulting  from  the 
deposition  of  a  large  amount  of  energy  in  a  very  short  period  of  time 
(microseconds  or  less)  into  a  piece  of  metal  of  small  cross  sectional 


area. 


A  rapidly  increasing  current,  initially  controlled  by  the  RLC 


circuit  parameters,  passes  through  the  wire.   After  a  brief  heating 

interval  the  wire  resistance  is  much  greater  than  the  resistance  in 
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the  rest  of  the  circuit.   The  heavy  current  density  (-10  -10   A/in.  ) 

heats  the  wire  material  to  its  melting  point  (See  Fig.  1).   (The  atmos- 
pheric melting  point  is  usually  taken  to  be  the  temperature  at  which 
the  wire  melts.   )   The  molten  wire  is  heated  further  by  the  current 
which  may  be  increasing  or  decreasing  at  this  time.   To  discuss  ener- 
getic wire  explosions  (simulating  a  nuclear  explosion) ,  let  us  consider 
only  the  increasing  current  case  (Fig.  la).   Since  the  entire  exploding 
wire  process  occurs  in  microseconds  or  less  the  physical  shape  of  the 
wire  is  maintained  by  inertial  and  magnetic  pressures.   The  wire  con- 
tinues to  be  heated  by  the  current  and  reaches  the  normal  boiling  point 
near  the  current  maximum.   However,  pressure  very  significantly  affects 
the  boiling  point  and  equilibrium  boiling  is  not  expected.   Superheating 
occurs  as  the  temperature  of  the  wire  rises   past  the  atmospheric 
boiling  point  while  the  wire  material  is  still  in  the  liquid  state. 

At  some  time  t  ,   with  the  wire  temperature  above  the  normal  boil- 
ing  point,  the  voltage  across  the  wire  rapidly  increases  to  a  peak  value 
and  the  wire  (wholly  or  partially  in  a  vapor  state)  actually  explodes. 
The  controlling  mechanism  near  the  time  of  burst  is  poorly  understood, 
however  several  mechanisms  have  been  proposed  including  (1)  an  expan- 
sion process,  (2)  a  shock  process,  and  (3)  an  explosive  phase  change. 
During  a  wire  explosion  a  considerable  amount  of  noise  and  an  intense 
flash  of  light  are  observed.   One  or  more  strong  shock  fronts  are  gen- 
erated during  this  same  time  interval. 

Conn  had  proposed  the  following  six  steps  for  the  mechanism  of  a 
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wire  explosion. 
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Figure  la.   Fast  Wire  Explosion 


Figure  lb.   Slow  Wire  Explosion 


Figure  1.   Curves  for  typical  fast  and  slow  wire  explosions 


1.  Heating  of  the  wire  by  the  heavy  current  passing  through  it. 

2.  Formation  of  a  liquid  column  of  metal  replacing  the  wire. 

3.  Development  of  instability  in  the  wire  which  forms  unduloids 
(figures  of  rotation)  and  causes  the  appearance  of  striation 
in  the  metal  vapor. 

4.  Upon  disruption  of  the  arcs  formed  between  the  remaining 
unduloids,  the  so  called  'dark  interval'  is  obtained  during 
which  no  current  flows  through  the  wire,  the  voltage  remaining 
constant. 

5.  A  sudden  flash  of  light  is  next  observed,  the  spectrum  of 
which  is  continuous  and  independent  of  the  material  used. 
Thereafter  emission  and  absorption  lines  are  observed  depending 
on  the  materials  used. 

6.  One  or  several  shock  fronts  are  observed  during  the  same  time 
interval. 

It  is  known,  however,  that  certain  types  of  wire  explosions 
exhibit  definite  deviation  from  the  above  behavior.   Steps  3  and  4  of 
Conn's  mechanism  are  strongly  dependent  on  the  rate  of  energy  delivery 
to  the  wire.   Thus  it  has  been  suggested  by  Chace  and  Levine  that 
rather  than  describing  a  general  wire  explosion  as  Conn  has  done  it  is 
preferable  to  classify   the  explosion  based  on  the  rate  of  energy  de- 
livery to  the  wire.   On  this  basis  Chace  and  Levine  have  suggested  the 

following  four  classes  of  wire  explosions:   (1)  melting,   (2)  slow 

13 
explosive,   (3)  fast  explosive,  and   (4)  explosive  ablation.    Webb 

14 
has  also  classified  exploding  wires  by  the  rate  of  energy  delivery. 

The  energy  deposited  in  the  wire  for  the  melting  classification  is 
less  than  that  necessary  to  vaporize  the  wire  completely.   The  wire 
simply  behaves  as  a  common  fuse  by  melting  and  falling  apart. 

In  a  slow  explosion  instabilities  have  time  to  develop  before 
sufficient  energy  to  vaporize  the  wire  has  been  deposited.   The  explo- 
sion occurs  slowly  enough  that  these  physical  distortions  affect  the 
mechanism  appreciably. 


For  a  fast  explosion  the  energy  to  vaporize  the  wire  is  deposited 
before  significant  changes  occur  in  the  physical  shape  of  the  wire. 
Unduloids  do  not  form  before  the  wire  explodes  and  thus  the  explosion 
mechanism  is  not  affected  by  physical  distortions.   If  the  explosion 
is  fast  enough,  inertial  effects  hold  the  wire  in  the  same  physical 
shape  until  after  enough  energy  has  been  deposited  to  vaporize  the 
wire. 

For  an  explosive  ablative  type  explosion  the  energy  to  vaporize 
the  wire  is  deposited  in  a  very  short  time  interval.   In  this  case  the 
skin  effect  described  by  Bennett   and  the  magnetic  pressure  described 
by  Katzenstein   control  the  explosion  mechanism.   Appendix  G  presents 
a  discussion  of  skin  effect  and  magnetic  pressure. 

In  a  slow  wire  explosion  segments  of  the  wire  between  the  undu- 
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loids  vaporize  to  produce  a  rapid  rise  in  resistance.   '     For  a  fast 

explosion  an  expansion  wave  proceeds  into  the  center  of  the  wire  which 

18  ^   , 
also  causes  a  large  increase  in  resistance.    The  increase  in  resis- 
tance in  either  case  can  result  in  a  'dark  interval'  during  the  explo- 
sion process.   The  dense  wire  vapor  is  then  nonconducting  until  it  has 
expanded  enough  for  the  voltage  remaining  across  the  wire  to  cause 
breakdown  by  arcing  through  the  vapor  and/or  air.   If  the  voltage  remain- 
ing across  the  gap  when  the  resistance  approaches  its  peak  is  large 
enough   the  arc  breakdown  will  occur  immediately  or  even  before  the  wire 
has  completely  vaporized  and  no  'dark  interval'  will  be  observed.   If 
the  voltage  remaining  across  the  wire  after  the  first  pulse  is  very  low 
the  wire  will  not  become  conducting  again  or  in  the  vernacular  of  explo- 
ding wire  technology,  there  will  be  no  restrike  after  the  onset  of  cur- 
rent dwell. 


2.2  Parameters  Affecting  a  Wire  Explosion 

Parameters  associated  with  the  energy  source,  the  exploding  wire, 
the  circuit  switch,  the  transmission  lines,  and  the  ambient  medium  are 
all  important  in  the  mechanism  of  a  wire  explosion  and  in  determining 
the  energy  deposited  in  the  wire.   A  capacitor  bank  is  usually  used  as 

the  energy  source  but  a  coaxial  cable  has  been  used  by  a  Sandia 

19 

group.     The  advantages  of  the  coaxial  cable  as  an  energy  source  are: 

(1)  the  current  approximates  a  step  function,  (2)  the  discharge  circuit 
equations  are  algebraic  rather  than  nonlinear  differential  equations, 

and  (3)  the  coaxial  symmetry  provides  stabilizing  effects  on  the  wire 

19 
and  also  provides  good  frequency  response.    Since  a  capacitor  bank 

was  readily  available  and  a  suitable  coaxial  cable  was  not,  the  capac- 
itor bank  was  used  in  this  experimental  investigation. 

The  parameters  of  interest  in  a  capacitor  bank  energy  storage 

system  are  the  charging  voltage  and  the  capacitance.   Since  the  energy 

2 
stored  in  a  capacitor  bank  is  equal  to  1/2  CV  ,  both  the  charging 

voltage  and  the  capacitance  should  be  large  to  obtain  the  desired  energy 

storage.   The  internal  resistance  and  inductance  of  the  capacitors 

should  be  small  in  order  to  produce  a  fast  frequency  response  as  well  as 

to  reduce  energy  loss  to  a  minimum. 

Since  transmission  lines  also  affect  exploding  wires,  specification 

of  transmission  characteristics  is  necessary  to  interpret  results  and 

20 

to  make  possible  accurate  duplication  of  experiments.     Therefore,  the 

transmission  characteristics  should  be  controlled  and  specified  in  each 
exploding  wire  experiment.   Now  if  the  transmission  lines  are  long  the 
time  for  energy  to  propagate  from  the  capacitor  bank  to  the  wire  can  be 
of  the  same  magnitude  as  the  time  for  wire  explosions  and  interpretation 


of  results  and  accurate  reproducibility  of  wire  explosions  are  not 
easily  obtained.   If  transmission  lines  are  short  (the  preferred 
case),  the  circuit  inductance,  resistance,  and  capacitance  can  also 
be  treated  as  lumped  parameteres.   The  lumped  inductance  and  circuit 
capacitance  produce  a  natural  'ringing  frequency'  response  of  the 
system. 

Special  consideration  must  be  given  to  the  design  of  the 
switching  circuitry.   Large  energies  (approaching  a  megajoule)  and 
short  delivery  times  (less  than  a  microsecond)  are  required.  Switches 

which  have  been  used  in  exploding  wire  circuits  have  been  both  complex 

8  21  22 

and  simple.  '   '    The  switch  must  be  able  to  withstand  both  the  high 

voltages  applied  to  it  and  high  current  pulses  flowing  through  the 
switch.   In  addition  the  rise  and  jitter  times  of  the  switch  must  be 
less  than  a  microsecond. 

The  effect  of  wire  parameters  on  the  wire  explosion  have  been 
discussed  by  several  people.       It  is  generally  believed  that  a 
high  intensity  short  duration  energy  pulse  is  better  for  depositing 
maximum  energy  in  the  wire  than  a  long  duration  low  intensity  pulse. 
However,  Zinke  et  al.  suggest  that  in  certain  types  of  wires  a  rela- 
tively slow  mechanism  may  result  in  high  energy  densities  and  high  tem- 
peratures.  In  any  case  the  wire  parameters  must  be  matched  to  the  cir- 
cuit parameters  to  obtain  maximum  energy  deposit. 

The  wire  material  and  the  dimensions  of  the  wire  are  important  in 
the  mechanism  and  characteristics  of  a  wire  explosion.   Practical  con- 
siderations of  wire  material  include  ease  of  attachment,  mechanical 
strength,  toxicity,  and  compatibility  with  other  circuit  materials.   In- 
trinsic considerations  of  the  material  include  resistivity,  temperature 
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coefficient  of  resistivity,  density,  specific  heat,  melting  point, 
boiling  point,  heat  of  fusion,  heat  of  vaporization,  and  ionization 
energy. 

The  resistivity  of  an  exploding  wire  of  a  given  material  has  been 

shown  to  be  independent  of  the  current  density  (up  to  at  least 

8      2 
10  A/cm  )  or  energy  deposition  rate  up  to  the  energy  necessary  to  melt 

27 
the  wire.     The  measured  resistance  of  an  exploding  wire  agrees  with 

the  resistance  calculated  from  tabulated  data  for  temperatures  below  the 

28 
vaporization  temperature.     However,  the  resistance  is  not  a  unique 

function  of  energy  deposited  above  the  vaporization  energy  of  the  wire. 
An  increase  in  the  rate  of  electrical  energy  deposition  results  in  pro- 
gressively decreasing  slopes  of  the  resistance  versus  energy  curves.  The 

dependence  of  resistance  on  the  rate  of  energy  deposited  is  called  the 

29 

'resistance  anomaly'.     Bennett  has  shown  that  an  expansion  wave,  which 

converts  the  material  behind  it  into  a  nonconducting  state,  can  be  used 

18 
to  explain  the  resistance  anomaly  and  the  dark  interval.    The  expansion 

wave  proceeds  from  the  wire  surface  to  the  wire  center;  the  increase  in 
resistance  caused  by  the  expansion  wave  is  characteristic  of  the  resis- 
tance anomaly  and  the  dark  interval. 

Above  the  melting  temperature  the  resistance  increases  with  the 
introduction  of  energy  for  nearly  all  metals.   At  the  beginning  of  melt- 
ing there  is  a  rapid  rise  in  resistance  which  is  noticeable  on  the  vol- 
tage-current traces.   In  aluminum,  copper,  and  silver  wires   the  resis- 
tance increases  during  the  liquid  phase  with  increasing  energy  deposited. 
The  liquid  phase  of  platinum,  tungsten,  and  iron  have  small  temperature 
coefficients  of  resistivity.   Thus  the  liquid  phases  of  exploding  wires 
of  platinum,  tungsten,  and  iron  are  characterized  by  voltage  and  current 
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curves  which  give  a  resistance  almost  independent  of  temperature. 

Webb   et  al.  in  their  classification  have  proposed  that  metals  used 

14 
for  wires  fall  into  two  classifications.     Low  boiling  points  -  low  heat 

of  vaporization  materials  are  in  class  I.   High  boiling  point  -  high 

23 
heat  of  vaporization  materials  are  in  class  II.   Leopold   found,  in 

the  materials  he  studied,  an  approximately  constant  amount  of  energy  was 

deposited  in  the  optimum  length  of  wire  for  each  material.   If  this  is 

true  in  general,  class  I  materials  would  be  raised  to  higher  temperatures 

since  the  excess  energy  above  that  required  to  vaporize  the  wire  would 

be  used  to  heat  the  plasma.   Also  the  time  required  to  deposit  sufficient 

energy  to  vaporize  the  wire   is  less  for  class  I  materials   than  for 

class  II  materials.    With  this  shorter   energy  deposit   time  heat  losses 

by  conduction,   convection,   and   radiation  would  be  less   for  class   I 

materials. 

Energy  required  to  ionize  the  wire  material  does  not  increase  the 
temperature  of  the  wire  material.    Thus,   to  obtain  high  temperatures, 
a  wire  material  which  has  low  numbers  of  electrons  to  ionize  and   low 
ionization  energies  is  desirable. 

Leopold  has   shown   that   the   dimensions   of   the  wire  are   im- 
portant in  determining  the  amount  of  energy  deposited  in  the  wire,  the  peak 

23 

power,  and  the  time  of  burst.     Bennett  also  considered  optimum  size  wires 

in  matching  wire  resistance  to  the  circuit   parameters  to  obtain  maximum 

25 

energy  deposition.    The  time  to  burst  (time  interval  from  initiation  of 

discharge  to  the  voltage  peak)  decreases  with  an  increase  in  charg- 
ing voltage  of  the  capacitor  bank,  and  increases  with  increasing 
wire   diameter   if   other  parameters  are  constant.   Too    thin  wires   are 
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poorly  matched  to  the  discharge  circuit  because  they  explode  in  short 
times  and  use  little  of  the  stored  energy.   If  the  diameter  of  the 
wire  is  too  large,  the  energy  absorbed  by  the  wire  will  not  be  suffic- 
ient to  vaporize  the  wire.   The  best  time  reproducibility  is  obtained 
if  the  wire  explodes  on  the  initial  portion  of  the  first  current 
pulse.   For  maximum  power  and  energy  deposition  the  wire  should  ex- 
plode near  the  first  current  peak  (See  Fig.  la). 

The  power  and  the  energy  deposited  in  the  wire  also  depend  on 
the  wire  length.  There  is  an  optimum  length  for  the  maximum  energy 
deposition  and  peak  power.   However,  the  peak  power  per  unit  length 

and  the  energy  density  increase  with  decreasing  wire  length  if  end 

23 

effects  are  not  important. 

The  ambient  medium  in  which  a  wire  explodes  has  an  effect  on 
the  behavior  of  an  exploding  wire.   A  striking  example  is  the  sensi- 
tive dependence  of  the  pause  phenomena  on  the  identity  or  existence 
of  the  surrounding  gas.   The  burst  spike  is  flattened  and  broadened 

when  the  wire  is  confined  by  tubes  which  exceed  the  wire  diameter  by 

30 
a  factor  of  up  to  3  to  5. 

Energy  is  rapidly  transferred  to  material  surrounding  the  wire 

for  large  temperature  differences  between  the  wire  plasma  and  the 

surrounding  medium.   If  the  discharge  circuit  is  designed  so  that  the 

return  current  from  the  wire  to  the  capacitor  bank  is  distributed 

symmetrically  about  the  wire,  the  magnetic  pinch  should  confine  the 

■I  c. 

wire  material  for  a  longer  period  of  time.     (Appendix  G  discusses 
the  effect  of  magnetic  pinch.)   Longer  confinement  times  allow  the 
wire  to  be  driven  to  higher  temperatures.   Thus,  if  the  magnetic  pinch 
confines  the  wire  material,  energy  would  be  transferred  more  rapidly 
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to  the  material  surrounding  the  wire  than  if  the  wire  material  is 

not  confined.   Results  of  Sakurai  and  Takao  indicate  that  applied 

31 
axial  magnetic  fields  may  confine  exploding  wires. 

Energy  is  deposited  in  material  surrounding  exploding  wires 

during  wire  explosions.   Spectral  lines  of  the  confining  material 

have  been  observed  in  the  emission  spectra  obtained  during  the  time 

30 

interval  of  wire  bursts.     Lower  rates  of  blackbody  cooling  than  ex- 
pected were  observed  unless  a  fairly  large  amount  of  deposited  energy 

had  been  transferred  to  a  more  massive  body  than  the  wire  (such  as 

30 

the  surrounding  material) .     Thus  blackbody  cooling  at  lower  rates 

than  expected   indicates  that  energy  is  deposited  in  the  material 
surrounding  the  wire. 
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3.0  Theory 

A  quantitative  analysis  of  an  exploding  wire  must  consider  the 
electrical  circuit,  the  wire,  and  the  ambient  medium.   For  a  capacitor 
discharge  into  a  wire  the  electrical  circuit  can  be  described  by 
Kirchoff's  laws.   The  law  of  energy  conservation  can  be  used  to  describe 
the  thermal  heating  of  the  wire  and  the  energy  losses  due  to  conduction, 
convection,  radiation,  hydrodynamic,  and  other  forms  of  energy. 

Assuming  the  capacitance,  the  resistance,  and  the  inductance  of 
the  circuit  can  be  treated  as  lumped  parameters,  the   circuit   equation 
can  be  written  as 


L  (di/dt)  +  R  i  +  V  =  0  (1) 

c 

The  inductance  of  the  wire  (L  )  is  much  less  than  the  inductance  of 

w 

the  remaining  discharge  circuit  (L  ).   Therefore,  the  inductance 

(L  =  L  +  L  )  is  approximately  equal  to  L  .   After  a  brief  heating  in- 
cw  c 

terval  the  resistance  of  the  circuit  is  dominated  by  the  wire  resistance, 
Therefore,  R  can  be  considered  to  be  the  wire  resistance. 
The  law  of  energy  conservation  can  be  written  as 

(dU/dt)  =  P  -  a£A  T4  -  (dW/dt)  -  Q  (2) 


where: 

U  is  the  internal  energy  in  the  wire, 

P  is  the  rate  at  which  energy  is  deposited  in  the  wire, 

A 
aeAT  is  the  radiative  energy  losses, 

dW/dt  is  the  hydrodynamic  energy  losses,  and 
Q  is  the  conductive  and  convective  heat  losses. 
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Thus,  the  amount  of  internal  energy  in  the  wire  depends  on  the  power 
(i.e.,  the  rate  at  which  energy  is  deposited  in  the  wire)  and  the  en- 
ergy losses.   At  some  time  after  the  energy  for  vaporization  has  been 
deposited  in  the  wire  the  liquid  metal  expands  and  forms  vapor.   It 

has  been  estimated  that  metallic  conduction  should  cease  for  expansions 

32 
of  2-5  times.    With  no  current  flowing  there  is  no  power  input  to 

the  wire.   Thus,  in  a  wire  explosion  it  is  only  the  energy  added  in  the 
short  time  interval  after  the  wire  has  exceeded  the  vaporization  tem- 
perature and  before  it  begins  to  expand  that  is  effective  in  producing 
extremely  high  temperatures  (~  10   K) .   The  rate  of  radiative  energy 
losses  at  high  temperatures  is  of  course  very  large  (~  10   watts) . 
Thus,  high  energy  densities  and  high  temperatures  in  the  wire  can  only 
be  obtained  by  depositing  the  energy  from  the  electrical  circuit  very 
rapidly  (in  microseconds  or  less) . 

The  maximum  rate  of  energy  delivery  to  the  wire  occurs  if  the 
wire  is  well  matched  to  the  circuit  parameters.   The  wire  is  said  to  be 
well  matched  to  the  circuit  parameters  if  the  short  time  interval  after 

the  wire  has  exceeded  the  vaporization  temperature  and  before  it  begins 

32 
to  expand,  occurs  near  the  current  peak  (See  Fig.  la).   Bennett    has 

shown  that  a  wire  is  well  matched  to  the  circuit  parameters  for  an 

optimum  mean  value  of  wire  resistance  which  satisfies  the  relation: 

1.1   <  R   ..  (C/L)1/2     <  1.3  .  (3) 

—  opt  — 

Before  the  current  peak  the  exploding  wire  circuit  is  character- 
ized by  a  simple  RLC  circuit  and  is  described  by  equation  (1) .   During 
the  time  period  from  the  start  of  the  discharge  until  the  current  peak, 
the  energy  which  was  originally  stored  in  the  capacitor  bank  is  trans- 
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ferred  to  ohmic  heating  of  the  wire  and  to  storage  in  the  circuit  in- 
ductance.  In  a  well  designed  exploding  wire  circuit  at  the  time  of 
the  current  peak  the  voltage  on  the  capacitor  bank  should  drop  to 
approximately  zero  and  the  energy  originally  in  the  capacitor  bank 
should  be  completely  transferred  to  ohmic  heating  of  the  wire  and  to 
storage  in  the  magnetic  field.   Since  the  voltage  on  the  capacitor 
bank  is  approximately  zero  near  the  current  peak,  equation  (1)  can  be 
approximated  by 

■  L  (di/dt)   +   R  i  =  0  .  (4) 

If  the  wire  is  well  matched  to  the  circuit  parameters,  the  wire 
material  exceeds  the  vaporization  temperature  near  the  current  peak  at 
which  time  the  energy  stored  in  the  magnetic  field  is  rapidly  deposited 
in  the  wire.   During  the  time  interval  when  the  wire  exceeds  the  vapor- 
ization temperature  to  the  time  of  wire  burst,  the  wire  resistance  is 
rapidly  increasing.   An  RL  circuit  described  by  equation  (4)  has  a  time 
constant  of  (L/R) .   Therefore  since  the  wire  resistance  is  large  just 
before  the  wire  burst,  the  energy  stored  in  the  magnetic  field  is  rapidly 
deposited  in  the  wire.   Therefore,  to  obtain  high  energy  densities  and 
high  temperatures  in  an  exploding  wire   it  is  desirable  to  have  some  of 
the  energy  stored  initially  in  the  capacitor  bank  used  for  ohmic  heating 
of  the  wire,  and  the  majority  of  the  energy  from  the  capacitor  bank 
transferred  to  the  magnetic  field  of  the  circuit  inductance. 

The  equations  describing  the  exploding  wire  circuit  could  be 
solved  analytically  if  the  resistance,  the  cross  sectional  area,  and  the 
specific  heat  of  the  wire  were  constant.   However,  these  factors  depend 
on  the  physical  state  of  the  wire,  and  for  wire  temperatures  above   the 
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vaporization  temperature   they  also  depend  on  the  rate  of  energy  deliv- 
ery to  the  wire.   An  equation  of  state  is  required  for  each  of  the  dif- 
ferent types  of  wire  explosions  if  the  equations  are  to  be  solved  to 
obtain  definite  results.   In  addition,  solutions  of  equations  require 
as  a  corallary  a  relationship  between  resistance  and  the  physical  state 
of  the  wire.   Neither  an  equation  of  state  nor  a  quantitative  relation- 
ship for  the  wire  resistance  is  available.   Thus,  a  combined  analytic 
and  experimental  approach  must  be  used  over  limited  temperature  ranges 
and  limited  time  intervals  to  solve  the  equations. 

The  wire  resistance  is  calculated  using  the  relation 

_        V   -   L  (di/dt)  ,_. 

R   =     w      w  (5) 

w     -. 


where: 


R  is  the  wire  resistance, 
w 

V  is  the  voltage  across  the  wire, 

w  ° 

L   is  the  inductance  of  the  wire, 
w 


i  is  the  current  through  the  wire,  and 

(di/dt)  is  the  time  rate  of  change  of  current. 

V  ,  i,  and  (di/dt)  are  measured  experimentally.   L   can  be  calculated 
w  r  J  w 

33 

from  the  physical  properties  of  the  wire. 

The  energy  deposited  in  the  wire  material  at  any  time  t  after  the 
start  of  the  discharge  of  the  capacitor  bank  can  be  computed  by  inte- 
grating the  instantaneous  power  over  the  time  interval  from  the  start  of 
the  discharge  until  time  t. 
t 


E   = 
w 


R  i2  dt'   .  (6) 

w 
o 
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Assuming  the  inductance  of  the  wire  is  constant:  and  substituting 
equation  (5)  into  equal  ion  (G)  gives 

r   t 


v 


( 


V  i  dt.1   )   -  1/2  L  i  (t)   . 

W  W 


(7) 


Thus,  the  energy  placed  in  the  wire  can  be  computed  by  integrating 
the  product  of  the  voltage  and  current,  and  subtracting  the  energy 
stored  in  the  inductance  of  the  wire  gap  at  time  t. 

The  internal  energy  of  the  wire  is  related  to  the  temperature  of 
the  wire  by  the  relation. 


U(T) 

dU« 
U(To) 


T*  -  6 


j-1   6->0 


r  xj 


(dU'/dT)  dT  +   AU.   u(T-T.)  ]  +  [lim 


T.  -+5 
J-1 


6->0 


(dU'/dT)  dT  ] 


T  +  6 
n 


(8) 


where: 

T   is  the  initial  temperature  of  the  wire, 
o 

T.;  j=  1,2,...,  n  are  transition  temperatures, 

U(T)  is  the  internal  energy  and  is  a  function  of  temperature, 

•th      .  _.      j 
AU  is  the  internal  energy  change  at  the  j    transition,  and 

J 
u(T-T.)  is  the  unit  step  function. 

An  upper  limit  of  the  temperature  reached  by  an  exploding  wire 
can  be  determined  by  equating  the  energy  deposited  in  the  wire  with 
the  change  in  internal  energy  of  the  wire.   The  energy  deposited  in 
the  wire  can  be  computed  using  equation  (7).   Substitution  of  this  cal- 
culated energy  into  equation  (8)  permits  a  determination  of  the  upper 
temperature  limit.   (This  upper  limit  is  calculated  assuming  no  energy 
losses  in  equation  (2).) 

If  radiation  losses  are  considered  in  equation  (2),  another 
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estimate  of  the  maximum  possible  wire  temperature  is  obtained.   For  a 
fast  wire  explosion  the  energy  is  "dumped"  in  the  wire  in  less  time 
than  the  instability  time  constant.   Since  the  maximum  temperature  is 
reached  before  the  wire  has  expanded,  the  hydrodynamic  losses  (dW/dt) 
are  negligible.   However,  at  extremely  high  temperatures  the  radiation 
energy  losses  are  much  greater  than  the  conductive  and  convective 
energy  losses  (See  Appendix  F) .   Thus  when  large  amounts  of  energy  are 
deposited  in  fast  wire  explosions,  high  temperatures  are  produced.   For 
such  explosions,  equation  (2)  can  be  approximated  by 

(dU/dt)   =  P  -ea  A  T4   .  (9) 

The  maximum  internal  energy  and  the  maximum  temperature  are 
reached  when  the  rate  of  change  of  internal  energy  is  zero  (dU/dt  =  0). 
Thus,  from  equation  (9)  the  maximum  power  is 

P    =  eaA  T4     .  (10) 

max         max 

Thus,  to  obtain  high  temperatures  in  an  exploding  wire  it  is  necessary 

to  have  high  peak  powers  and  small  radiating  areas.   An  estimate  of  the 

maximum  power,  P    ,  delivered  to  a  constant  resistive  load  from  an  RLC 
max 

._  .   34 
circuit  is 

Pmax=  0/2/2) (C/L)1/2   .  (11) 

Substituting  equation  (11)  into  equation  (10)  and  solving  for  T    gives 

IH.3.X 

Tmax  =  t  (C/L)1/2   (V*/2  eaA)  ]lM   .  (12) 

Therefore,  to  obtain  high  temperatures  in  an  exploding  wire,  it  is 
necessary  to  have  high  initial  charging  voltage  on  the  capacitor, 
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appreciable  capacitance,  and  low  inductance. 

To  study  the  high  temperature  chemistry  of  the  material 
surrounding  an  exploding  wire,  it  is  necessary  to  transfer  the  energy 
from  the  wire  to  the  surrounding  material.   Little  energy  is  lost 
from  the  wire  material  by  any  method  other  than  radiation  until  after 
the  wire  bursts,  and  little  energy  is  deposited  in  the  wire  after 
about  one  microsecond  following  the  wire  explosion.   If  energy  is  de- 
posited instantaneously  along  a  line,  a  thermal  spike  may  be  used  to 
represent  the  system.   Thus,  the  energy    deposited  in  a  system  con- 
sisting of  a  thin  central  wire  and  a  low  thermal  conductivity  surround- 
ing medium,  e.g.>  earth,  can  be  represented  by  a  thermal  spike.   If  Q' 

is  the  energy  released  per  unit  length  of  wire,  the  temperature  distri- 
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bution  in  the  surrounding  medium  is  as  follows 

T(r,t)  =  T   +   (Q'/4ttC  p)(l/at)  exp(-r2/4at)  (13) 

where : 

T(r,t)  is  the  time  dependent  temperature  distribution, 
T  is  the  temperature  before  energy  is  deposited, 
C  is  the  heat  capacity  of  the  surrounding  material, 
p  is  the  density  of  the  surrounding  material,  and 
a  is  the  thermal  diffusivity  of  the  surrounding  material. 
Thus,  the  time  dependent  temperature  distribution  is  just  a  function 
of  Q'  and  the  physical  properties  of  the  surrounding  material. 

Typically,  materials  such  as  those  in  the  earth's  crust  have 
small  thermal  conductivities.   Therefore,  essentially  all  of  the  energy 
would  be  deposited  in  a  wire  surrounded  by  earth  before  any  energy 
would  be  conducted  into  the  surrounding  material.   Hence,  the  thermal 
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spike  can  be  used  to  predict  the  time  dependent  temperature  distri- 
bution in  a  sample  of  earth's  crust  surrounding  an  exploding  wire. 

The  temperature  or  the  energy  density  necessary  to  study  high 
temperature  chemistry  or  the  formation  of  a  cavity  in  a  sample  ma- 
terial is  usually  known  or  can  be  estimated.   Knowing  the  temperature 
necessary  to  achieve  the  desired  effect,  the  necessary  parameters  of 
the  system  can  be  estimated  as  follows. 

1)  Calculate  the  energy  which  must  be  deposited  in  the  wire  to 
raise  the  desired  amount  of  sample  material  to  the  desired  temperature. 
To  complete  this  calculation  handbook  or  estimated  values  of  thermal 
diffusivity,  density,  and  heat  capacity  for  the  sample  material  are 
used  in  equation  (13).   Typical  results  are  shown  in  Fig.  2. 

2)  Calculate,  using  equation  (8),  the  wire  temperature,  at  which 
the  wire  material  will  contain  the  energy  calculated  in  step  (1). 
Handbook  values  of  the  physical  properties  are  used  for  the  calculation 
when  they  are  available.   For  the  extremely  high  temperatures  typical 
of  wire  explosions,  the  specific  heat  and  latent  heats  must  usually  be 
estimated  since  they  are  considerably  beyond  handbook  values.   The  re- 
sults shown  in  Figure  3  were  calculated  assuming  an  ideal  plasma  gas 
and  that  the  i   electron  was  ionized  from  each  atom  when  the  kinetic 
energy  (3kT/2)  of  the  gas  particles  was  equal  to  the  ionization  energy 
I.,  i  =  1,2,*«',Z.   Ionization  energies  tabulated  by  Moore   were  used 

for  the  calculations.   (Rouse  has  calculated  more  accurate  results  by 
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using  the  Saha  equation.   ) 

3)  Determine  the  rate  of  energy  input  to  the  wire  (power)  which 
will  balance  the  radiation  losses  from  the  wire  at  the  necessary  wire 
temperature.   The  power  corresponding  to  a  given  temperature  can  be  cal- 
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Figure  2.  Temperature  distributions  in  a  cylindrical  thermal  spike  at 
various  times.  Parameters  typical  of  the  earth  crust  were„ 
used.   C  =.2  cal/g«°K  ;  k=.005  cal/seccm-  K;   p=2.24  g/cm  , 
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Figure  3.   Energy  per  unit  length  as  a  function  of  wire  temperature 
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culated  using  equation  (10);  typical  results  are  shown  in  Fig.  4. 
(The  product  of  the  wire  diameter  and  the  emissivity  of  the  wire  ma- 
terial is  used  as  a  parameter.) 

4)  Determine  the  electrical  circuit  parameters  necessary  to 
obtain  the  power  calculated  in  step  (3).   The  maximum  theoretical 
power  available  from  an  RLC  network  to  a  constant  resistive  load  can 
be  estimated  using  equation  (11).   Figure  5  shows  the  estimated  max- 
imum theoretical  power  available  from  an  RLC  network  as  a  function  of 
charging  voltage  with  (C/L)  as  a  parameter.   This  maximum  power  is 
available  to  an  exploding  wire  if  the  wire  is  well  matched  to  the  cir- 
cuit parameters.   In  determining  the  electrical  circuit  parameters,  it 

is  necessary  that  the  initial  energy  stored  by  the  capacitor  bank 

2  j 

(E   =  1/2  CV  )  be  greater  than  the  energy  calculated  in  step  (1) . 
co         o 
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4.0   Experimental 

4.1  Apparatus 

A  schematic  diagram  of  the  discharge  circuit,  the  voltage 
measuring  circuit,  and  the  current  measuring  circuit  is  shown  in 
Fig.  6.   The  discharge  circuit  consisted  of  the  specimen  wire,  the 
air-gap  switch,  the  capacitor  bank,  the  transmission  lines,  and  the 
high  voltage  power  supply.   The  capacitor  bank  consisted  of  8  capac- 
itors (General  Electric  -  Pyranol)  of  nominal  1  microfarad  capaci- 
tance and  a  maximum  charging  voltage  of  25  kilovolts. 

The  air-gap  switch  was  similar  to  the  one  described  by  Harvey 

22 
and  Yevick.     (Nichrome  wire  was  used  as  the  trigger  electrode, 

which  differs  from  the  molybdenum  wire  used  by  Harvey  and  Yevick.) 
Harvey  and  Yevick  give  the  following  characteristics  of  their  spark 
gap  switch:   operational  voltage  range  -  20-35  kV,  peak  operational 
current  -  100  kA  at  500  kHz,  maximum  peak  current  -  greater  than  200 
kA,  closing  and  jitter  times  -  less  than  0.1   usee,  and  switch  induc- 
tance -  50  nH  with  parallel  plate  return  path.   During  an  experiment 
after  the  high  voltage  power  supply  (North  American  Philips  Co.,  Inc.  - 
X-ray  Supply  Unit)  had  been  used  to  charge  the  capacitors,  the  power 
supply  was  disconnected  from  the  basic  discharge  circuit.   The  power 
supply  was  then  used  to  supply  a  voltage  pulse  to  the  trigger  elec- 
trode of  the  air-gap  switch.   This  voltage  pulse,  applied  through  small 
blocking   capacitors,  rendered  the  air-gap  conducting  and  initiated  the 
discharge  of  the  capacitor  bank  through  the  v/ire.   (The  capacitor  bank 
and  x-ray  supply  unit  were  borrowed  from  the  KSU  Physics  Department.) 
Massive  brass  and  aluminum  leads  were  used  as  transmission  lines 
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Figure  6.   Schematic  diagram  of  measuring  circuits 
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to  keep  the  extraneous  circuit  resistance  and  inductance  low.  Specimen 
wires  were  placed  in  small  holes,  drilled  in  the  ends  of  5/16"  diameter 
leads;  the  wires  were  fastened  securely  to  the  leads  with  set  screws. 

The  circuits  used  to  measure  the  voltage  across  the  wire  and  the 
current  through  the  wire  are  shown  schematically  in  Fig.  6.   The  char- 
acteristics of  the  scope  and  plug  in  amplifiers  obtained  from  the 
Tektronix  manuals  showed  that  the  frequency  response  of  the  scope  and 
plug  in  amplifiers  was  adequate  in  the  range  from  DC  to  18  MHz  and  the 
risetime  was  20  nsec.   High  voltage  probes  used  in  this  study 
(Tektronix  P6013)  had  a  maximum  peak  pulse  limit  of  12  kV  below  100  kHz 
and  a  risetime  of  14  nsec.   It  should  be  noted  that  the  voltage  measur- 
ing circuits  ground  had  no  physical  connection  with  the  discharge  cir- 
cuit.  Operating  under  this  condition,  ground  loops  in  the  voltage 
measuring  circuit  were  prevented.   The  high  voltage  probes  were  attached 
to  the  brass  leads  close  to  the  wire.   The  signals  from  the  high  vol- 
tage probes  were  then  led  to  the  input  terminals  of  a  differential  am- 
plifier G   (Tektronix  type  G)  powered  by  channel  one  of  a  P-127 
Tektronix  power  supply.   The  output  of  channel  one  of  the  power  supply 
was  then  applied  to  amplifier  G   (Tektronix  type  G)  which  was  used  as  a 
variable  gain  amplifier.   Amplifier  G  was  in  the  upper  channel  of  a 
dual-beam  oscilloscope  (Tektronix  type  551) . 

The  current  was  measured  using  a  small  coil  placed  near  a  current 
lead  of  the  discharge  circuit.   The  voltage  induced  in  this  coil  was 
proportional  to  the  time  derivative  of  the  current  through  the  lead.  The 
voltage  signal  from  the  coil  was  integrated  by  a  simple  RC  integrating 
network  with  a  time  constant  (RC  =  5  microseconds)  large  with  respect 
to  the  most  rapid  oscillation.   The  integrated  signal,  which  was  pro- 
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portional  to  the  current  in  the  discharge  circuit,  was  amplified  by 
amplifier  L   (Tektronix  type  L)  placed  in  the  lower  channel  of  a  dual 
beam  scope. 

The  time  derivative  of  the  current  was  measured  using  a  sep- 
arate pickup  coil.   The  voltage  on  the  coil  was  attenuated  by  ampli- 
fier L   in  the  second  channel  of  the  P-127  power  supply.   The  output 
of  channel  two  was  then  applied  to  amplifier  L  which  was  in  the  lower 
channel  of  the  dual  beam  scope. 

Oscilloscope  curves  of  either  voltage  and  current  or  voltage  and 

current  derivative  were  measured  simultaneously.   Note  that  amplifier 

L  was  used  in  both  the  current  measuring  circuit  and  the  current 
s  ° 

derivative  measuring  circuit.   The  curves  were  recorded  on  3000  speed/ 
type  47  Polaroid  film.   A  Hewlett-Packard  oscilloscope  camera  (Model 
196A)  was  used  to  take  the  photographs.   The  instantaneous  values  of 
voltage,  current,  and  current  derivative  waveforms  were  determined 
from  the  photographs.    A  Leitz   Coordinate  Comparator, 
available  in  the  Applied  Mechanics  Department,  was  used  to  measure 
distances  on  the   photographs.   The  current  and  current  derivative 
were  related  to  the  same  time  base  by  aligning  prominent  voltage  peaks 
or  discontinuities  of  the  separate  voltage  records. 

4.2  Calibration  Procedure 

The  parameters  of  interest  in  an  exploding  wire  experiment  have 
been  discussed  in  Section  2.0;  the  power  and  the  energy  deposited  in 
the  wire  were  the  parameters  of  specific  interest  in  this  study.   The 
power  and  the  energy  were  calculated  with  the  equations  in  Section  3.0 
using  measured  values  of  voltage,  current  and  current  derivative.   Be- 
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fore  taking  the  necessary  exploding  wire  data,  the  measuring  circuits 
were  calibrated  and  the  effect  of  the  air-gap  switch  determined.  The 
experimental  sequence  was  as  follows: 

1)  The  voltage  measuring  circuit  and  the  static  voltmeter  were 
calibrated. 

2)  The  circuit  capacitance  was  measured. 

3)  The  circuit  inductance  and  resistance  were  measured. 

4)  The  current  measuring  circuit  was  calibrated. 

5)  The  effect  of  the  air-gap  switch  was  determined. 

6)  The  current  derivative  measuring  circuit  was  calibrated. 

The  vertical  deflection  of  the  oscilloscope  photographic  records 

for  the  voltage  measuring  circuit  was  calibrated  using  two  high  voltage 

power  supplies  in  series  (Fluke,  Model  412B  and  Fluke,  Model  400BDA) . 

The  maximum  available  calibration  voltage  was  7  kilovolts.   The  Fluke 

power  supplies  had  been  previously  calibrated  using  a  digital  voltmeter 

(Electronic  Associates,  Inc.  TDVM  5001).   The  deviation  of  the  vertical 

deflection  of  the  voltage  measuring  circuit  was  within  the  reading 

accuracy  of  the  oscilloscope  traces  for  each  combination  of  amplifier 

settings  (G  and  G  )  used  in  succeeding  measurements.   The  oscilloscope 
P      s 

was  then  used  to  calibrate  the  electrostatic  voltmeter  (Sensitive 
Research  Instrument  Corporation,  Model  ESH)  used  to  measure  the  charging 
voltage  on  the  capacitor  bank.   The  fractional  standard  deviation  of  the 
electrostatic  voltmeter  was  .011  volts  per  volt. 

The  capacitor  bank  with  either  eight  or  four  capacitors  in  parallel 
was  used  to  explode  wires.   The  capacitance  of  each  combination  was 
determined  by  measurements  of  voltage-time  oscillograms  obtained  by  dis- 
charging the  capacitor  bank  through  calibrated  resistors.   The  calibrated 
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resistors  were  placed  in  the  location  which  was  used  for  the  specimen 
wires.   The  detailed  procedure  used  in  measuring  and  calculating  the 
capacitance  is  given  in  Appendix  A.   The  following  values  were  calcu- 
lated for  the  capacitance  at  an  initial  charging  voltage  of  8  kilovolts. 

Number  of  Total  Capacitance 

Capacitors  (microfarads) 

4  3.587  +  .042 

8  7.712  +  .051 

The  extraneous  resistance  (resistance  of  the  circuit  with  no  wire 
resistance)  and  the  circuit  inductance  were  determined  by  placing  a 
heavy  shorting  bar  (5/16"  diameter  brass  rod)  in  the  position  of  the 
specimen  wire.   Photographic  records  of  the  current  through  the  circuit 
and  the  voltage  on  the  capacitor  bank  were  obtained.   These  records 
were  measured  using   the  Leitz  Coordinate   Comparator.    From 
these  measurements  the  inductance  and  the  extraneous  circuit  resistance 
were  calculated  using  the  equations  of  an  underdamped  RLC  circuit  (See 
Appendix  A) .   The  same  records  were  used  to  calibrate  the  vertical  de- 
flection on  the  oscilloscope  records  for  the  current  measuring  circuit 
as  described  in  Appendix  A.   The  pickup  coil  which  measured  current 
through  the  discharge  circuit  was  placed  near  a  lead  of  the  discharge 
circuit.   After  calibration,  this  pickup  coil  of  the  current  measuring 
circuit  was  fixed  in  position. 

The  vertical  deflection  of  the  current  derivative  measuring  cir- 
cuit on  the  oscilloscope  traces  was  also  calibrated  with  the  shorting 
bar  in  place.   Underdamped  RLC  circuit  equations  similar  to  the  current 
waveform  calibration  were  used  for  the  calibration  of  the  current 
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derivative  waveform  (See  Appendix  A) . 

The  effect  of  the  air-gap  switch  was  determined  by  measuring 
the  circuit  resistance.   Initial  charging  voltages  of  5,  8,  11,  and 
15  kilovolts  on  the  capacitors  were  used.   The  distance  of  the  air- 
gap  switch  was  also  varied,  the  extraneous  circuit  resistance  in- 
creasing with  an  increase  in  the  air-gap  distance  for  constant  initial 
voltage.   The  circuit  resistance  decreased  for  the  same  air-gap  set- 
ting with  an  increase  in  initial  charging  voltage.   Variation  in  circuit 
resistance  due  to  atmospheric  conditions  was  of  the  same  order  of  mag- 
nitude as  the  variation  due  to  the  switch  setting.   Since  the  circuit 
resistance  was  small  (~  0.035  ohms)  compared  to  the  wire  resistance 
(except  for  a  short  time  interval  after  the  initiation  of  the  discharge) 
neither  the  air-gap  setting  nor  atmospheric  conditions  significantly 
affected  the  circuit  resistance  during  the  explosion  of  the  wire. 

Moses  and  Korneff   have  suggested  that  induced  voltage  in  the 
high  voltage  probes  may  contribute  significantly  to  the  voltage  measure- 
ments.  This  was  checked  by  shorting  the  tips  of  the  high  voltage  probes 
together  with  a  shorting  bar.   The  probes  were  then  placed  near  their 
actual  position  when  measuring  the  voltage  across  a  wire,  and  a  wire 
was  exploded.   (The  probes  had  no  physical  connection  with  the  discharge 
circuit  while  checking  for  induced  voltages.)   It  was  determined  that  the 
induced  voltage  was  not  significant.   However,  the  position  of  the  oscil- 
loscope with  respect  to  the  discharge  circuit  was  important.   If  the 
scope  was  close  to  the  discharge  circuit,  significant  induced  voltages 
were  noted.   Thus,  for  exploding  wire  measurements  the  scope  was  position- 
ed so  that  stray  signals  were  not  monitored. 
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4.3   Exploding  Wire  Measurements 

Voltage,  current,  and  current  derivative  were  measured  for  bare 
wires  exploded  at  atmospheric  conditions.   The  calibrated  circuits  shown 
in  Fig.  6  were  used  for  the  measurements.   The  wire  material,  wire  size, 
circuit  capacitance,  and  circuit  charging  voltage  were  varied.   Table  I 
shows  the  wire  and  circuit  parameters  used  in  this  study.   The  experi- 
mental sequence  used  to  obtain  exploding  wire  data  was  as  follows: 

1)  A  wire  of  interest  was  inserted  in  the  holes  in  the  ends  of  the 
brass  leads  and  fastened  securely  with  set  screws.   The  brass  leads 
with  high  voltage  probes  attached  were  then  adjusted  so  the  \\dre  was 
straight  between  the  leads. 

2)  Leads  to  amplifier  L  were  attached  for  either  the  current  measur- 
ing  circuit  or  the  current  derivative  measuring  circuit.   The  vertical 
scale  of  the  voltage  measuring  circuit  and  the  vertical  scale  of  either 
current  or  current  derivative  circuit  were  set  to  desired  values.   The 
time  scale  on  the  oscilloscope  was  set  to  a  desired  value.   The  scope 
was  set  for  single  sweep  operation.   (Scale  settings  were  determined  by 
a  trial  procedure  until  the  magnitudes  of  a  particular  wire  explosion 
were  determined.) 

3)  The  safety  ground  on  the  capacitor  bank  was  disconnected  and  the 
capacitor  bank  was  charged  to  the  desired  voltage  using  the  high  voltage 
power  supply.   (If  the  capacitor  bank  discharged  through  the  switch  be- 
fore reaching  the  desired  voltage,  the  electrodes  of  the  air-gap  switch 
were  set  farther  apart.) 

4)  The  camera  shutter  was  opened  for  a  time  exposure. 

5)  Circuit  discharge  was  initiated  by  a  voltage  pulse  to  the  trigger 
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electrode  of  the  air-gap  switch.   (If  the  capacitor  bank  did  not  dis- 
charge through  the  switch  when  a  voltage  pulse  was  applied,  the  capac- 
itor bank  was  grounded  and  the  switch  air-gap  distance  was  decreased.) 
The  scope  was  triggered  by  the  signal  from  the  voltage  measuring  cir- 
cuit. 

6)  The  safety  ground  was  allowed  to  spring  back  and  ground  the  capac- 
itor bank.   The  safety  ground  prevented  a  charge  from  remaining  or 
accumulating  on  the  capacitor  bank. 

7)  The  oscilloscope  was  retriggered  to  obtain  a  zeroline  reference 
for  each  photograph. 

8)  The  camera  shutter  was  closed. 

9)  Portions  of  the  wire  still  connected  to  brass  leads  were  removed. 

10)  The  camera  film  was  developed  if  desired.   (More  than  one  set  of 
traces  were  obtained  on  one  picture  by  adjusting  the  multiple  picture 
knob  on  the  camera.) 

The  consistency  of  results  was  determined  by  exploding  eleven 
standard  wires  on  two  successive  days.   The  wire  and  circuit  parameters 
were  not  varied  for  the  standard  wire  explosions.   Iron  wire  14  mils  in 
diameter  and  one  inch  long  was  used  as  the  standard.   Circuit  conditions 
used  for  the  standard  wire  were  8  microfarads  and  11  kilovolts.   The 
following  values  were  calculated  for  the  standard  wire  explosions  (The 
error  limits  reported  are  estimated  standard   deviations   for   single 
observations . ) : 
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Table  I.   Data  taken  for  exploding  wires 

Wire  Parameters  Circuit  Parameters* 
Wire  Material        Wire  Size 

Diameter   Length  Capacitance    Voltage 

(in.)      (in.)  (yF)          (kV) 

(Wires  exploded  at  ambient  atmospheric  conditions) 


Aluminum 

.032 

1 

Aluminum 

.032 

1 

Copper 

.010 

1 

. Copper 

.010 

1 

Copper 

.010 

1 

Copper 

.010 

1,3 

Copper 

.010 

1,3 

Copper 

.020 

1/2,1,2 

Copper 

.020 

1/2,1,2,3 

Iron 

.014 

1 

Iron 

.014 

1 

Iron 

.014 

1 

Iron 

.014 

1,3 

Iron 

.014 

1 

Nickel- 

-  Chromium 

.028 

1 

Nickel- 

•Chromium 

.028 

1 

8 

15 

4 

11 

4 

15 

8 

8 

8 

11 

8 

15 

8 

11 

8 

15 

4 

11 

4 

15 

8 

8 

8 

11 

8 

15 

8 

11 

8 

15 

(Wires  exploded  in  sample  cell  geometry) 

Copper  .010       1  8  11 

Copper  .020       1  8  11 

Iron  .014       1  8  11 

^Circuit  inductance  for  all  cases  was  1.33  microhenrys . 
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Time  to  burst  3.04  +   .22  microseconds, 

Peak  Voltage  12.91  +   .55  kilovolts, 
Peak  Resistance  1.90  +   .21  ohms, 

Peak  Power  91.   +  11.  megawatts,  and 

Energy  Deposited  119.    +  15.  joules. 

Data  for  other  experiments  conducted  at  ambient  atmospheric 
conditions  as  listed  in  Table  I  were  then  measured.   For  each  set  of 
conditions  listed  in  Table  I  a  minimum  of  two  sets  of  voltage- 
current  curves  and  one  set  of  voltage-current  derivative  curves  were 
determined  (3  wire  explosions).   The  records  taken  with  the  oscillo- 
scope camera  were  measured  using  the  coordinate  comparator.  To  deter- 
mine if  the  discharge  and  measuring  circuits  were  operating  consistent- 
ly a  standard  wire  was  exploded  periodically. 

After  analyzing  the  bare  wire  results  data  were  obtained  for 
wires  exploded  in  the  sample  cell  geometry.   The  experimental  sequence 
used  for  wires  exploded  in  the  sample  cell  geometry  was  the  same  as 
for  the  bare  wires  except  the  sample  cell  was  included  in  the  discharge 
circuit.   Figure  7  shows  a  drawing  of  the  cross-section  of  the  sample 
cell.   (As  noted  earlier  this  cell  can  be  used  to  hold  material  in 
close  proximity  around  the  wire  during  the  explosion  and  as  discussed 
in  Section  2.0   the  ambient  medium  has  an  effect  on  exploding  wire 
characteristics.)   The  sample  cell  was  designed  to  use  the  magnetic 
pinch  principle.   Since  the  return  current  from  the  wire  to  the  capaci- 
tor bank  was  symmetrical  about  the  wire,  magnetic  pressure  tended  to 
confine  the  wire  material  after  vaporization  (See  Appendix  G) .   Longer 
confinement  times  should  result  in  more  energy  deposition  in  the  wires; 
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Figure  7.  Schematic  of  Sample  Cell  Geometry 
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the  attainment  of  high  energy  depositions  and  consequently  high  tem- 
peratures was  an  objective  of  this  research. 

Since  voltages  in  the  kilovolt  range  were  used  for  wire  explo- 
sions in  the  sample  cell,  voltage  breakdown  or  arcing  tended  to  occur 
between  points  of  the  current  path  in  the  sample  cell.   This  voltage 
breakdown  allowed  the  current  to  short  circuit.   Therefore,  the  pointed 
brass  lead  and  the  outside  of  the  cell  were  covered  with  electrical 
tape  to  prevent  arcing;  in  addition  an  insulator  cap  was  used  to  pre- 
vent arcing  from  the  brass  lead  to  the  cell  wall.   Three  mylar  disks, 
each  .010  inch  thick,  made  up  the  insulator  cap. 

Several  materials  (mylar,  mica,  polyethylene,  plexiglas,  teflon, 
and  compressed  CaC0„)  were  tested  for  use  as  an  insulator  cap.   However, 
the  energy  density  and  pressure  of  the  wire  explosions  either  melted 
or  shattered  all  tested  materials  except  mylar.   The  melting  or  shatter- 
ing of  the  insulator  cap  resulted  in  arcing  from  the  brass  lead  to  the 
cell  wall  and  thus  if  the  insulator  cap  melted  or  shattered  before  the 
time  of  burst,  the  amount  of  energy  deposited  in  the  wire  was  signifi- 
cantly reduced.   Since  maximizing  the  energy  trapped  in  the  material 
surrounding  the  wire  was  of  interest,  the  melting  or  shattering  of  the 
insulator  cap  after  the  time  of  burst  was  also  undesirable. 

One  inch  long  stainless  steel  tubes  were  used  for  the  sample  cell 
walls.   The  inside  diameter  of  the  tubes  used  were  1/16",  1/8",   and 
3/8".   To  provide  insulation  between  the  wire  and  the  cell  wall,   a 
silicon  based  caulk  was  placed  in  the  sample  volume  and  allowed  to  dry 
at  atmospheric  conditions.   If  sufficiently  high  temperatures  had  been 
produced  by  the  wire,  the  silicon  caulk  would  also  have  simulated  an 
environment  of  sand  around  an  underground  nuclear  explosion.   The  resid- 
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ual  acetone,  after  drying  the  caulk  at  atmospheric  conditions,  would 
have  simulated  hydrocarbon  vapor  formed  after  a  nuclear  explosion. 
The  caulk,  however,  was  insufficient  insulation  to  prevent  arcing 
from  the  central  wire  to  the  cell  v/all  for  the  two  smaller  diameter 
tubes.   Thus,  only  the  3/8"  diameter  tube  was  used  for  the  sample  cell 
measurements . 

Analysis  of  the  wire  explosions  at  atmospheric  conditions  showed 
that  the  two  sizes  of  copper  wire  and  the  iron  wire  were  better  matched 
to  the  circuit  parameters  than  the  other  wires  noted  in  Table  I.   There- 
fore, these  wires  were  exploded  in  the  sample  cell  geometry.   The  wire 
and  circuit  parameters  used  for  the  sample  cell  explosions  are  also 
listed  in  Table  I. 
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5.0  Results  and  Discussion 
5.1   General  Considerations  of  Exploding  Wire  Results 

The  parameters  which  significantly  affect  the  results  of  an  ex- 
ploding wire  experiment  are  discussed  in  Section  2.0.   Table  I  lists 
the  values  of  parameters  used  in  exploding  wire  experiments  conducted 
during  the  course  of  this  study.   (Wire  material,  wire  diameter,  wire 
length,  ambient  environment,  circuit  capacitance,  and  circuit  charging 
voltage  were  the  parameters  varied.)   Equations  developed  in  Section 
3.0  were  used  to  calculate  the  quantities  of  interest  for  the  analysis 
of  the  exploding  wire  data;  the  details  of  the  calculations  are  given 
in  Appendix  E.   Figures  8  through  33  show  voltage,  current,  resistance, 
power,  and  energy  as  functions  of  time  for  the  parameters  listed  in 
Table  I.   Error  bars  were  determined  but  were  not  indicated  on  the  al- 
ready complex  graphs.   (Errors  are  discussed  in  Section  5.2.) 

The  characteristics  of  the  voltage,  current,  resistance,  power, 
and  energy  curves  for  an  exploding  wire  have  been  examined  by  other 
investigators.   The  curves  of  Figs.  8  through  33  were  examined  to  de- 
termine if  the  wire  experiments  conducted  in  this  study  and  the  results 
of  other  investigators  (discussed  in  Section  2.0)  were  consistent. 
Quantitative  agreement  could  not  be  checked  since  the  discharge  circuits 
of  other  investigators  did  not  have  the  same  capacitance,  inductance, 
and  charging  voltage  of  the  circuits  used  in  the  present  study.   Thus, 
only  qualitative  agreement  with  the  results  in  the  literature  was  pos- 
sible for  the  voltage,  current,  resistance,  power,  and  energy  curves  of 
Figs.  8  through  33. 
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Figure  8.   Curves  for  32  mil  x  1  inch  aluminum  wire  exploded  at 
8  yF  and  11  kV. 
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Figure  10.   Curves  for  20  mil  x  1/2  inch  copper  wire  exploded  at 
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Figure  11.   Curves  for  20  mil  x  1  inch  copper  wire  exploded  at  8  yF 
and  11  kV. 
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Figure  12.   Curves  for  20  mil  x  2  inch  copper  wire  exploded  at  8  yF 
and  11  kV. 
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Figure  13.   Curves  for  20  mil  x  1/2  inch  copper  wire  exploded  at 
8  yF  and  15  kV. 
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Figure  14.   Curves  for  20  mil  x  1  inch  copper  wire  exploded  at 
8  yF  and  15  kV. 
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Figure  15.   Curves  for  20  mil  x  2  inch  copper  wire  exploded  at 
8  uF  and  15  kV. 
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Figure  16.   Curves  for  20  mil  x  3  inch  copper  wire  exploded  at 
8  uF  and  15  kV. 
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Figure  17.   Curves  for  10  mil  x  1  inch  copper  wire  exploded  at 
4  yF  and  11  kV. 
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Figure  18.   Curves  for  10  mil  x  1  inch  copper  wire  exploded  at 
4  yF  and  15  kV.  . 


E 
(joules) 


P 
(MW) 


I 
(kA) 


V 
(kV) 


R 
(ohms) 


500 


400 


300 


200 


100 


A  Energy 
o  Pov7er 


©  A  A  A  A 


C:: 


©  A 


0 


•_g ^a*_ 


45 


30 


15 


A  Current 

O  Voltage 
D  Resistance 


O 


O 


A 


A 


A 


A 


A 


AA 


A 


aaaaaaaaaaAaa 


'o  Q 


A 


w 


^JVo  Q  Q  Q  G)  QQ  OQ  QQ  3  Q  <p 


0 


oo°Q° 


A      "O 


p  r  r  -1  -  ^  5 


E^^B 


12  3 

Time  (microseconds) 


53 


Figure  19.   Curves  for  10  mil  x  1  inch  copper  wire  exploded  at 
8  yF  and  8  kV. 
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Figure  20.   Curves  for  10  mil  x  1  inch  copper  wire  exploded  at 
8  yF  and  11  kV. 
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Figure  21.   Curves  for  10  mil  x  3  inch  copper  wire  exploded  at 
8  uF  and  11  kV. 
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Figure   22.      Curves   for   10  mil  x   1   inch   copper  v/ire   exploded   at 
8  yF  and   15  kV. 
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Figure  23.   Curves  for  10  mil  x  3  inch  copper  wire  exploded  at 
8  uF  and  15  kV. 
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Figure  24.   Curves  for  14  mil  x  1  inch  iron  wire  exploded  at  4  pF 
and  11  kV. 
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Figure   25.      Curves   for   14  mil  x   1   inch   iron  wire   exploded   at   4  uF 
and   15  kV. 


60 


250 


200 


E 
(joules) 


P 
(MW) 


100 


A   Energy 
q      Power 


A  A 


A  A  A  A__ 


A 


A  A 


A 


°   A       -O 


^_A 


..'A 


3VQ 


Ofia 


c:  I.-';;-- 


o 


i  . 


> 


(kA)     15 


V 
(kV) 


R 

(ohms) 


10 


A   Current 
0   Voltage 

„   Resistance 


O 


A 


A    AAA 


O 


A 


A 


AA 


'AA 


AAA 


A       AAAA    oOOQn  A     ^^ 


O 


oo0 


Oq 


^A 


'O 


A 


Oq, 


AAAO    ° 


°ooooooo°A 


o 


AA 

Ogaooo0ooo 


A         °au    nLcocanDaGI]!n0acCDaDcaDGaGacC     I 


A 


fee 

1 


2  3 

Time  (microseconds) 


Figure  26.   Curves  for  14  mil  x  1  inch  iron  wire  exploded  at  8  yF  and 
and  8  kV. 
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material  for  longer  periods  of  time  after  the  vaporization  energy  is 
deposited.   Toward  this  end  a  sample  cell  employing  the  magnetic 
pinch  principle  (See  Appendix  G)  was  constructed  and  tested  to  determine 
its  effect  on  the  energy  density  of  wire  explosions. 

5.8  Results  of  Wire  Explosions  in  Sample  Cell 

Because  higher  temperatures  were  desired  than  obtained  from 
bare  wire  explosions,  a  sample  cell  was  designed  to  confine  the  wire 
material  after  vaporization  in  the  hope  this  would  lead  to  more  effec- 
tive and  uniform  heating  of  the  wire  material  with  consequently  higher 
energy  deposition.   In  the  bare  wire  studies  highest  energy  densities 
had  been  obtained  for  the  wires  well  matched  to  the  circuit  parameters; 
these  wires  included  the  14  mil  iron  wire,  the  10  mil  copper  wire,  and 
the  20  mil  copper  wire.   Since  deposited  energy  density  increased  with 
a  decrease  in  wire  length  and  end  effects  were  not  significant  for  one 
inch  wires  but  were  for  half-inch  wires,  one  inch  lengths  of  the  above 
wires  were  selected  for  study  in  the  sample  cell.   In  the  final  sample 
cell  design,  3/8  inch  diameter  stainless  steel  tubes  were  used  for  the 
cell  walls.   The  exploding  wire  was  located  at  the  center  of  the  tube 
and  insulated   from  the  walls  with  a  silicon  based  caulk.   Data,  taken 
with  the  10  mil  copper  wire  exploded  at  8  microfarads  and  11  kilovolts  in 
the  sample  cell,  indicated  that  arcing  between  the  wire  and  cell  walls 
had  occurred;  therefore  data  for  the  10  mil  copper  wire  were  not  reported, 

Figure  32  shows  the  results  of  a  14  mil  iron  wire  exploded  at 
8  microfarads  and  11  kilovolts  in  the  sample  cell.   Comparison  of  Fig.  32 
with  Fig.  27  shows  the  effect  of  the  sample  cell  geometry  on  a  wire  ex- 
plosion.  The  current  curves  differ  during  the  initial  portion  of  the 
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Figure   27.      Curves    for   14  nil   x   1   inch   iron  wire   exploded   at   8   uF  and 
11  kV. 
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Figure  28.   Curves  for  14  mil  x  3  inch  iron  wire  exploded  at  8  pF 
and  11  kV. 
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Figure  29.   Curves  for  14  mil  x  1  inch  iron  wire  exploded  at  8  uF  and 
and  15  kV. 
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Figure  30.   Curves  for  28  mil  x  1  inch  nickel-chromium  wire  exploded 
at  8  yF  and  11  kV. 
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Figure  31.   Curves  for  28  mil  x  1  inch  nickel-chromium  wire  exploded 
at  8  yF  and  15  kV. 
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Figure  32.   Curves  for  14  mil  x  1  inch  iron  wire  exploded  in  sample 
cell  geometry  at  8  yF  and  11  kV. 
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Figure  33.    Curves  for  20  nil  x  1  inch  copper  wire  exploded 
in  sample  cell  geometry  at  8  yF  and  11  kV. 
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The  resistance  curves  were  examined  to  see  if  they  were  consis- 
tent with  the  known  physical  properties  of  the  wire  and  the  relative 
energetics  of  the  wire  explosions.   Wire  resistance  increases  with  an 
increase  in  energy  content  if  the  temperature  coefficient  of  resis- 
tivity is  positive.   Melting  of  a  wire  is  characterized  by  small 
rapid  changes  in  the  resistance  curve  while  vaporization  is  indicated 
by  a  large  rapid  increase  in  the  wire  resistance  to  a  peak  value. 

The  resistance  anomaly  is  characteristic  of  wire  explosions  when 
the  deposited  energy  is  greater  than  the  energy  necessary  to  vaporize 
the  wire  (i.e.,  The  x^ire  resistance  depends  on  the  rate  of  energy  de- 
posited.)  The  resistance  reaches  larger  magnitudes  for  slower  (less 
energetic)  explosions.   Near  a  resistance  peak  a  dark  pause  is  observed 
if  the  voltage  across  the  wire  gap  does  not  cause  arc  breakdown  of  the 
surrounding  material.   The  voltage  and  current  also  have  characteristic 
responses  for  a  dark  interval;  the  current  is  approximately  zero  and 
the  voltage  remains  relatively  constant. 

The  voltage  and  current  curves  of  Figs.  8  through  33  were  examined 
for  the  characteristic  differences  between  low  heat  of  vaporization 
material  (class  I)  and  high  heat  of  vaporization  material  (class  II) .  The 
voltage  and  current  curves  for  copper  (class  I  material)  were  signifi- 
cantly different  than  the  voltage  and  current  curves  of  iron  (class  II 
material).   Typically,  class  II  materials  melt  near  the  peak  current 
(Fig.  lb).   The  current  rapidly  decreases  during  melting  and  then  re- 
mains nearly  constant  until  vaporization.   Thus,  class  II  materials  ex- 
hibit a  relatively  low  current  level  when  the  energy  to  vaporize  is 
reached.   For  class  I  materials  the  current  peaks  near  the  end  of  vapor- 
ization (Fig.  la)  and  rapidly  decreases  near  the  voltage  peak.   Thus, 
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class  I  materials  exhibit  a  high  current  level  when  the  energy  to 
vaporize  is  attained. 

Class  II  materials  require  a  larger  energy  deposit  to  vaporize 
a  wire.   Thus,  the  time  to  burst  is  longer  for  class  II  material 
while  the  voltage  and  power  peaks  are  broader.   Peaks  of  voltage  and 
power   curves  correspond  to  explosive  vaporization.   Since  class  II 
materials  require  longer  time  periods  to  deposit  the  vaporization  en- 
ergy, the  explosion  mechanism  tends  to  be  affected  by  unduloids. 
Typically  the  entire  length  of  a  class  II  material  wire  does  not  vapor- 
ize. 

The  voltage,  current,  power,  and  energy  curves  were  examined  to 
determine  the  effect  of  wire  dimensions  on  wire  explosions.   The  vol- 
tage and  current  curves  determined  if  wires  were  well  matched  to  the 
circuit  parameters.   As  discussed  in  Section  2.0,  a  wire  is  well 
matched  to  the  circuit  parameters  if  the  wire  explodes  near  the  first 
current  peak.   The  time  of  burst,  (time  at  the  voltage  peak)  which 
depends  on  wire  parameters  and  circuit  parameters,  was  examined  and 
compared  to  results  in  the  literature.   The  energy  curves  and  the  power 
curves  of  Figs.  8-33  were  compared  to  determine  if  burst  energy  and  peak 

power  per  unit  length  increased  for  decreasing  wire  length  as  reported 

23 
by  Leopold.     (Specific  values  of  the  energy  and  power  curves  obtained 

from  the  calculated  results  are  given  in  Table  II.)   Since  Leopold 
found  that  energy  deposited  up  to  one  microsecond  after  the  wire  burst 
is  effective  in  transferring  energy  to  material  surrounding  the  wire, 
the  energy  plateau  after  the  wire  burst  was  also  considered  in  the  anal- 
ysis. 

Transition  energies  of  specific  interest  in  the  analysis  of  ex- 
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ploding  wires  were  the  energies  necessary  (1)  to  melt  the  wire, 
(2)  to  vaporize  the  wire,  and  (3)  to  ionize  electrons  from  the  wire. 
Handbook  values  of  specific  heats,  latent  heats,  and  ionization  en- 
ergies  were  used  to  calculate  the  energy  densities  in  Table  III.   ' 
Assuming  thermal  equilibrium  and  energy  densities  greater  than 

the  vaporization  energy,  the  particle  energy  and  percent  ionization 

37 
can  be  calculated  with  the  Saha  equation.   Rouse   calculated  ioni- 
zation equilibrium  for  a  number  of  elements  (including  aluminum  and 
iron)  by  the  use  of  the  Saha  equation.   Rouse  considered  the  ionization 
of  any  number  of  electrons  from  the  atoms  as  a  function  of  temperature 
(up  to  10   K)  and  density  of  the  gas.   Since  energy  densities  obtained 
in  this  experimental  investigation  were  less  than  or  only  slightly 
greater  than  the  ionization  energy  of  the  first  electron,  the  Saha 
equation  was  not  used.   Instead,  approximate  results  were  calculated 
(e.g.,  Fig.  3)  assuming  that  ionization  took  place  if  kinetic  energies 
of  the  plasma  particles  was  the  same  magnitude  as  the  ionization  en- 
ergy. 

5.2   Estimate  of  Error  in  Exploding  Wire  Results 

The  errors  present  in  the  measured  discharge  circuit  parameters 
(charging  voltage,  extraneous  circuit  resistance,  and  circuit  induc- 
tance) ,  the  calibration  of  the  measuring  circuits  (voltage,  current, 
and  current  derivative) ,  and  the  measurements  made  with  the  coordinate 
comparator  are  discussed  in  Appendix  A.   Errors  present  in  the  cali- 
bration of  the  measuring  circuits  and  those  present  in  the  measure- 
ments of  the  oscilloscope  records  by  the  coordinate  comparator  deter- 
mined the  errors  in  the  final  measured  and  calculated  quantities.  (The 
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propagation  of  variance  formula  was  used  in  the  error  analysis  of 
calculated  quantities.) 

Since  the  accuracy  with  which  the  various  photographic  records 
could  be  read  determined  the  error  associated  with  the  voltage  meas- 
urements, the  error  of  the  photographic  measurements  was  taken  to  be 
the  only  source  of  error  in  the  voltage  curves.   In  general,  the 
standard  deviation  in  the  voltage  measurements  was  less  than  3%  except 
when  the  voltage  was  near  the  zeroline.    The  inductive  voltage  drop 

across  a  wire  was  corrected  using  the  relation,  L  (di/dt) ,  where   the 

w 

initial  value  of  the  wire  inductance  was  used  for  L  .   Since  the  in- 

w 

ductance  of  a  cylinder  decreases  with  an  increase  in  radius  and  a  wire 

expands  during  an  explosion,  smaller  values  of  L  should  have  been 

w 

used.   (L  was  ,  in  general,  a  function  of  time.)   The  inductive 
w 

corrections  were  generally  less  than  2%  except  near  the  beginning  of 
the  discharge,  where  the  current  varied  rapidly  (di/dt  large)  and  the 
voltage  measured  was  relatively  low.   Note  also  that  the  voltage  was 
always  near  the  peak  value  when  the  wires  were  significantly  expanded 
(i.e.,  near  or  above  the  vaporization  point).   Thus,  the  inductive 
corrections  were  much  less  than  2%  when  the  wires  were  significantly 
expanded.   As  discussed  in  Section  4.2,  the  possibility  of  induced  vol- 
tages in  the  high  voltage  probes  was  checked  and  found  to  be  small. 
Thus,  the  standard  deviation  in  the  measured  resistive  voltage  drop 
across  a  wire  (total  voltage  minus  the  induced  voltage)  was  less  than 
5%  except  near  the  initial  portion  of  a  discharge  or  where  the  voltage 
was  low. 

Errors  in  the  calibration  of  the  current  measuring  circuit  and 
in  reading  of  the  oscilloscope  current  records  determined  the  errors 
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characteristic  of  current  curves  in  Figs.  8  through  33.   The  fractional 
standard  deviation  of  the  calibration  of  the  current  measuring  circuit 
was  .033  amperes  per  ampere.   The  standard  deviation  of  distance  meas- 
ured using  the  coordinate  comparator  was  .00497  inches.  The  standard 
deviation  in  the  current  measurements  due  to  both  current  calibration 

and  coordinate  comparator  measuring  errors  was  less  than  5%  except  when 
the  current  approached  zero.  Contact  bounce  in  the  closing  of  the  air- 
gap  switch,  which  was  operated  at  lower  voltage  than  the  operational 

22 

voltage  range  given  by  Harvey  and  Yevick,   was  observed  on  the  initial 

portion  of  the  current  curves  (one  microsecond  or  less) .   Because  of 
contact  bounce,  the  current  was  not  uniquely  determined  for  the  initial 
portion  of  a  discharge.   Thus,  the  standard  deviation  in  current  meas- 
urements during  the  initial  portion  of  the  discharge  may  have  exceeded 
five  percent. 

The  energy  curves,  calculated  by  the  integration  of  the  power 
curves,  were  of  course  determined  in  part  by  the  initial  current  values. 
(Other  calculated  curves  were  not  affected  at  later  times  by  the  in- 
ability to  obtain  accurate  initial  current  measurements.)   However, 
since  the  power  curves  also  depend  on  the  voltage  curves  and  since  vol- 
tage was  low  during  the  initial  time  interval  (first  microsecond  of  a 
discharge)  in  which  current  was  not  uniquely  determined,  the  power  was 
at  a  low  level  during  the  initial  portion  of  each  discharge.   Thus,  the 
portion  of  the  power  curves  with  the  greatest  fractional  standard  de- 
viation had  small  amplitude  and  large  errors  were  not  introduced  in  the 
overall  energy  calculations. 

To  calculate  the  energy  curves  the  power  curves  were  integrated 
numerically  by  Simpson's  rule.   In  carrying  out  the  integration  process 
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for  the  energy  curves,  a  sufficient  number  of  experimental  points 
were  measured  to  keep  the  estimated  integration  error  less  than  one 

percent.   Accuracy  of  the  integration  of  the  power  curves  was  esti- 

40 
mated  by  the  relation 

E  =  -(At/90)  ["-^W^VV^l+Vl^^V---4-^ 
s 

+8(y3+y5+...+yn_3)  ]     , 

where: 

At  is  the  time  interval  between  measured  points, 

n  is  the  number  of  intervals  used  in  the  numerical  integration, 

and 
y.  is  the  instantaneous  power  after  the  i   time  interval. 

By  applying  the  propagation  of  variance  formula  to  Simpson's  rule,  it 
was  also  noted  that  the  fractional  standard  deviation  of  calculated 
energy  was  less  than  the  fractional  standard  deviation  of  the  power  at 
any  time,  if  a  sufficient  number  of  intervals  were  used  in  the  numer- 
ical integration.   Thus,  in  the  energy  calculations  a  sufficient  num- 
ber of  intervals  were  used  to  keep  the  integration  error  less  than  one 
per  cent  and  to  also  insure  that  the  fractional  standard  deviation  of 
the  calculated  energy  was  less  than  .10  joules  per  joule  at  the  time 
of  burst.   Therefore,  the  estimate  of  standard  deviation  for  the  curves 
of  the  calculated  quantities  (resistance,  power,  energy)  was  generally 
less  than  10%.   Exceptions  occurred  if  the  current  curves  or  the  vol- 
tage curves  were  near  the  zeroline,  and  during  a  short  time  interval 
when  contact  bounce  of  the  switch  was  important. 

In  determining  the  estimates  of  the  maximum  standard  deviations 
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reported  above,  standard  deviations  were  estimated  for  all  the  points 
of  reported  curves  using  the  propagation  of  variance  formula  and  the 
known  estimated  standard  deviations  of  the  measuring  circuits;  these 
maximum  estimated  standard  deviations  were  as  follows:   voltage  -  5%, 
current  -  5%,  and  calculated  values  -  10%.   To  determine  if  significant 
differences  existed  between  wire  explosions  with  identical  circuit  and 
wire  parameters,  eleven  standard  iron  wires  were  exploded  as  discussed 
in  Section  4.3.   (The  estimated  standard  deviation  of  quantities  near 
the  time  of  burst  should  have  been  significantly  less  than  the  maximum 
standard  deviation  estimated  by  the  propagation  of  variance  formula, 
since  neither  the  current  nor  the  voltage  curves  were  near  the  zeroline 
and  contact  bounce  was  also  unimportant  at  the  time  of  wire  burst.) 
Experimental  standard  deviations  noted  for  the  standard  iron  wires  were 
as  follows:   time  to  burst  -  7%,  peak  voltage  -  4%,  peak  resistance  - 
11%,  peak  power  -  12%,  and  energy  to  burst  -  13%.   Since  the  estimated 
standard  deviations  for  values  obtained  near  the  time  of  burst  of  the 
standard  wires  were  not  significantly  less  than  the  maximum  standard 
deviation  estimated  by  the  propagation  of  variance  formula,  it  was  con- 
cluded that  actual  differences  existed  in  the  wire  explosions.    The 
existence  of  such  differences  in  the  standard  wire  explosions  was  also 
indicated  by  the  estimated  standard  deviation  of  the  time  of  burst, 
since  the  sweep  rate  of  the  oscilloscope  was  more  consistent  than  the 
estimated  standard  deviation.   It  is  concluded  that  reproducibility  of 
wire  explosions  was  not  within  the  calibration  accuracy  of  the  measur- 
ing circuits  because  of  inhomogeneities  in  the  wires  and  slight  vari- 
ation in  the  degree  of  contact  between  the  wires  and  the  brass  leads 
of  the  discharge  circuit. 
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Since  only  one  size  of  wire  with  one  set  of  circuit  conditions 
was  used  in  determining  the  reproducibility  of  wire  explosions,  error 
for  other  conditions  can  only  be  discussed  qualitatively.   The  iron 
wire  results  (voltage,  current,  resistance,  power,  and  energy)  for  all 
cases  reported  (Figs.  24  -  29)  should  be  characterized  by  approximately 
the  same  fractional  standard  deviation  as  determined  for  the  standard 
wire  explosions.   However,  time  reproducibility  was  expected  to  be 

better  for  wires  exploded  near  the  beginning  of  the  capacitor  dis- 

23 
charge.     (e.g.  ,  The  results  of  Fig.  29  should  have  better  time  re- 
producibility for  the  time  of  burst  than  the  results  of  Fig.  28.) 

In  general,  the  10  mil  copper  wire  explosions  had  better  time  re- 
producibility than  the  iron  wire  explosions  since  the  copper  wire  ex- 
plosions occurred  earlier  in  the  discharge  than  the  iron  wire  explosions. 
However,  since  the  voltage  for  the  10  mil  copper  wire  explosions  in- 
creased rapidly  to  high  values,  the  response  of  the  high  voltage  probes 
could  have  limited  the  voltage  measurement  accuracy.   Thus,  the  vertical 
error  of  curves  for  the  10  mil  copper  wire  may  have  been  larger  near  the 
voltage  peaks  than  the  error  noted  for  the  standard  iron  wire. 

The  wires  which  exhibited  oscillatory  behavior  exploded  after  long 
time  intervals  (10-50  microseconds).   Thus,  inhomogeneities  in  the  wires 
had  long  time  intervals  to  affect  the  wire  explosions.   Inhomogeneities 
in  the  wires  affected  the  formation  of  unduloids  and  the  hydrodynamic 
motions  which  were  present  in  slow  wire  explosions.   Thus,  slight  initial 
inhomogeneity  in  a  wire  could  have  caused  significant  differences   in 
measured  and  calculated  quantities  for  an  exploding  wire.   (e.g.,   The 
wires  which  exploded  on  the  second  current  pulse  did  not  show  good  time 
reproducibility  for  the  time  to  burst.)   Therefore,  larger  errors   than 
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observed  for  the  standard  wire  explosions  were  expected  and  observed 
for  wires  which  exhibited  oscillatory  behavior  in  current  and  voltage 
waveforms . 

5.3  Aluminum  Wire  Results 

Figures  8  and  9  show  the  results  for  32  mil  diameter  (28 
gauge)  one  inch  long  aluminum  wires.   The  capacitance  used  in  both 
cases  was  8  microfarads,  and  the  charging  voltages  were  respectively 
11  and  15  kilovolts.   Aluminum  wires  did  not  explode  for  capacitor 
discharges  at  11  kilovolts  but  remained  clamped  in  the  circuit  after 
capacitor  bank  discharge,  although  they  were  weakened  at  the  point  of 
attachment  to  the  massive  brass  leads.   The  wires'  physical  appear- 
ance indicated  surface  melting  and  possibly  some  surface  oxidation  had 
occurred.   The  voltage  and  current  traces  had  the  characteristics  of 
an  RLC  underdamped  circuit  for  the  entire  discharge. 

When  the  capacitor  bank  was  initially  charged  to  15  kilovolts 
and  then  discharged  through  an  aluminum  wire,  the  voltage  peaked  dur- 
ing the  second  current  pulse.   Sufficient  energy  to  vaporize  the  wires 
was  deposited  before  the  time  of  burst.   However,  Table  II  shows  that 
deposited  energy  densities  were  not  sufficient  to  cause  appreciable 
ionization,  since  the  energy  density  at  which  electrons  are  ionized 
from  aluminum  atoms  is  9.74  eV/atom  (See  Table  III).   Since  the  vol- 
tage peaks  were  not  relatively  high  and  the  voltage  increases  were  not 
rapid,  the  aluminum  wire  explosions  were  classified  as  slow  explosions, 
Data  were  not  taken  at  4  microfarads  capacitance  for  aluminum  wires 
since  low  energy  oscillatory  responses  do  not  generate  the  high  temper- 
atures and  high  energy  densities  necessary  to  simulate  a  nuclear  ex- 
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plosion. 

5.4   Copper  Wire  Results 

Figures  10,  11,  and  12  show  the  results  for  20  mil  diameter 
copper  wires  exploded  at  8  microfarads  and  11  kilovolts.   The  voltage 
and  current  curves  resemble  the  oscillatory  response  of  an  underdamped 
RLC  circuit  discharge.   Sufficient  energy  to  vaporize  the  wire  was  not 
deposited  until  after  the  fourth  current  pulse.   Thus,  physical  dis- 
tortions dominated  the  explosion  mechanism;  the  wire  melted  and  fell 
apart  like  a  fuse.   Arc  breakdown  of  the  surrounding  air  allowed  the 
discharge  to  continue  the  oscillatory  response.   As  in  the  case   of 
32  mil  aluminum  wire,  data  were  not  taken  at  4  microfarads  capacitance 
for  20  mil  copper  wire. 

Figures  13  through  16  show  the  results  for  20  mil  diameter 
copper  wires  exploded  at  8  microfarads  and  15  kilovolts.   These  curves 
are  particularly  interesting  since  the  explosion  occurs  during  the 
second  current  pulse.   Most  of  the  reported  exploding  wire  results  are 
for  cases  of  explosion  during  the  first  current  pulse  or  slow  discharges 
with  oscillatory  current  and  voltage  curves. 

The  voltage  peak  (time  of  burst)  was  not  completely  reproducible 
when  20  mil  copper  wires  were  exploded  using  the  same  wire  and  circuit 
parameters.   Since  literature  results  indicate  that  best  time  reproduc- 
ibility is  noted  for  wires  exploding  near  the  beginning  of  the  first 

23 
current  pulse,     good  time  reproducibility  was  not  expected  for  wires 

exploding  on  the  second  current  pulse.   Factors  which  affected  the 

time  to  burst  were  the  wire-circuit  contact  and  inhomogeneities  of  the 

wires.   For  example,  small  initial  differences  in  a  wire  would   be 
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greatly  increased  by  physical  distortions  over  extended  discharge  in- 
tervals (10-15  microseconds  to  the  second  current  pulse) . 

Figure  34  shows  the  peak  power  as  a  function  of  wire  length  for 
the  20  mil  copper  wires  exploded  at  8  microfarads  and  15  kilovolts. 
The  peak  power  was  approximately  constant  for  the  wire  lengths  of  1", 
2",  and  3"  but  definitely  lower  for  the  1/2"  wire.   The  peak  power 
per  unit  length  increased  with  decreasing  wire  length.   Since  the 
measured  peak  power  per  unit  length  of  the  1/2"  wire  was  not  as  large 
as  predicted  by  extrapolating  the  results  for  the  other  three  wire 
lengths  to  1/2  inch,  end  effects  were  particularly  important  for  the 
1/2"  long  20  mil  diameter  copper  wire. 

The  results  shown  in  Fig.  34  indicate  that  the  peak  power  per 
unit  length  of  wire  for  wires  which  explode  on  the  second  current 
pulse  depends  on  wire  length  in  the  same  way  as  wires  which  explode 
on  the  first  current  pulse.   The  increase  in  peak  power  per  unit  wire 

length  with  a  decrease  in  wire  length  observed  for  the  20  mil  copper 

23 

wire  agrees  with  the  results  of  Leopold.    Leopolds 's  results,  however, 

are  for  wire  explosions  during  the  first  current  pulse  and  different 
circuit  and  wire  parameters  than  used  in  this  work.   The  peak  power 
curve  of  Leopold  shows  a  maximum,  but  only  a  broad  maximum  was  indi- 
cated in  results  reported  here.   Figure  34  indicates  that  the  maximum 
occurs  for  wires  longer  than  1/2"  but  insufficient  data  were  available 
to  determine  at  exactly  what  wire  length  the  maximum  occurred. 

The  energy  deposited  in  the  various  wire  lengths  are  compared  to 
the  amount  of  energy  necessary  to  vaporize  the  wire  in  Fig.  35  with 
the  time  interval  after  burst  used  as  a  parameter.   The  amount  of  energy 
above  the  vaporization  energy  was  greatest  for  the  1/2"  wire  for  all 
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Figure  34.   Peak  power  as  a  function  of  wire  length. 
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Figure  35.   Energy  deposition  profiles  for  20  mil  diameter  copper 
wire  of  various  length. 
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time  intervals  of  interest.   The  amount  of  energy  above  the  vapori- 
zation energy  was  about  the  same  for  the  1"  and  2"  wires  for  the  time 
intervals  of  interest.   The  vaporization  energy  was  not  deposited  in 
the  3"  wire  before  the  end  of  the  second  current  pulse. 

Figure  36  shows  the  energy  density  deposited  in  a  20  mil  copper 
wire  as  a  function  of  wire  length  with  the  time  after  burst  used  as  a 
parameter.   Here  the  energy  density  increased  with  decreasing  wire 
length.   The  more  time  considered  after  burst  the  sharper  the  increase 
in  energy  density  with  decrease  in  wire  length.   The  results,  however, 
do  not  determine  if  the  energy  deposited  immediately  after  burst  raised 
the  temperature  of  the  wire  higher  than  the  temperature  at  burst. 
Leopold's  results  showed  that  the  energy  deposited  in  a  short  time  in- 
terval (approximately  one  microsecond)  after  the  time  of  burst  is  im- 
portant in  obtaining  high  energy  densities  or  at  least  keeping  the  en- 

23 

ergy  densities  at  high  levels.     Thus,  explosions  of  shorter  wires 

have  high  energy  densities  for  longer  periods  of  time  than  explosions 
of  longer  wires  if  end  effects  are  not  important. 

The  curves  of  Figs.  17  through  23,  showing  the  results  obtained 

for  10  mil  diameter  copper  wire,  have  the  same  shape  as  those  reported 

25 

by  Bennett.     The  voltage  traces  exhibit  the  step  change  in  voltage 

which  corresponds  to  the  melting  of  the  wire.   The  magnitude  of  the 
voltage  peaks  were  about  twice  the  charging  voltage  of  the  capacitor 
bank;  the  peak  power  was  also  significantly  greater  than  the  theoretical 
estimate  for  an  RLC  discharge  with  constant  resistance.   These  voltage 

and  peak  power  characteristics  are  also  observed  in  results  reported  by 

25 

Bennett. 

As  discussed  in  Section  2.2,  the  temperature  coefficient  of  wire 
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Figure  36.   Deposited  energy  density  for  20  mil  diameter  copper  wire 
of  various  length.   Time  after  burst  used  as  parameter. 
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resistance  determines  the  shape  of  the  resistance  curve  between  melt- 
ing and  vaporization.   The  temperature  coefficient  of  resistivity  for 
the  liquid  phase  of  copper  is  positive  and  thus,  the  resistance  curves 
for  the  10  mil  copper  wire  explosions  increase  between  the  melting 
and  the  beginning  of  vaporization.   Comparison  of  the  resistance  curves 
for  10  mil  copper  wire  explosions  (Figs.  17  through  23)  show  that  the 
curves  corresponding  to  the  more  energetic  explosions  had  lower  re- 
sistance peaks;  thus  the  resistance  anomaly  was  observed.   The  dark 
pause  which  is  related  to  the  resistance  anomaly  (discussed  in  Section 
2.2)  was  observed  f jr  two  of  the  less  energetic  wire  explosions  (i.e., 
the  3  inch  wire  which  exploded  at  8  microfarads  and  11  kilovolts  and 
the  1  inch  wire  which  exploded  at  8  microfarads  and  8  kilovolts) . 

The  time  to  burst,  which  was  taken  to  be  the  time  interval  from 
the  initiation  of  discharge  to  the  voltage  peak,  decreased  with  a 
decrease  in  wire  length  while  an  increase  in  the  initial  charging  vol- 
tage decreased  the  time  to  burst.   The  data  taken  for  the  10  mil  copper 
wire  did  not  show  a  significant  dependence  of  the  time  to  burst  on 
capacitance.   The  explosion  of  the  10  mil  copper  wire  during  the  first 

current  pulse  and  the  explosion  of  the  20  mil  copper  wire  during  the 

23  24 
second  current  pulse  were  consistent  with  reported  results.       The 

literature  results  in  general  show  that  the  time  to  burst  increases 

with  an  increase  in  wire  diameter  if  other  parameters  are  held  constant. 

14 
Since  copper  has  a  low  heat  of  vaporization   it  would  be  ex- 
pected that  the  time  to  deposit  sufficient  energy  to  vaporize  the  copper 
is  short  compared  to  the  time  necessary  for  instabilities  to  develop. 
The  current  curves  of  Figs.  17  through  23  have  maximum  values  near  the 
end   of  vaporization  which  agrees  with  the  characteristics  of  a  class  I 
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14 
material  described  by  Webb.     The  10  mil  copper  wire  results  in 

23 

Table  II  are  also  consistent  with  results  reported  by  Leopold, 

i.e.,  (1)  peak  power  and  deposited  energy  increase  with  an  increase 
in  circuit  capacitance  or  circuit  charging  voltage,  and  (2)  energy 
density  and  peak  power  per  unit  length  increase  with  decreasing  wire 
length. 

The  energy  density  necessary  to  completely  vaporize  a  one  inch 
copper  wire  10  mils  in  diameter  is  4.13  eV/atom  (See  Table  III).   In 
contrast  with  the  iron  wire  results  sufficient  energy  was  deposited 
to  vaporize  completely  the  copper  wires  before  the  time  to  burst,  ex- 
cept for  the  three  inch  wire  exploded  at  8  microfarads   and  11  kilo- 
volts.   However,  within  the  reliability  of  the  handbook  data  used  in  the 
the  theoretical  calculations  and  the  experimental  standard  deviation, 
the  vaporization  energy  for  the  3  inch  copper  wire  exploded  at  8  micro- 
farads and  15  kilovolts  was  deposited  before  the  energy  plateau  was 
obtained. 

Although,  in  general,  sufficient  energy  was  deposited  to  vapor- 
ize completely  the  copper  wires  before  the  time  of  burst,  unduloid 
characteristics  were  observed  in  the  explosions.   The  voltage  curves 
for  the  3  inch  wires  show  vaporization  peaks  even  though,  as  noted 
above,  sufficient  energy  was  not  deposited  to  vaporize  entirely  the 
wire  material  at  the  time  the  voltage  peaks  were  noted.   Unduloids  were 
more  prominent  in  the  explosive  mechanism  for  longer  wires,  since  the 
burst  times  were  longer.   (These  unduloids  were  probably  formed  in 
regions  of  the  wire  with  initially  higher  resistivity  where  more  energy 
was  deposited  and  where  the  positive  coefficient  of  resistivity  of  cop- 
per further  increased  inhomogeneity  to  form  the  unduloids.   Thus,   the 
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specific  energy  deposited  along  the  length  of  v/ire  was  not  constant  and 
vaporization  peaks  in  the  voltage  curves  were  observed  even  though  the 
total  energy  deposited  was  not  sufficient  to  vaporize  completely  the 
wire. 

If  the  energy  deposited  in  the  10  mil  copper  wires  were  greater 
than  necessary  to  vaporize  entirely  the  wire,  it  would  be  expected  that 
the  remaining  energy  would  be  used  to  heat  the  wire  material  above  the 
vaporization  temperature.   If  the  energy  density  were  raised  higher  than 
11.85  eV/atom,  electrons  would  be  ionized  from  the  copper  atoms.   At  the 
time  of  burst  the  *  ighest  energy  density  obtained  in  the  10  mil  diameter 
copper  wire  explosions  in  this  study  (i.e.,  for  a  one  inch  wire  exploded 
at  8  microfarads  and  15  kilovolts)  was  8.68  eV/atom. 

Although  the  energy  densities  deposited  at  the  time  of  burst  were 
less  than  the  ionization  energy  (The  energy  density  necessary  for 
ionization  of  copper  atoms  is  approximately  11.85  eV/atom.),  LeopoldTs 
results  show  that  the  energy  deposited  in  a  short  time  interval  after 
the  time  of  burst  is  important  in  obtaining  high  energy  densities.   Thus, 
deposited  energy  densities  in  the  first  current  pulse  (last  column  of 
Table  II)  were  considered  to  determine  if  electron  ionization  was  achiev- 
ed.  Considering  the  first  pulse  energies  of  the  two  most  energetic   10 
mil  copper  wire  explosions  (12.70  eV/atom  and  11.60  eV/atom) ,  it  is  seen 
that  some  electron  ionization  \>7as  probably  achieved.   However,  even  for 
these  energetic  cases,  the  energy  densities  were  not  great  enough   to 
completely  ionize  the  first  electron  from  all  copper  atoms  in  the  wire. 
Therefore,  the  temperature  needed  for  the  ionization  of  the  first  elec- 
tron from  the  copper  atoms  was  an  upper  limit  of  the  temperatures  reach- 
ed by  copper  wires  exploded  during  this  study. 


88 


5.5   Iron  Wire  Results 

The  results  for  iron  wires  exploded  at  atmospheric  conditions 
are  shown  in  Figs.  24  through  29;  the  characteristics  of  the  resistance 
curves  were  determined  by  the  physical  condition  of  the  wire  at  various 
times  during  the  discharge.   The  small  jump  in  resistance  noted  between 
one  and  two  microseconds  after  the  initiation  of  the  discharge  was 
caused  by  melting  of  the  wire.   Since  the  temperature  coefficient   of 
resistance  for  the  liquid  phase  of  iron  is  small,  the  resistance  curves 

have  a  nearly  constant  value  between  the  end  of  melting  and  the  begin- 

28 
ning  of  vaporization.     The  resistance  reached  larger  magnitudes  for 

the  less  energetic  explosions,  a  characteristic  of  resistance  anomaly 

of  exploding  wires.   From  these  curves  it  is  noted  that  when  energy 

available  from  the  capacitor  bank  was  decreased,  or  when  the  length  of 

the  wire  was  increased,  the  restrike  occurred  later  in  the  explosion. 

It  is  also  apparent  that  arc  breakdown  of  the  air  surrounding  the  wire 

occurred  earlier  in  the  explosion  for  the  more  energetic  cases. 

The  time  to  burst  (time  interval  from  the  initiation  of  the  dis- 
charge to  the  voltage  peak)  decreased  with  a  decrease  in  wire  length; 
an  increase  in  the  initial  charging  voltage  also  decreased  the  time  to 
burst.   The  data  taken  for  the  iron  wire  did  not  show  a  significant 
dependence  of  the  time  to  burst  on  capacitance.   The  time  to  burst  of 
the  14  mil  iron  wire  was  greater  than  the  10  mil  copper  wire  with  other 
parameters  constant;  this  was  expected  since  iron  has  a  higher  boiling 
point  than  copper  and  also  the  iron  wires  were  of  larger  diameter  than 
the  copper  wires. 

A  comparison  of  results  noted  for  one  inch  and  three  inch  iron 
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wires  exploded  at  8  microfarads  and  11  kilovolts  were  consistent  with 

23 
the  results  of  Leopold;    the  energy  density  and  the  peak  power  per 

unit  length  were  greater  for  the  shorter  wires.   Iron,  with  its  large 

14 
heat  of  vaporization  is  classified  as  a  class  II  material  by  Webb. 

The  current  curves  in  Figs.  24  through  29  have  peaks  near  the  time  of 

melting  (i.e.,  near  the  resistance  jump)  and  the  current  curves  are 

at  low  values  or  rapidly  decreasing  at  the  time  of  burst.   Thus  the 

current  curves  for  the  iron  wire  explosions  are  characteristic  of 

those  for  a  class  II  material. 

The  energy  to  raise  an  iron  wire  from  room  temperature  to  the 
atmospheric  boiling  point  and  completely  vaporize  the  iron  was  calcu- 
lated from  heat  capacity  and  latent  heat  values  (See  Table  III) .   The 
results  of  the  iron  wire  explosions  show  that  except  for  the  more 
energetic   explosions,  sufficient  energy  to  vaporize  the  wires  com- 
pletely was  not  deposited  before  burst.   (Compare  values  in  column  9, 
Table  II,  with  values  in  column  7,  Table  III.)   Even  for  the  most 
energetic  cases,  corresponding  to  Figs.  25  and  29,  the  deposited  energy 
before  burst  was  not  much  greater  than  the  vaporization  energy. 
Sufficient  energy  to  vaporize  the  wire  for  the  case  corresponding  to 
Fig.  27  was  deposited  during  the  first  current  pulse  but  after  the  time 
of  burst.   In  remaining  iron  wire  experiments,  sufficient  energy  to 
vaporize  the  wire  completely  was  not  deposited. 

The  time  to  deposit  sufficient  energy  for  vaporization  of  class 
II  materials  is  typically  long  compared  to  the  instability  time  constant, 
Thus,  the  mechanism  of  the  iron  wire  explosions  was  controlled  by  physi- 
cal distortions  of  the  wire  (unduloids),  and  only  segments  of  the  wire 
were  vaporized.   Fewer  unduloids  were  formed  and  larger  portions  of  the 
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wire  were  vaporized  in  the  more  energetic  wire  explosions  with  higher 
peak  powers. 

The  ionization  energy  at  which  the  first  electron  is  ionized 
from  an  iron  atom  is  12.70  eV/atom.   Table  II  shows  that  energy  den- 
sities obtained  experimentally  in  this  investigation  were  not  sufficient 
to  ionize  electrons  from  iron  wires.   The  results  in  Table  II  also  show 
that  peak  power  and  deposited  energy  increase  with  an  increase  in  cir- 
cuit capacitance  or  circuit  charging  voltage.   Thus,  an  increase  in 
charging  voltage  or  circuit  capacitance  would  increase  the  deposited 
energy  and  produce  higher  temperatures.   However,  all  eight  available 
capacitors  were  used  to  obtain  data  presented  in  Table  II  and  an  in- 
crease in  voltage  was  not  practical  during  the  present  work  because  the 
high  voltage  probes  were  already  being  used  at  their  upper  operating 
limit. 

5.6  Nickel-Chromium  Wire  Results 

Figure  30  shows  the  results  for  a  28.5  mil  (21  gauge)  one  inch 
long  nickel-chromium  wire  exploded  at  8  microfarads  and  11  kilovolts. 
Hot  metal  was  scattered  from  the  original  wire  position  in  the  discharge 
circuit  indicating  that  the  wire  simply  melted  and  was  blown  apart.  The 
voltage  and  current  were  oscillatory. 

Figure  31  shows  the  results  obtained  when  a  nickel-chromium  wire 
was  exploded  at  8  microfarads  and  15  kilovolts.   Sufficient  energy  was 
deposited  during  the  first  pulse  of  this  explosion  to  reach  the  boiling 
point  of  chromium  (2495  K)  and  to  vaporize  the  chromium.   However,  the 
power  was  low  when  the  vaporization  energy  had  been  deposited  and  heat 
losses  kept  the  temperature  below  the  boiling  point  of  chromium  until 
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the  second  current  pulse.   The  voltage  peak  occurring  on  the  second 
current  pulse  exhibited  a  rapid  rise  but  did  not  have  large  magnitude. 
At  the  time  of  burst  sufficient  energy  had  been  deposited  to  reach  the 
boiling  point  of  iron  (3008  K)  or  nickel  (3010  K) ,   but  sufficient 
energy  had  not  been  deposited  to  vaporize  the  wire  material  completely. 
The  vaporization  of  either  chromium  alone  or  only  segments  of  the  wire 
completely  vaporizing  caused  the  voltage  spike.   Since  long   times 
(~  15  microseconds)  were  available  for  physical  distortions  to  affect 
the  explosion,  the  most  probable  mechanism  was  complete  vaporization  of 
segments  of  the  wire. 

As  noted,  sufficient  energy  was  apparently  deposited  before  the 
end  of  the  second  current  pulse  to  vaporize  the  nickel-chromium  wire 
completely.   However,  the  wire  did  not  remain  together  more  than  about 
a  microsecond  after  the  voltage  peak  which  was  early  in  the  second 
pulse.   Thus,  the  explosion  was  dominated  by  arc  discharge  near  the  end 
of  the  second  current  pulse  and  the  wire  did  not  actually  absorb  energy 
during  the  last  half  of  the  second  current  pulse. 

5.7   Summary  of  Bare  Wire  Results 

The  characteristics  of  the  voltage,  current,  resistance,  power, 
and  energy  curves  of  Figs.  8  through  31  are  examined  in  Sections  5.3 
through  5.6.   These  curves  show  good  qualitative  agreement  with  the 
literature;  quantitative  agreement  could  not  be  ascertained,  however, 
since  the  circuit  parameters  were  generally  different  from  those  de- 
scribed in  the  exploding  wire  literature.   Significant  variation  between 
separate  wire  explosions  was  noted  for  wire  explosions  with  identical 
circuit  and  wire  parameters  (See  Section  5.2).   It  is  expected  that  the 
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experimental  error  in  future  experiments  can  be  reduced  if  a  more 
consistent  method  is  devised  for  attaching  the  wire  to  the  brass 
leads  and  if  additional  improvements  in  the  triggering  circuit  are 
made. 

The  characteristics  of  the  curves  shown  in  Figs.  8  through 
31  are  consistent  with  the  physical  properties  of  the  wires.   Melting 
of  a  wire  was  characterized  by  a  small  jump  in  the  resistance  curve 
while  the  temperature  coefficient  of  resistivity  was  in  agreement  with 
the  behavior  of  the  resistance  curves  in  the  time  interval  between 
phase  changes.   The  voltage  curves  had  vaporization  spikes  for  cases 
which  had  deposited  energy  densities  less  than  the  vaporization  energy 
at  the  time  of  burst;  this  is  characteristic  of  a  typical  unduloid 
mechanism.   The  physical  characteristics  of  low  heat  of  vaporization 
and  high  heat  of  vaporization  materials  were  observed  in  the  10  mil 
copper  wire  and  the  14  mil  iron  wire  results  respectively.   The  re- 
sistance anomaly,  which  is  dependent  on  the  physical  condition  of  the 
wire,  existed  for  wire  explosions  which  had  sufficient  energy  deposited 
for  vaporization. 

The  characteristics  of  the  explosions  noted  in  this  study  may  be 
summarized  as  follows:   (1)  oscillatory  voltage  and  current  response, 
(2)  explosion  of  the  wire  on  the  second  current  pulse,  and  (3)  explosion 
of  the  wire  on  the  first  current  pulse.   Oscillatory  responses  are  known 
to  be  characteristic  of  exploding  wires  x^ith  slow  energy  deposit;  they 
produce  low  temperatures  and  are  thus  not  of  interest  in  simulating  a 
nuclear  explosion.   The  wire  explosions  noted  during  the  second  current 
pulse  were  unique  to  this  study;  literature  results  showed  either 
oscillatory  responses  or  responses  with  explosion  during  the  first  current 
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pulse  only. 

The  energetics  of  this  study's  wire  explosions  were  consistent 
with  exploding  wire  literature,  i.e.,  energy  density  increased  with 
an  increase  in  charging  voltage  and  a  decrease  in  wire  length.   (End 
effects  were  important  for  wires  shorter  than  about  one-half  inch.) 
The  relation  between  copper  wire  diameter  and  time  to  burst  was  con- 
sistent with  literature  results,  i.e.,  the  time  to  burst  decreased 
with  a  decrease  in  wire  diameter.   Copper,  a  low  heat  of  vaporization 
material,  exhibited  more  energetic  wire  explosions  than  iron,  a  high 
heat  of  vaporization  material.   The  copper  wires  which  exploded  near 
the  peak  of  a  current  pulse  were  well  matched  to  the  circuit  parameters 
and  produced  the  more  energetic  wire  explosions. 

If  an  appreciable  amount  of  material  is  used  to  simulate  the 
chemical  environment  of  an  underground  nuclear  explosion,  electrons 
must  be  ionized  not  only  from  the  wire  material  but  also  from  the 
material  surrounding  the  wire.   To  obtain  appreciable  ionization  of 
material   surrounding  the  wire,  the  wire  material  must  have  much  greater 
energy  densities  initially  than  required  for  the  surrounding  material. 
Tables  II  and  III  show  that  the  energy  deposited  in  the  wires  before 
burst  was  not  sufficient  to  ionize  electrons  from  the  wire  material. 

Although  energy  deposited  in  a  short  time  interval  (approximately  one 

23 
microsecond)  after  the  time  of  burst  can  actually  be  added  to  the  system, 

energy  densities  obtained  during  these  experiments  were  still  not 
sufficient  to  ionize  material  surrounding  the  wires.   Thus,  energy  den- 
sities greater  than  obtained  in  the  bare  wire  explosions  of  this  study 
must  be  achieved.   Increased  energy  deposit  in  the  wires  can  be  attained 
by  increasing  the  charging  voltage  or  capacitance;  as  an  alternative, 
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containment  of  wire  material  for   a   longer   period   of   time 
would  also  result  in  increased  energy  deposit. 

The  upper  temperature  limit  attained  by  wires  exploded  during 
this  study  was  estimated  by  the  two  methods  developed  in  Section  3.0. 
The  two  methods  assumed  respectively  (1)  no  energy  losses  from  the  wire 
before  the  time  of  wire  burst,  and  (2)  radiation  losses  were  the  only 
significant  energy  losses  during  the  discharge.   (The  assumption  of 
the  second  method  results  in  radiation  losses  balancing  the  power  in- 
put to  the  wire  at  the  time  of  the  peak  power.)   The  results  calculated 
by  the  two  methods  are  reported  in  Table  IV  and  show  that  the  maximum 
temperatures  estimated  by  the  second  method  are  an  order  of  magnitude 
larger  than  the  maximum  temperatures  estimated  by  the  first  method. 
From  this  observation  it  appears  that  a  discrepancy  has  resulted  since 
the  second  method  considered  energy  losses  while  the  first  method  did 
not.   However,  the  cases  reported  in  Table  IV  have  power  curves  which 
increase  rapidly  to  their  peak  values;  this  means  that  the  power  (rate 
of  energy  deposit)  was  more  than  sufficient  to  balance  the  radiant 
energy  losses  at  the  time  of  the  peak  power.   Therefore  it  is  concluded 
that  enough  energy  was  not  actually  deposited  in  the  wire  before  the 
peak  power  to  account  for  the  maximum  temperatures  estimated  by  the 
second  method. 

The  time  of  burst  of  an  exploding  wire  is  near  the  peak  power 
(i.e.,  at  the  peak  voltage);  a  wire  is  therefore  not  confined  after  the 
peak  power.   If  it  were  possible  to  keep  wire  intact  a  longer  time, 
sufficient  energy  would  be  deposited  to  produce  temperatures  that  better 
correspond  to  an  equality  between  radiation  losses  and  peak  power. 
Therefore  it  is   desirable  to  develop  a  technique  to  contain  the  wire 
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current  pulse;  the  current  was  initially  lower  for  the  sample  cell 
geometry  than  for  the  bare  iron  wire  explosion.   This  result  would  be 
expected  if  adding  the  sample  cell  arrangement  to  the  circuit  increased 
the  circuit  resistance.   (Increased  circuit  resistance  would  lower  the 
current  for  the  initial  portion  of  the  discharge.)   Note  also  that  the 
current  during  the  initial  portion  of  the  discharge  was  not  accurately 
determined  because  of  experimental  difficulties  (See  Section  5.2). 

The  current  noted  for  both  the  bare  wire  and  the  wire  in  the 
sample  cell  were  nearly  the  same  magnitude  shortly  after  the  wire 
melted  but  before  the  vaporization  voltage  peak.   In  this  region  the 
observed  currents  for  the  two  cases  were  within  the  variation  observed 
for  the  currents  of  the  eleven  standard  one  inch  iron  wire  explosions. 
Since  the  wire  resistance  dominated  the  circuit  resistance  during  the 
time  interval  after  the  wire  melted  (for  both  bare  wire  explosions  and 
sample  cell  explosions) ,  current  and  voltage  were  not  expected  to  be 
affected  by  the  sample  cell  during  the  time  interval  after  the  wire 
melted  but  before  the  vaporization  voltage  peak. 

The  current  rapidly  decreased  to  low  values  immediately  after 
the  time  of  burst  for  an  iron  wire  in  the  cell  geometry;  thus  increased 
wire  resistance  was  observed  in  the  cell  geometry.   This  increase  in 

wire  resistance  in  the  cell  geometry  was  expected  based  on  the  work  of 

29 

Reithel  and  Blackburn  with  confined  wires. 

The  explosion  of  a  14  mil  iron  wire  one  inch  long  at  8  microfarads 
and  11  kilovolts  was  not  as  rapid  in  the  cell  arrangement  as  under  ambient 
atmospheric  conditions;  the  time  to  burst  in  the  cell  geometry  was  3.42 
microseconds  while  the  time  to  burst  for  a  bare  iron  wire  was  3.04  micro- 
seconds.   The  increased  time  for  the  iron  wire  to  explode  in  the  sample 
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cell  arrangement  agrees  with  the  results  of  Nash  et.  al.,  who  reported 
that  pulses  were  broadened  and  flattened  for  confined  wire  explosions. 

The  energy  deposited  in  the  iron  wire  in  the  sample  cell  geometry 
was  not  significantly  different  from  the  deposited  energy  in  the  bare 
iron  wire  experiments.   The  average  energy  deposited  at  the  time  of 
burst   for  11  bare  wire  explosions  was  119  joules  with  an  estimated 
standard  deviation  of  15  joules.   The  energy  deposited  in  the  iron  wire 
in  the  cell  geometry  was  100  joules,  which  is  lower  than  the  average  of 
the  bare  wire  explosions  by  about  one  standard  deviation. 

The  peak  power  was  definitely  lower  for  the  sample  cell  geometry 
than  for  the  bare  wire  geometry.   The  average  of  the  11  standard  one  inch 
wires  was  91  megawatts  with  an  estimated  standard  deviation  of  11  mega- 
watts.  The  peak  power  for  the  sample  cell  geometry,  49  megawatts,  was 
lower  because  the  voltage  pulse  of  the  confined  wire  was  not  as  high  as 
those  for  the  bare  wire.   The  voltage  peak  in  the  cell  geometry  was 
apparently  limited  by  sparkover  from  the  wire  to  the  cell  wall  during 
the  time  interval  of  the  peak  voltage  drop  across  the  wire.   Although  the 
peak  power  was  lower  for  the  sample  cell  geometry,  the  amount  of  energy 
deposited  was  not  significantly  different.   Thus,  if  sparkover  had  not 
limited  the  voltage  peak,  higher  energy  densities  would  have  been  deposit- 
ed for  the  sample  cell  geometry  than  for  bare  wire  explosions.   It  is 
therefore  evident  that  more  attention  must  be  given  to  sample  cell  design 
with  improved  short  circuit  prevention. 

One  inch  long,  20  mil  diameter  copper  wires  were  also  exploded  at 
8  microfarads  and  11  kilovolts  in  the  sample  cell.  Comparing  the  sample 
cell  results  of  Fig.  33  with  the  bare  20  mil  copper  wire  results  of  Fig. 
11,  it  is  noted  that  the  explosion  in  the  sample  cell  showed  a  definite 
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voltage  peak  during  the  second  current  pulse  while  the  peak  was  not 
observed  for  the  bare  wire  explosions.   The  sample  cell's  voltage 
peak  indicated  that  the  sample  cell  stabilized  and  confined  the  wire 
after  melting.   This  conclusion  is  based  on  the  fact  that  the  period 
of  the  waveforms  was  shorter  for  confined  wires  than  for  bare  wires. 
(In  fact,  the  bare  wire  period  was  nearly  the  same  as  the  period  of 
the  free  response  with  a  shorting  bar  used  in  place  of  the  wire.)   The 
shorter  period  of  the  sample  cell  wire  explosion  indicated  that  the 
coaxial  feature  of  the  sample  cell  geometry  had  reduced  the  inductance  of 
the  circuit.   The  second  pulse  peak  power  levels  were  about  5  times 
larger  using  the  sample  cell  geometry  as  compared  with  bare  wire  results 
for  the  same  voltage  and  capacitance.   The  deposited  energy  was  also 
greater  for  the  sample  cell  geometry.   The  calculated  energy  required  to 
vaporize  the  20  mil  copper  wire  completely  at  the  atmospheric  boiling 
point  was  282  joules.   Since  the  sample  cell  voltage  peak  was  near  the 
end  of  the  second  current  pulse,  most  of  the  second  current  pulse  en- 
ergy (226  joules)  was  deposited  in  the  wire.   Thus,  the  deposited  ener- 
gy was  less  than  the  energy  to  completely  vaporize  the  wire.   Since  the 
voltage  peak  was  broad  and  flat,  the  explosion  in  the  cell  geometry  was 
not  expected  to  completely  vaporize  the  20  mil  copper  wire. 

The  voltage  peak  for  the  20  mil  copper  wire  was  less  than  3 
kilovolts  in  the  cell  geometry.   The  peak  voltage  for  the  iron  wire  in 
the  cell  geometry  was  7.7  kilovolts,  while  the  voltage  peak  for  bare 
iron  wire  was  about  13  kilovolts.   Thus  voltage  drop  across  the  iron  wire 
in  the  cell  geometry  was  probably  limited  by  sparkover.   Consistent  with 
this  interpretation  was  the  observation  that  when  attempting  to  use 
sample  cells  with  diameters  less  than  3/8"   it  was  not  possible  to  obtain 
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meaningful  results  because  of  voltage  sparkover.   The  10  mil  diameter 
copper  wire  exploded  at  8  microfarads  and  11  kilovolts  in  air  had  a 
voltage  peak  of  approximately  27  kilovolts  (See  Fig.  20).   A  voltage 
trace  taken  for  a  10  mil  diameter  copper  wire  exploded  at  8  micro- 
farads and  11  kilovolts  in  the  sample  cell  geometry  had  the  character- 
istic jump  associated  with  the  wire  melting,  but  the  voltage  suddenly 
dropped  to  zero  a  short  time  after  the  wire  melted.   Sparkover  from 
the  wire  to  the  tube  wall  would  account  for  the  sudden  voltage  drop. 
Thus,  although  the  sample  cell  geometry  stabilized  and  confined  wires 
after  melting,  breakdown  of  the  insulation  between  the  wire  and  the 
tube  wall  caused  sparkover  to  the  cell  wall,  limiting  the  voltage  peak 
and  precluding  the  attainment  of  more  energetic  explosions.   Therefore, 
future  work  should  consider  the  design  of  an  improved  sample  cell. 
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6.0   Conclusions 

The  basic  considerations  in  the  simulation  of  the  chemical  en- 
vironment of  a  nuclear  explosion  are  the  generation  of  high  tempera- 
tures, high  pressures,  strong  shock  waves,  nonequilibrium  conditions, 
and  a  material  environment  of  the  correct  composition.   Knowing  the 
material  environment  of  a  nuclear  explosion  and  the  necessary  tempera- 
ture distribution  and/or  energy  densities  about  30  seconds  after  the 
shot  (the  approximate  time  at  which  the  chemistry  of  the  fallout  is 
determined)  the  required  exploding  wire  circuit  parameters  to  simulate 
a  nuclear  environment  can  be  calculated  theoretically  as  discussed  in 
Section  3.0.   However,  for  the  theoretically  calculated  circuit 
parameters  to  produce  the  desired  temperatures,  wires  must  be  well 
matched  to  the  circuit  parameters. 

Since  radiation  energy  losses  are  large  at  high  temperatures, 
energy  must  be  deposited  rapidly  in  an  exploding  wire  to  simulate  a 
nuclear  explosion.   Thus,  it  would  seem  that  peak  power  would  be  the 
important  factor  in  describing  the  energetics  of  a  wire  explosion. 
However,  the  experimental  results  for  wire  explosions  during  the  second 
current  pulse  (See  Section  5.4)  show  that  these  wire  explosions  had 
high  energy  densities  and  thus  give  experimental  evidence  that  extremely 
fast  explosions  with  extremely  high  peak  powers  may  not  be  absolutely 
necessary  to  obtain  temperatures  above  10    K.     It  was  also  concluded 
from  these  experimental  results  that  parameters  other  than  peak  power 
are  important  in  obtaining  the  greatest  energy  deposit  and  highest  tem- 
peratures . 

Since  enough  energy  was  not  deposited  before  the  peak  power  to 
account  for  the  maximum  temperatures  estimated  by  the  method  of  balancing 
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radiation  losses  with  peak  power  (See  Section  5.7),  it  was  concluded 
that  deposited  energy  must  be  considered  as  well  as  the  peak  power. 
It  was  noted  for  the  wire  explosions  of  this  study  that  only  15-40% 
of  the  energy  initially  stored  in  the  capacitor  bank  was  deposited 
in  the  wires  before  burst.   Therefore,  when  using  the  design  method 
proposed  in  Section  3.0  one  must  insure  that  about  6  times  the  energy 
needed  to  raise  the  wire  material  to  the  desired  temperature  is  avail- 
able from  the  capacitor  bank. 

If  sufficient  energy  is  available   the  theoretical  design  method 
proposed  in  Section  3.0  is  applicable  to  fast  wire  explosions;  further 
the  experimental  results  of  this  study  indicate  that  temperatures 
within  an  order  of  magnitude  of  design  temperature  are  obtainable  from 
a  circuit  designed  by  the  theoretical  method  proposed  in  Section  3.0. 
The  experimental  results  show  that  peak  power  (theoretically  calcu- 
lated using  equation  11)   predicted  from  the  circuit  parameters  can  be 
exceeded  with  copper  wire  if  the  wire  size  is  well  matched  to  the  cir- 
cuit parameters.   (This  situation  results  because  the  wire  resistance 
can  vary  causing  an  experimental  peak  power  greater  than  that  theoreti- 
cally calculated  assuming  constant  wire  resistance.)   However,  determin- 
ing the  temperature  of  the  wire  based  on  the  peak  power  leads  to   an 
overestimation  of  the  temperature  by  about  an  order  of  magnitude.   Thus 
experimental  temperatures  attainable  for  fast  wire  explosions  may  be  an 
order  of  magnitude  lox^er  than  that  predicted  by  the  design  method  pro- 
posed in  Section  3.0. 

Although  the  wire  explosions  of  this  study  were  too  slow  to  produce 
explosive  ablative  type  wire  explosions  (See  Appendix  G) ,  more  energetic 
explosions  may  produce  this  type  of  explosion.   The  skin  effect  present 
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in  explosive  ablative  explosions  would  result  in  lower  temperatures 
than  could  be  attainable  if  the  skin  effect  were  not  present.   Since 
the  proposed  design  method  was  shown  to  underestimate  the  experimental 
temperature   for  fast  wire  explosions,  experimentally  attainable  tem- 
peratures for  explosive  ablative  type  wire  explosions  would  be  more 
than  an  order  of  magnitude  lower  than  temperatures  predicted  by  the 
theoretical  method  proposed  in  Section  3.0. 

With  respect  to  the  design  of  exploding  wire  circuits  and  experi- 
ments to  obtain  particular  desired  temperatures  it  should  be  noted 
that  several  different  combinations  of  circuit  and  wire  conditions  can 
be  used  to  obtain  the  desired  temperatures  and  energy  densities.  However, 
in  designing  the  experiments   the  circuit  parameters  determine  the  cost 
of  necessary  experimental  equipment.    Thus,  the  optimum  design  (minimum 
cost  consistent  with  experimental  objectives)  would  be  obtained  by 
selecting  the  circuit  parameters  with  theoretical  expressions  while  also 
considering  economic  factors  and  experimentally  determining  the  wire 
parameters  best  matched  to  the  selected  circuit  parameters.   Rapid  match- 
ing of  the  wire  parameters  to  the  circuit  parameters  is  possible  since 
the  qualitative  effects  of  wire  parameters  (discussed  in  Section  2.2)  are 
known . 

Experimental  results  show  that  broadening  the  peak  of  the  power 
curve  (resulting  from  longer  confinement  of  the  wire  material  after  the 
vaporization  energy  was  deposited)  would  produce  higher  energy  densities 
even  though  the  peak  may  be  slightly  lowered.   The  sample  cell  was  design- 
ed to  accomplish  this  confinement.   Although  only  limited  experimental 
results  were  obtained  in  the  sample  cell  geometry  some  useful  results  were 
noted  (See  Section  5.8).   Energy  deposited  in  a  20  mil  copper  wire  in  the 
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cell  geometry  was  greater  than  that  deposited  in  a  bare  copper  wire; 
however,  energy  deposited  in  a  confined  iron  wire  was  not  greater  than 
deposited  in  a  bare  iron  wire.   Sparkover  was  believed  to  have  limited 
both  the  peak  power  and  energy  deposited  in  the  iron  wire.   Although 
the  sample  cell  geometry  will  apparently  confine  the  wire  material  for 
longer  periods  of  time  than  wires  exploded  at  atmospheric  conditions, 
the  sparkover  problem  must  be  eliminated  for  the  sample  cell  geometry 
to  be  effective  in  obtaining  high  energy  densities. 

From  the  results  of  this  study   it  has  been  found  that  additional 
work  is  needed  before  quantitative  data  for  simulation  of  the  chemical 
environment  in  a  nuclear  explosion  can  be  obtained  by  exploding  wire 
experiments.   Considerations  in  this  additional  work  should  include: 
1)  The  conditions  which  will  produce  the  highest  energy  densities  in 
the  medium  surrounding  an  exploding  wire  should  be  determined  by  a 
parameter  study.   Specific  parameters  of  interest  would  include  wire 
size,  wire  material,  sample  cell  size,  and  material  in  the  sample  cell. 
Constraints  on  the  system  would  include:   (a)  adequate  sample  cell  de- 
sign to  prevent  sparkover,  and  (b)  the  total  material  present  in  the 
cell  would  have  to  simulate  the  environment  of  a  nuclear  explosion.   Pre- 
vention of  sparkover  was  a  major  problem  in  sample  cell  design,  particu- 
larly for  wire  explosions  with  high  peak  voltages  and  high  peak  powers . 
Therefore,  it  would  be  desirable  to  experimentally  determine  if  a  sample 
cell  which  alters  the  magnetic  pressure  on  an  exploding  wire  is  signifi- 
cantly better  than  a  sample  holder  simply  made  of  an  insulating  material. 
If  a  sample  cell  incorporating  the  magnetic  pinch  is  found  to  produce 
significantly  higher  experimental  temperatures  than  a  simple  insulating 
sample  holder,  a  sample  cell  with  improved  insulating  characteristics  will 
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be  required. 

2)  A  study  of  the  upper  limiting  energy  density  that  can  be  deposited 
in  an  exploding  wire  is  required.   Although  the  general  qualitative 
effects  (e.g. ,  An  increase  in  charging  voltage  of  the  capacitor  bank 
increases  the  energy  density  in  an  exploding  wire.)   are  known  for  fast 
explosions,  practical  limits  may  occur  in  explosive  ablative  type  wire 
explosions.   Thus,  an  analog   computer  model  to  simulate  the  mechanism 
of  an  exploding  wire  should  be  developed.   Model  design  considerations 
would  include:   (a)  an  equation  of  state  for  the  wire,  (b)  theoretical 
and  experimental  results  describing  the  energy  distribution  in  an  ex- 
ploding wire,  (c)  the  skin  effect,  (d)  the  application  of  magnetic 
fields  on  an  exploding  wire,  and  (e)  confinement  of  wire  material  by 
surrounding  material.   The  Saha  equation  should  be  an  adequate  first 
approximation  to  calculate  the  amount  of  ionization  of  the  material. 
Also,  equations  which  describe  the  application  of  magnetic  fields  to 
plasmas  for  controlled  fusion  research  should  be  applicable  to  exploding 

wires . 

41 

3)  Direct  experimental  methods  (such  as  spectrometry    or  calorimetry) 

for  measuring  the  temperature  or  energy  densities  in  the  material  of  an 
exploding  wire  experiment  should  be  used  to  check  the  accuracy  of  the 
voltage-current  method.   The  voltage-current  method  of  analysis  may  be 
limited  by  the  following  factors:   (a)  Only  deposited  energy  was 
measured  and  one  must  rely  on  incomplete  theoretical  models  to  determine 
the  energy  transferred  from  the  wire  to  the  sample  material,   (b)  Specific 
heat  and  latent  heat  data  are  incomplete  for  the  high  temperatures  in- 
volved,  and  (c)   The  energy  density  or  temperature  distribution  in  the 
sample  material  can  not  be  measured  directly  by  the  voltage-current  method, 
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Spectrometry  and  calorimetry  can  measure  more  directly  the  temperature 
or  the  energy  densities  in  the  sample  volume.   The  measurement  of 
pressure  in  the  sample  volume  would  also  be  necessary  to  completely 
specify  the  thermodynamic  conditions;  if  the  temperature  and  pressure 
are  reproducible  the  chemical  effects  should  also  be  reproducible. 
4)  Since  the  chemical  effects  to  be  simulated  in  a  nuclear  explosion 
are  far  removed  from  generally  known  chemistry  and  since  literature  on 
chemical  effects  induced  by  exploding  wires  is  very  limited,  it  would 
be  desirable  to  simulate  known  high  temperature  chemical  processes 
with  exploding  wires.   In  this  way   experimental  results  with  exploding 
wires  could  be  verified  and  the  consistency  improved.   The  high  tempera- 
ture chemistry  effects  produced  by  shock  tubes,  plasma  arcs,  and  laser 
beams  have  been  studied.   However,  results  with  these  methods  are  limit- 
ed to  low  pressures  and  are  also  limited  in  the  amount  and  type  of  chem- 
ical material  present.   Thus,  exploding  wires  with  high  pressures  and 
with  the  introduction  of  metal  in  the.  high  temperature  chemical  en- 
vironment are  not  completely  analogous  to  the  chemistry  of  shock  tubes, 
plasma  arcs,  or  laser  beams.   Therefore,  some  interpretation  of  results 
obtained  by  these  methods  in  terms  of  the  experimental  conditions  present 
in  the  exploding  wire  experiment  will  still  be  required  to  determine  the 
experimental  accuracy  obtainable  for  chemical  effects  produced  by  explod- 
ing wires. 


107 


Acknowledgments 

The  author  wishes  to  express  his  gratitude  to  Dr.  Walter  Meyer 
for  his  advice  and  constructive  criticism  during  this  study  and  to 
the  Kansas  State  University  Nuclear  Engineering  faculty  and  staff  for 
their  help.   In  particular   thanks  are  due  Mr.  R.  E.  Hightower  and  Mr. 
R.  L.  Wild  for  their  assistance  in  assembling  the  apparatus  and  their 
helpful  suggestions  concerning  the  operation  of  the  apparatus.   Thanks 
are  also  due  Dr.  Basil  Curnutte,  Jr.,  Department  of  Physics,  for  the 
loan  of  the  capacitor  bank  and  x-ray  power  supply  unit.   The  use  of  the 
Department  of  Applied  Mechanics'  traveling  microscope  is  also  appreciated , 
The  author  wishes  to  express  his  sincere  appreciation  to  both  the  Kansas 
State  University  Engineering  Experiment  Station  and  the  National  Science 
Foundation  for  their  financial  support. 


108 


LITERATURE   CITED 


W.  G.  Chace,  "A  Survey  of  Exploding  Wire  Progress,"  in  Exploding 
Wires,  Vol.  3,  edited  by  W.  G.  Chace  and  II.  K.  Moore,  pp.  1-7, 
Plenum  Press,  New  York,  (1964). 


Plowshare  Series:   Report  No.  2,  Proceedings  of  the  Second 
Plowshare  Symposium  (Held  at  San  Francisco,  May  13-15,  1959) 
Part  I:  Phenomenology  of  Underground  Nuclear  Explosions, 
UCRL-5675. 


Plowshare  Series:   Report  No.  2,  Proceedings  of  the  Second 
Plowshare  Symposium  (Held  at  San  Francisco,  May  13-15,  1959): 
Part  II:   Excavation,  UCRL-5676. 


Plowshare  Series:   Report  No.  2,  Proceedings  of  the  Second 
Plowshare  Symposium  (Held  at  San  Francisco,  May  13-15,  1959): 
Part  III:   Recovery  of  Power  and  Isotopes  from  Contained 
Underground  Nuclear  Explosions,  UCRL-5677. 


Plowshare  Series:   Report  No.  2,  Proceedings  of  the  Second 
Plowshare  Symposium  (Held  at  San  Francisco,  May  13-15,  1959): 
Part  IV:   Industrial  Uses  of  Nuclear  Explosives  in  the  Fields 
of  Water  Resources,  Mining,  Chemical  Production,  and  Petroleum 
Recovery,  UCRL-5678. 


Plowshare  Series:   Report  No.  2,  Proceedings  of  the  Second 
Plowshare  Symposium  (Held  at  San  Francisco,  May  13-15,  1959): 
Part  V:   Scientific  Applications  of  Nuclear  Explosives  in  the 
Fields  of  Nuclear  Physics,  Seismology,  Meteorology,  and  Space, 
UCRL-5679. 


Proceedings  of  the  Third  Plowshare  Symposium  (Held  at  Lawrence 
Radiation  Laboratories,  April  21-23,  1964)   Engineering  with 
Nuclear  Explosives,  TID-7695. 


Walter  Meyer  and  F.  G.  Svensson,  August  Monthly  Report,  Chemical 
Engineering  Division,  Lawrence  Radiation  Laboratory,  (1964). 


10. 


W.  G.  Chace  and  E.  M.  Watson,  "A  Bibliography  of  Electrically 
Exploded  Conductor  Phenomenon,"  Optical  Physics  Laboratory 
Project  8608,  Air  Force  Cambridge  Research  Laboratories,  L.  G. 
Hanscom  Field,  Mass.,  (1962). 

W.  H.  Richardson,  "Exploding  Wire  Phenomena,"  SCR-53,  (1958). 


109 


11.    R.  C.  Maninger,  "Preburst  Resistance  and  Temperature  of  Exploding 
Wires,"  in  Exploding  Wires,  Vol.  3,  edited  by  W.  G.  Chace  and 
H.  K.  Moore,  pp.  47-64,  Plenum  Press,  Nev;  York,  (1964). 


12.     Friingel,  Frank  B.  A.,  High  Speed  Pulse  Technology,  Vol.  1, 
pp.  379-391,  Academic  Press,  New  York,  (1965). 


13.    W.  G.  Chace  and  M.  A.  Levine,  "Classification  of  Wire  Explosions," 
J.  Appl.  Phys.  31,  1298   (1960). 


14.     F.  H.  Webb,  H.  H.  Hilton,  P.  H.  Levine,  and  A.  V.  Tollestrup, 
"The  Electrical  and  Optical  Properties  of  Rapidly  Exploded 
Wires,"  in  Exploding  Wires,  Vol.  2,  edited  by  W.  G.  Chace  and 
H.  K.  Moore,  pp.  37-75,  Plenum  Press,  New  York,  (1962). 


15.     F.  D.  Bennett  and  J.  W.  Marvin,  "Current  Measurement  and  Transient 
Skin  Effects  in  Exploding  Wire  Circuits,"  Rev.  Sci.  Inst.  33, 
1218-1226  (1962). 


16.     Jack  Katzenstein,  "The  Pinch  Effect  in  the  Exploding  Wire 

Phenomena,"  in  Exploding  Wires,  Vol.  1,  edited  by  W.  G.  Chace  and 
H.  K.  Moore,  pp.  135-144,  Plenum  Press,  New  York,  (1959). 


17.     Jan  Nasilowski,  "Unduloids  and  Striated  Disintegrations  of  Wires," 
in  Exploding  Wires,  Vol.  3,  edited  by  W.  G.  Chace  and  H.  K. 
Moore,  pp.  295-313,  Plenum  Press,  New  York,  (1964). 


18.     F.  D.  Bennett,  G.  D.  Kahl,  and  E.  H.  Wedemeyer,  "Resistance 
Changes  Caused  by  Vaporization  Waves  in  Exploding  Wires,"  in 
Exploding  Wires,  Vol.  3,  edited  by  W.  G.  Chace  and  H.  K.  Moore, 
pp.  65-88,  Plenum  Press,  New  York,  (1964). 


19.    T.  J.  Tucker  and  F.  W.  Neilson,  "The  Electrical  Behavior  of  Fine 
Wires  Exploded  by  a  Coaxial  Cable  Discharge  System,"  in  Exploding 
Wires,  Vol.  1,  edited  by  W.  G.  Chace  and  H.  K.  Moore,  pp.  73-82, 
Plenum  Press,  New  York,  (1959). 


20.    R.  C.  Maninger,  "Effects  of  Transmission  Lines  in  Applications  of 
Exploding  Wires,"  in  Exploding  Wires,  Vol.  2,  edited  by  W.  G. 
Chace  and  H.  K.  Moore,  pp.  109-126,  Plenum  Press,  New  York,  (1962) 


110 


21.     J.  K.  Trolan,  F.  M. ,  Charbonnier,  F.  II.  Collins,  and  A.  H. 
Guenther,  "A  Versatile  Ultrafast  Pulsed  Power  System,"  in 
Exploding  Wires,  Vol.  3,  edited  by  W.  G.  Chace  and  H.  K.  Moore, 
pp.  361-389,  Plenum  Press,  New  York,  (1964). 


22.     R.  Harvey  and  G.  J.  Yevick,  "Quiet,  Simple,  Inexpensive  Air-Gap 
Switch,"  Rev.  Sci.  Inst.  34,  1444-45  (1963). 


23.     H.  S.  Leopold,  "Effect  of  Bridgewire  Parameters  on  Explosive 

Initiation,"  in  Exploding  Wires,  Vol.  3,  edited  by  W.  G.  Chace 
and  H.  K.  Moore,  pp.  125-152,  Plenum  Press,  New  York,  (1964). 


24.  C.  P.  Nash  and  C.  W.  Olsen,  "Factors  Affecting  the  Time  to  Burst 
in  Exploding  Wires,"  in  Exploding  Wires,  Vol.  2,  edited  by  W.  G. 
Chace  and  H.  K.  Moore,  pp.  5-13,  Plenum  Press,  New  York,  (1962). 


25.     F.  D.  Bennett,  H.  S.  Burden,  and  D.  D.  Shear,  "Correlated 

Electrical  and  Optical  Measurements  of  Exploding  Wires,"  Phys . 
Fluids  5,  102-113   (1962). 


26.     0.  H.  Zinke,  R.  B.  Owen,  and  C.  K.  Manka,  "Trans lational 

Equilibrium  of  Wires  Exploded  in  Vacuo,"  in  Exploding  Wires, 
Vol.  3,  edited  by  W.  G.  Chace  and  H.  K.  Moore,  pp.  103-113, 
Plenum  Press,  New  York,  (1964). 


27.    M.  E.  Meeker,  "Solid  and  Solid-Liquid  Phases  in  Wires  at  High 
Current  Densities,"  SCTM  314-58(51),  (1958). 


28.     V.  V.  Bondarenko,  I.  F.  ,  Kvartskhava,  A.  A.  Pliutto,  and  A.  A. 
Chernov,  "Resistance  of  Metals  at  High  Current  Densities," 
JETP  1,  221-226  (1955). 


29.     R.  J.  Reithel  and  J.  H.  Blackburn,  "A  Hydrodynamic  Explanation 
for  the  Anomalous  Resistance  of  Exploding  Wires,"  in  Exploding 
Wires,  Vol.  2,  edited  by  W.  G.  Chace  and  H.  K.  Moore,  pp.  21-31, 
Plenum  Press,  New  York,  (1962). 


30.     C.  P.  Nash,  R.  P.  DeSieno,  and  C.  W.  Olsen,  "Electrical  and 

Emissive  Properties  of  Exploded  Confined  Conductors,"  in  Exploding 
Wires,  Vol.  3,  edited  by  W.  G.  Chace  and  H.  K.  Moore,  pp.  231-246, 
Plenum  Press,  New  York,  (1964). 


Ill 


31.    A.  Sakurai  and  T.  Takao,  "Effect  of  Applied  Axial  Magnetic 
Field  on  the  Exploding  Wire  Phenomenon,"  in  Exploding  Wires 
Vol.  3,  edited  by  W.  G.  Chace  and  H.  K.  Moore,  pp.  247-256, 
Plenum  Press,  New  York,  (1964). 


32.     F.  D.  Bennett,  "Energy  Partition  in  Exploding  Wire  Phenomena, 
Phys.  Fluids  1,  515-522  (1958). 


33.    Handbook  of  Chemistry  and  Physics,  C.  D.  Hodgman,  Editor  in 

Chief,  (Chemical  Pub.  Co.  , Cleveland ,  1959),  4Qth  ed. ,  p.  3257. 


34.    V.  E.  Scherrer,  "The  NRL-AFSWP  Exploded  Wire  Research  Program," 
in  Exploding  Wires,  Vol.  1,  edited  by  W.  G.  Chace  and  H.  K. 
Moore,  pp.  118-134,  Plenum  Press,  New  York,  (1959). 


35.     G.  J.  Dienes  and  G.  H.  Vineyard,  Radiation  Effects  in  Solids 

(Interscience  Publishers,  Inc.,  New  York,  1957),  Chap.  2,  p.  35 


36.     C.  E.  Moore,  Atomic  Energy  Levels,  U.  S.  National  Bureau  of 
Standards  Circular,  No.  467,  Vol.  Ill,  Washington,  (1958). 


37.     C.  A.  Rouse,  "Ionization  Equilibrium  Equation  of  State," 
Astrophys.  J.  134,  435-446  (1961). 


38.    K.  G.  Moses  and  T.  Korneff ,  "Voltage  Measurements  in  the  Presence 
of  Strong  Fields,"  Rev.  Sci.  Inst.  34,  849-853  (1963). 


39.    Handbook  of  Chemistry  and  Physics,  R.  C.  Weast,   Editor  in 

Chief,  (Chemical  Pub.  Co.  ,  Cleveland,  1967),  48th  ed.,  p.  D-32 


40.    K.  L.  Nielsen,  Methods  in  Numerical  Analysis  (Macmillan  Company, 
New  York,  1964),  2nd  ed.,  Chap.  4,  p.  137. 


41.    R.  Eastwood,  R.  Hysom,  and  R.  Woodruff,  "Determination  of  the 

Ion  Density  in  an  Exploding  Wire  Plasma  by  Spectroscopic  Means," 
unpublished  report  Kansas  State  University  Department  of  Physics 


42.    L.  Baker,  Jr.,   R.  L.  Warchal,  R.  C.  Vogel  ,  and  M.  Kilpatrick, 
"Studies  of  Metal-Water  Reactions  at  High  Temperatures:   I.  The 
Condenser  Discharge  Experiment:   Preliminary  Results  with 
Zirccnium,"  ANL-6257,  (May.  1961). 


112 


43.     N.  R.  Draper  and  H.  Smith,  Applied  Regression  Analysis  (John 
Wiley  &  Sons,  Inc.,  New  York,  1966),  Chap.  2,  pp.  77-81. 


44.  E.  J.  Kobetich,  "Plotting  Subprogram  for  use  on  the  IBM-1410 
Computer  and  the  IBM-1403  Printer,"  unpublished  report  Kansas 
State  University  Department  of  Physics,  (June  20,  1965). 


45.     J.  D.  Jackson,  Classical  Electrodynamics  (John  Wiley  &  Sons, 
Inc.,  New  York,  1962),  Chap.  10,  pp.  320-322. 


113 


APPENDIX  A 
Measurement  and  Calibration  of  Circuit   Parameters 

In  the  calibration  of  the  voltage  measuring  circuit   calibrated 
voltage  amplitudes  up  to  7  kilovolts  were  used;  an  analysis  of  the 
oscilloscope  photographic  records  of  voltage  measurements  showed  that 
the  accuracy  of  these  measurements  was  limited  by  the  reading  of  the 
records   with  a  Leitz   Coordinate   Comparator    with   scale 
divisions   of     10   inches.   The  precision  of  measured  distances 
was  limited  by  the  ability  to  set  the  microscope  cross  hairs  on   the 
oscillogram  traces.   By  measuring  55  distances  5  times  each,  the 
standard  deviation  of  each  reading  was  determined  to  be  .00497  inches. 
After  calibration  of  the  voltage  measuring  circuit   the  static  volt- 
meter was  calibrated  using  the  voltage  measuring  circuit;  the  frac- 
tional standard  deviation  of  the  static  voltmeter  was  found  to  be 
.011  volts  per  volt. 

Experimental  data  to  determine  the  capacitance  of  the  discharge 
circuit  were  obtained  by  discharging  the  capacitor  bank  through  cali- 
brated resistors  (Tru-ohm) .   Voltage-time  oscillogram  records  were  ob- 
tained  by  measuring  the  voltage  across  5  x  10   ohms  resistance.   The 
vertical  position  of  each  voltage  trace  was  measured  by  the  coordinate 

comparator  at  horizontal  intervals  of  At/R  ,  =  10   sec/ohm.   The  ver- 

cal 

tical  position  of  the  zeroline  reference  was  also  measured  for  each 
oscilloscope  record. 

A  weighted  least  squares  fit  in  the  logarithmic  domain  was  used 
to  calculate  the  capacitance  (See  Appendix  B) .   Since  the  positions 
measured  at  a  large  number  of  time  intervals  after  the  initiation  of 


114 


discharge  through  a  calibrated  resistor  approached  the  zeroline 
reference  position,  a  least-squares  fit  using  too  many  points  did  not 
give  accurate  results.   Switch  effects  were  also  significant  when  the 
voltage  remaining  on  the  capacitor  bank  was  small.   Therefore,  least- 
square  fitting  of  the  data  was   accomplished  using  3  to  11  data  points 
and  the  F  test  number  was  used  to  determine  the  number  of  points  which 
gave  the  best  fit  for  each  oscilloscope  record. 

Five  oscillograms  were  photographed  on  each  of  five  separate 
days  for  a  capacitor  bank  discharge  with  8  capacitors  (Nominally  1  yF 
each  and  initially  charged  to  8  kV)  in  parallel.   The  values  of  capac- 
itance which  gave  the  best  fit  (F  test  number)  for  each  of  the  25 
oscillograms  were  averaged  to  obtain  the  capacitance  of  the  capacitor 
bank.   With  4  capacitors  in  parallel,  7  oscillogram  records  were  ob- 
tained.  The  7  records  were  obtained  during  two  days.   The  same  pro- 
cedure was  used  to  determine  the  capacitance  which  gave  the  best  least- 
squares  fit  for  each  trace  and  the  results  of  7  traces  were  averaged. 
The  variance  of  the  capacitance  calculated  using  equation  (B-6)  was 
the  result  of  the  variance  in  the  fitting  of  one  curve,  but  significant 
variation  was  expected  to  possibly  occur  between  different  discharges. 
Thus,  the  variance  of  the  capacitance  was  estimated  using 

N 
a2(C)  =  (JL)   (|)     Z   (C  -C.)2   ,  (A-l) 

i=l 

and  also  using  equation  (B-6)  together  with  the  propagation  of  variance 

formula.   The  hypothesis  of  no  variation  of  the  capacitance  from  day 

to  day  was  accepted  since  the  Student's  t-test  did  not  show  differences 

greater  than  the  95%  confidence  range.   Also,  there  was  insufficient 

evidence  to  show  that  the  variance  calculated  by  equation  (A-l)  was  sig- 
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nificantly  (F  test  -  95%  confidence  level)  larger  than  the  variance 
calculated  by  equation  (B-6)  and  the  propagation  of  variance  formula. 
However,  the  estimated  variance  calculated  by  equation  (A-l)  was  the 
value  used  to  calculate  the  reported  standard  deviation  for  capaci- 
tance.  The  following  values  were  calculated  for  the  capacitance  at 
an  initial  charging  voltage  of  8  kV: 

Number  of  Total  Capacitance 

Capacitors  (microfarads) 

4  3.587  +  .042 

8  7.712  +  .051 

The  capacitance  was  also  determined  for  eight  capacitors  with  an 
initial  charging  voltage  of  11  kV;  the  calculated  value  of  capaci- 
tance for  each  of  these  discharges  was  within  the  standard  deviation 
of  a  single  observation  of  the  values  reported  above. 

The  current  waveform  was  used  to  determine  the  extraneous  cir- 
cuit resistance  and  the  inductance  of  the  exploding  wire  circuit.   The 

current  in  a  circuit  containing  only  resistance  and  inductance  when  a 

42 
charged  capacitor  is  connected  to  the  circuit  has  the  following  form: 

V 
i         °   [exp  (-Rt/2L)]  [sin  /(1/LC)  -  R2/4l2   t]   (A-2) 


(L/C)-(R^/4) 


2 

for  the  underdamped  case  of  R  <  4L/C 


where: 

i  is  the  current  in  kA, 


V  is  the  charging  voltage  in  kV, 

R  is  the  total  circuit  resistance  in  ohms, 
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L  is  the  total  circuit  inductance  in  henries, 
C  is  the  capacitance  in  farads,  and 
t  is  the  time  in  seconds. 

2 
For  the  case  of  R  <<4L/C   equation  (A-2)  can  be  expressed  approxi- 
mately as  follows: 

i  S  K[exp  (-Bt)][sin  (wt)]  (A-3) 


where: 

1/2 
K  =  V  (C/L)  '         , 

B  =  R/2L   ,  and 

a>  =  (1/LC)1/2   . 

For  the  exploding  wire  circuit  used  in  this  study  the  total 
circuit  resistance  was  the  sum  of  the  wire  resistance  and  the  ex- 
traneous circuit  resistance.   With  a  shorting  bar  substituted  for 

the  wire,  the  total  circuit  resistance  was  the  extraneous  circuit 

2 

resistance  and  R  was  much  less  than  4L/C.   Thus,  the  period  of  the 

oscillation  was 

t  =  2tt(LC)1/2   .  (A-4) 

Fifty-seven  oscillogram  records  were  used  to  measure  the  period  with 
the  capacitance  of  7.712  +  .051  microfarads.   The  arithmetic  average 
of  the  measured  periods  was  20.11  +  .03  microseconds.   (The  error 
limit  reported  is  the  estimated  standard  deviation.)   The  variance  of 
the  period  of  the  circuit  with  the  shorting  bar  in  place  was  estimated 
using 
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a2(T)  =  (^)   (|)     I  (t.  -a)2   .  (A-5) 

i=l 

where: 

N  -  57 

x  =  arithmetic  average  of  the  measured  periods  (20.11  psec) . 

The  circuit  inductance  was  calculated  from  the  measured  values  of 
the  period  and  the  previously  determined  value  of  capacitance  using 

L  =   (t/2t02   (1/C)    .  (A-6) 

The  variance  of  the  circuit  inductance  was  determined  using  the 
propagation  of  variance  formula: 


2/T,N         ,3F   .2        2.      v.    ,3F  N2      2.      N    ,  ,     ,  3F   .2      2.      N  ,.    _. 

o   (F)    =    (-3-)        a  (X]L)   +   (-£-)      a   (x2)  +   ...   +   (-£-)      a   (xm)      .      (A-7) 

12  m 


where: 


F  =   F(x    ,   x    , . . . ,x   ) ,    and 

LA  m 

x.;  i  =  l,2,...,m  are  independent  variables. 


Applying  equation  (A-7)  to  (A-6)  and  carrying  out  the  partial  differ- 
entiation yields 

a2(L)  =  (t/2tt2C)2  o2(t)  +  (t/2ttC)4  a2(C)   .  (A-8) 

Thus,  the  calculated  value  of  inductance  is  1.33  +  .01  uH. 

The  peak  values  of  current  are  when  t  =  i(2n  -  l)/4  where  n  is 
a  positive  integer.    The  value  of  the  current  at  the  n    peak  is   given 
by 
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t      .n+1    ,     .      . ,  -(2n-l)/4  /A  n. 

1   =  (-1)     K[exp  (Bt)J  .  (A-9) 

Therefore,  if  A  is  the  absolute  value  of  the  current  at  the  n 
n 

peak  and  is  measured  on  the  oscilloscope   photographic  record  using 
distance  units,  then 

A       T    A       -  ry 

X  =  a"     a   —  =  exP  (fiT)  =  exP  (2tt  RC/t)  .    (A-10) 
n    n+2  +  n+3 

The  effective  circuit  resistance  from  peak  n  to  peak  n+3  was  calcu- 
lated from  the  following  equation: 

R  =  (t/2tt2C)   In  X  (A-ll) 

n  n 

for  values  of  n  =  1,2,  ...,  7  for  each  recorded  oscillogram.  The 
variance  of  the  circuit  resistance  was  calculated  using  the  propa 
gation  of  variance  formula: 

a2(R  )  =  (R  /C)2  a2(C)  +  (R  /t)2  a2(x)  +  (R  /X  In  X  )2  a2(X  )   (A-12) 
n      n  n  n  n    n       n 


where: 


a2(Xn)  =  (  1  +  xn2)/(An+2+An+3)2  ^ {k  )    '  and 


2  2 

a  (A)  =  (.00497  in)   is  the  reading  variation. 


The  computer  program  described  in  Appendix  C  was  used  to  calcu- 
late values  from  equations  (A-ll)  and  (A-12)  and  the  output  of  the 
computer  program  was  used  to  calculate  a  weighted  average  resistance 
for  each  discharge  using 

7 

Z   w  R 

R.  =        n=l   n  n  (A-13) 
_ 

Z   w 
n=l 
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where : 

w  =  (1/  a2(R  )  )  . 
n  n 

Fifteen  oscillograms,  obtained  from  discharges  of  the  capacitor  bank 
through  a  shorting  bar,  were  measured  for  current  amplitudes  at  the 
current  peaks.   Results  for  each  oscillogram  were  calculated  as  de- 
scribed above  and  the  results  of  the  fifteen  oscillograms  were  aver- 
aged to  determine  the  extraneous  circuit  resistance; 

R  =  .0351  +  .0006  ohms.   The  standard  deviation  reported  is  the 
e         — 

square  root  of  the  variance  which  was  calculated  using 

1=1 

where: 

N  =  15,  and 

R  is  the  arithmetic  average  of  the  extraneous  circuit  resistance 

e  ° 

2  — 
of  fifteen  oscillograms.   The  variance,  o    (R  ) ,  was  also  calculated  by 

applying  the  propagation  of  variance  formula  to  equation  (A-13)  and  to 
the  averaging  of  the  fifteen  oscillograms.   The  F  test  showed  that  the 
variance  calculated  using  equation  (A-14)  was  significantly  larger 
than  the  variance  calculated  by  the  propagation  of  variance  formula. 
This  indicates  that  significant  variation  of  the  measured  current  occur- 
red between  different  discharges. 

The  output  of  the  current  measuring  circuit  was  calibrated  using 
the  same  oscillogram  measurements  which  were  used  to  calculate  the  cir- 
cuit inductance  and  the  extraneous  circuit  resistance.   Equation  (A-9) 
yields 
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v/   ,vn+l   ,,    ,_.  .,-(2n-l)/4    r    ,D  N1-(2n+l)/4  , 
i   -  i  +1  =  K(-l)      {[gxp(Bt)]  +  [exp(Bi)]  }  . 

(A-15) 


Substituting  from  (A-10)  yields 

i   -  i  x,  =  K  (-l)n+1   [  X-(2n-1)/4  +  X-(2n+1)M  ]  .  (A-16) 

n    n+1  ' 

The  measured  value  on  the  photographs  which  corresponded  to  the 

current,  (i  -  i  ,  ,),  was  the  distance  (A  +  A  . . )  which  was  the 
n    n+1  n    n+1 

vertical  distance  between  peaks  n  and  n  +  1  on  the  photographic 
records.   Thus,  the  vertical  scale  of  the  current  waveform  was 
calculated  using 

t  n  ma-     •  •   ^    r        K   %    r  v-(2n-l)/4,   -(2n+l)/4  , 
I   (kA/dxvision)  =  (.  ■  ■. )    [X         +  X  j  . 

n  n+1 

(A-17) 

for  n  =  1,2,.. .,7.   The  computer  program  described  in  Appendix  C 

used  equation  (A-17)  and  the  propagation  of  variance  formula  applied 

2 
to  equation  (A-17)  to  calculate  values  of  I   and  c  (I  )   for 

n         n 

n  =  1,2,.. .,7. 

The  output  of  the  computer  program  described  in  Appendix  C  was 

used  to  calculate  a  weighted  average  for  each  oscilloscope  trace  using 

7 

2  w  I 
I.(kA/division)  =     n=l  n  n  (A-18) 


i 


7 

I  w 

i   n 
n=l 


where 


n   a  (In) 
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Results  for  each  oscillogram  were  calculated  as  described  above  and 
the  results  of  the  fifteen  oscillograms  were  averaged.   From  these 
calculations  the  vertical  scale  of  the  current  waveforms  was  deter- 
mined to  be  I  =  19.44  +  .65  kA/photo  division  for  the  following 

circuit  parameters;  R.    =5  kohm,  C.    =  .001  uF,  and  switch  setting 
r  mt  int 

of  amplifier  L  =  20  volts /cm. 
r        s 

The  standard  deviation  of  a  single  observation  for  the  current 
measuring  circuit  was  calculated  by  the  propagation  of  variance  form- 
ula and  by  the  following  equation: 

a2 (I)  =  C—j)       I   (I  -  I.)2   ,  (A-19) 

i=l 

where: 

N  =  15,  and 

I  is  the  arithmetic  average  of  the  current  scale  which  was  cal- 
culated from  fifteen  oscillograms.   Since  the  standard  deviation  of 

each  current  measurement  was  of  interest,  the  standard  deviation  of  a 

2 

single  observation,  a    (I),  was  determined  instead  of  the  standard 

2  — 
deviation  of  the  mean,  a   (I).   Similar  to  the  results  of  the  extran- 
eous circuit  resistance  the  F  test  showed  that  the  variances  calculated 
by  the  two  methods  differed  significantly.   The  variance  calculated 
using  equation  (A-19)  gave  larger  variance  and  was  the  value  used  to 
calculate  the  reported  estimate  of  the  standard  deviation  of  a  single 
observation. 

The  vertical  scale  for  other  integrating  circuit  parameters  was 
given  by 

I'  (kA/photo  division)  =  (10  R.    C.    V  )  I     (A-20) 

s 
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where: 

R.    is  the  resistance  of  the  integrating  circuit  (kohrn)  , 

C.    is  the  capacitance  of  the  integrating  circuit  (jjF)  >  and 

V    is  the  switch  setting  of  amplifier  L  on  the  current 
s 

measuring  circuit. 

The  time  derivative  of  the  current  in  an  RLC  circuit  for  the 

underdamped  case  has  the  following  form: 

(di/dt)        *   -  (V  /L)  [ex?  -(Rt/2L)  ]  [cos  /T^"r2   t  ].  (A-21) 
/   R2C  °  LC   4l7 

4L 

where  the  notation  is  the  same  as  in  equation  (A-2) .    For  the  case 

2 
of  R   <<4L/C  equation  (A-21)  can  be  expressed  approximately  as  follows: 


(di/dt)  -  (V  /L)  [exp   (-Rt/2L)]  [cos  (at)]  .        (A-22 
o 

An  analogous  procedure  was  used  to  calibrate  the  (di/dt)  wave- 
form as  was  used  to  calibrate  the  current.   The  time  derivative  of  the 
current  at  the  m   peak  is  given  by 

(di/dt)   =  (-l)m_1   (V  /L)   [exp  -(rn-l)Bx/2]   .       (A-23) 
m  o 

The  quantity  Y  was  defined  by 

Y  =   Dm  +  Dm+1 

m+2    m+3 

where  D   is  the  absolute  value  of  the  current  derivative  at  the  m 
m 

peak  and  is  measured  on  the  oscilloscope  record  using  divisions  on  the 
photograph  as  units.   The  vertical  scale  of  the  (di/dt)  waveform  was 
determined  using  the  following  equation: 
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/  ,-*-  /  ,  \   /i  «     -1  ,  , .  .  .   N        \r  r„-  (m-1)  /2  .   -m/2n 

(dl/dt)   (kA. sec   /division)  =      VQ      [y         +  Y    ] 

m 


L(D  +  D    ) 
m    m+l 


(A-25) 


for  m  =  1,2,...,7.   A  weighted  average  of  the  vertical  scale  of  the 

current  derivative  on  each  oscilloscope  record  was  calculated  using 

7 

E  w   (dl/dt) 
(dl/dt).  (kA-sec   /division)  =   m-1  ™ m       (A-26) 


l 


7 

E  w 

i   m 
m=l 


where: 


w  = 
m 


a2 (dl/dt) 

m 


The  computer  program  described  in  Appendix  D  was  used  for 
calculations  in  calibrating  the  circuit  for  measuring  the  time  deriv- 
ative of  current.   Equations  (A-24) ,  (A-25),  and  (A-26)  were  used  in 
the  computer  program  described  in  Appendix  D.   The  propagation  of 
variance  formula  was  also  applied  to  the  above  equations  and  programmed 
in  the  computer  code.   Results  for  each  oscillogram  were  calculated 
as  described  above  and  the  results  of  ten  oscillograms  were  averaged. 
The  estimated  variance  of  a  single  observation  was  calculated  by  the 
propagation  of  variance  and  also  by  using 


2  1      10   2 

a  (dl/dt)  =  fc-)     I    [  (dl/dt)  -  (dl/dt).  T   .     (A-27) 

i=l  x 


-N-l' 


where : 

N  =  10,  and 
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(dl/dt)  is  the  arithmetic  average  of  the  current  derivative 
scale  of  ten  oscillograms. 
The  F  test  showed  that  the  variance  calculated  using  equation  (A-27) 
was  significantly  larger  than  the  variance  calculated  by  the  propa- 
gation of  variance  formula.   The  results  of  the  calibration  of  the 

current  derivative  measuring  circuit  were: 

9      -1 
dl/dt  =  (5.42  +  .12)  X  10  A«sec   /photo  division,  which  were  valid 

for  the  following  amplifier  settings;  V   =  20  V/cm,  and  VT   =1  V/cm. 

P  s 

The  vertical  scale  for  other  time  derivative  circuit  parameters 
for  the  measuring  circuit  was  given  by 
(dl'/dt)  =  (V  V   /20)   (dl/dt)  A. sec~1/photo  division.     (A-28) 

Ij    Ij 

p  s 
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APPENDIX  B 

Description  and  Explanation  of  IBM-1410  Computer  Program  for 
Calculating  Capacitance  of  Discharge  Circuit 

The  free  response  of  an  overdamped  resistance-capacitance  network 
has  the  form 


VR(t)   =  Vq  exp(-t/RC) 


(B-l) 


Transforming  to  the  natural  logarithmic  domain  gives 

z  =  z  +  (t/R)  a  (B-2) 

where  the  quantities  are  defined  as; 
z  =  In  V(t)  ,  and 
a  =  -(1/C)  . 

Written  in  matrix  form,  the  weighted  least  squares  fit  of  data  to 

43 
equation  (B-2)  is 


T   -1   T 
b  =  (T  WT)    II  z, 


(B-3) 


where 


z 

o 


z  = 


\  : 

z 


v.    j 


T  = 


1  t1/R 
1  t2/R 


lt/R 
n 


,  and 


W  = 


w   0   0 
0   w„  0 


0   .   . 


.  0 
.  0 

w 


n 


The  variance  of  the  estimate  is 
2 
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T 

z  Wz 


T   T 
b  (T  W)z 


est    //  of  degrees  of  freedom 
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(B-4) 


The  weights  of  z  used  in  the  calculations  were 
w.   =   (1/a  (z . ))   =     i 


(B-5) 


a2(V.) 


The  variance  of  the  voltage,  a  (V.),  was  determined  by  the  variance 
of  measuring  distances  on  the  photographs  with  the  traveling  micro- 
scope (See  Appendix  A) . 

The  variance  of  the  capacitance  for  the  least -squares  fit  was 

43 
calculated  from 


o2(C) 


°4  A22  aest    ' 


(B-6) 


where 


All 

A12 

A21 

A22 

J 

T    -1 
(T  WT) 


The  above  equations  were  programmed  in  FORTRAN  IV  language 
for  the  IBM  1410  computer.   The  source  program  is  listed  and  the  logic 
diagram  is  shown  in  this  appendix.   Table  B-l  defines  important  symbols 
used  in  the  program. 

The  input  to  the  program  was  arranged  as  follows: 
TIME  DATA:   FORMAT (8F 10 .0) 

T(J):   Time  parameter  data  (actual  time  from  the  start  of  the 
oscilloscope  trace  divided  by  the  calibrated  resistance 
in  the  circuit)  at  which  voltage  measurements  were  ob- 
tained. Data  were  read  in  with  eight  values  per  card. 
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Table  B-l:   Symbols  used  in  the  IBM-1410  computer  program  for  calcu- 
lating capacitance  of  discharge  circuit 


Symbol  Explanation 

CAP  Least  squares  fit  value  of  capacitance 

F  F  test  number 

I  I   time  the  photograph  was  measured 

J  J   position  on  the  photograph 

K  Number  to  identify  photograph  (data  set  for  sub-run) 

SIGC  Estimated  standard  deviation  of  capacitance  from  least 
squares  fit 

SIGV  Estimated  standard  deviation  of  initial  voltage  on  capaci- 
tor bank  from  least  -squares  fit 

SUMPC  Estimated  variance  removed  by  the  correlation 

T(J)  Time/calibrated  resistance  in  circuit 

V(I,J)  Voltage  on  photograph  at  time  T(J)  for  measurement  //  I 

VARES  Variance  of  the  estimate 

VO(I)  Zero  voltage  reference  for  the  I   measurement 

VZERO  Least  squares  value  of  initial  voltage 

WW  Weight  of  voltage  at  time  T(J) 

X  Reading  variance  of  voltages 

Z  Natural  logarithm  of  voltage  V(J) 
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CARDS  USED  FOR  SUB-RUNS 
FIRST  CARD:   FORMAT (212) 

IX:   Number  of  times  photograph  was  measured. 
JX:   Number  of  data  points  taken  on  each  photograph. 
SECOND  CARD:   FORMAT (5F10.0) 

VO(I):   Position  of  zero  line  for  I   measurement  of  the 
photograph. 
REMAINING  CARDS:   FORMAT (5F10 . 0) 

V(I,J):   Measured  position  of  voltage  at  time  corresponding 

to  T(J)  for  I   reading.   The  position  of  the  voltage 
at  T(l)  for  five  measurements  was  punched  on  the  first 
card.   The  following  cards  had  positions  of  voltage  at 
T(J)  consecutively  for  J  =  2, 3,..., 11.   Each  card 
had  five  punched  values  for  five  measurements  of  photo- 
graph. 
(Data  for  multiple  sub-runs  were  fed  in  order  with  one  set  of 
cards  for  each  of  the  sub-runs.   JX  for  each  sub-run  must  be 
less  than  or  equal  to  the  number  of  points  read  in  for  the  time 
data.) 
The  output  of  the  program  was  in  printed  form.   Values  printed 
were:   K,  J(number  of  points  used  in  least  squares  fit),  CAP,  SIGC, 
VZERO,  SIGV,  and  F. 
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Figure  37.   Logic  diagram  for  capacitance  of  discharge 
circuit  computer  program. 
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Invert 

Matrix: 

T    -1 
(TWT)  X 


Calculate: 
VZERO 
CAP 


Calculate: 
SIGV 
SIGC 


Output : 
K,J,CAP 
SIGC, VZERO 
SIGV,F 
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APPENDIX  C 

Description  and  Explanation  of  IBM-1410  Computer  Program  for 

Calculating  Extraneous  Circuit  Resistance  and  for 
Calculating  the  Calibration  of  the  Current  Measuring  Circuit 

This  computer  program  was  written  to  calculate  the  extraneous  cir- 
cuit resistance  of  the  discharge  circuit  and  to  calculate  the  calibration 
of  the  current  measuring  circuit.   The  experimental  procedure  used  to  ob- 
tain the  data  is  discussed  in  Section  4.2,  and  the  equations  used  in  the 
calculations  are  discussed  in  Appendix  A.   Equations  used  in  the  calcu- 
lations include  (A-17) , (A- 10) ,  (A-ll) ,  and  (A-12) .   These  equations  were 
programmed  in  FORTRAN  IV  language  for  the  IBM-1410  computer.   The  source 
program  is  listed  and  the  logic  diagram  is  shown  in  this  appendix. 
Table  C-l  defines  important  symbols  used  in  the  program. 

The  input  to  the  program  was  arranged  as  follows : 
FIRST  CARD:   FORMAT (8F10.0) 

1)  T:   Period  of  oscillation,  2)  C:   Capacitance  of  capacitor  bank, 
3)   VO:   Initial  charging  voltage  on  capacitor  bank,   4)  VART: 
Estimated  variance  of  measured  capacitance,   6)  VARV:   Estimated 
variance  of  static  voltmeter  measurements,   &)   VARA:   Estimated 
reading  variance  of  photograph  measurements,  and   8)  F:   Conversion 
factor  from  measuring  units  to  photograph  units. 
CARDS  USED  FOR  SUB-RUNS 
INITIAL  CARD  FOR  EACH  SUB-RUN:   FORMAT (212) 

1)  L:   Number  of  first  current  peak  entered  as  data,  and 

2)  N:   Number  of  last  current  peak  entered  as  data. 
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Table  C-l:   Symbols  used  in  the  IBM-1410  computer  program  for  calcu- 
lating the  extraneous  circuit  resistance  and  for  calculat- 
ing the  calibration  of  the  current  measuring  circuit. 


Symbol  Explanation 

A(I)  Vertical  position  of  1   current  peak  as  measured  by  the 
traveling  microscope 

AMIC  Current  per  microscope  division 

APH  Current  per  photo  division 

C  Capacitance  of  capacitor  bank 

DEVD  Estimated  standard  deviation  of  DN 

DEVM  Estimated  standard  deviation  of  AMIC 

DEVP  Estimated  standard  deviation  of  APH 

DEVR  Estimated  standard  deviation  of  RE 

DN  Difference  in  value  of  current  between  peaks  I  and  1+1 

F  Conversion  factor  (.363  microscope  divisions  per  photo 
division) 

I  Number  to  identify  current  peaks  of  current  waveform 

J  Number  to  identify  photograph  (data  set  for  sub-run) 

L  Number  of  first  current  peak  entered  as  data 

N  Number  of  last  current  peak  entered  as  data 

RE  Effective  extraneous  circuit  resistance  from  peak  I  to 
peak  1+3  for  SB  //  J 

T  Period  of  discharge  circuit  with  shorting  bar 

VARA  Estimated  reading  variance  of  photograph  measurements 

VARC  Estimated  variance  of  capacitance 

VARD  Estimated  variance  of  DN 

VARR  Estimated  variance  of  RE 

VART  Estimated  variance  of  period  of  discharge  circuit 

VARV  Estimated  variance  of  static  voltmeter  measurements 

VO  Initial  charging  voltage  on  capacitor  bank 
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REMAINING  CARDS  FOR  EACH  SUB-RUN:   FORMAT (8F10.0) 

A(I):   Measured  vertical  position  of  I   current  peak  as  measured 
by  the  traveling  microscope.   (Data  for  multiple  sub-runs  were 
fed  in  order  with  one  set  of  cards  for  each  of  the  sub-runs.) 

The  output  of  the  program  is  set  up  for  printed  form.   Values 
printed  were:   RE,  DEVR,  VARR,  DN,  DEVD,  VARD,  AMIC,  DEVM,  APH,  and 
DEVP.   Some  of  these  values  were  used  for  further  calculations  as  de- 
scribed in  Appendix  A. 


135 


THIS    Po       RAv    COMPUTES    THF    C':      EMI     CALIBRATION    Of      1  ?ENT 

■  •  -      t    ir  T  N G    ' ' r  "1  '    10     "I  • ' r     R  :  '   I  c.  1  •'  :  C  F    A  M    U N D E C        '     ' 

PLC     N'  E  T  'a'  0  P  K 


L-rkr  ,nm'C.ENGG.  ]  .*'■'!*'•  » 1  op  ACE  5 


>"-,v   | 

|-t 

r  u  r  r  A  L 

. .,  v  -  a 

r  z  ■■  j 

ol  I vf  ' 

•  "">  r  ■  i  " 

A  C  '-,  •  l 

WJR,1 ? 

..  -  .  ^g 

'■  s  G  N 

MGC» 16 

WON ':  f 

'••'ODE 

G  0  ♦  T  F  S  T 

Vp:V 

FXFQ 

FCRTR/ ' 

OIVf;^clC^ 

1  M  ? . 

) 

,  ,  .  ,  ,  ,CURCAL 


or-    FORMAT  (  8F1C.0  ) 
^i     FORMAT (2  12) 

oq    FORMAT ( 1HL»10X »32HSHGRTIN<  '° 

219HPH~TO.     NUMBER  S8»I3) 

1  r  n    F  OR  M  A  t  (  1HK»I2»5X»22HPEC  I S  T  A  f*  C  F 

?r  1  2.8  . -'  X  »  4HV/  R=»F1  2.8  ) 
j  P  i     p  or  v  At(]X»I2»6X»2-,HPFAK     TO    P  r  A  y    CURRF\iT  =  »F12.7*5X»4HDfr\/: 
2^X»  4H\/AR  =  ,F12.8  ) 
?    FOP'"AT  (  i  X  ,  I  2  *  ?7uri  irrfivt    ppR    MICRO    D  I V I  S  I  OM=  >  F  1 2  •  7  »  5  X  »^HDEV=  » F 12  •  «  ) 


CURRENT    CALI  bRAT  I  0N»  10X» 

JIT=»F12.7,5X,4H  FV=, 

»  c  1  ?  .  8  , 


1 


F12.7»5X »4HDEV= »F12.8 ) 


FORMAt( IX » I2»27HCMRRFNT  PER  PHPTO  DIVISION: 

READ( i ? 90 )  T»C»V  »VART»VARC »VARV»VARA ,F 

D  0 1  0  J  =  ]  ,  1 0  C 

R  E A  D (1 , 9 1  )  L  ,  N 

READ{ i ,90 ) ( A( r ) , I=L»NJ 

•/■.  -.   "_o 

PI =3. 1415926536 

p ;-D!¥DJ 

D010I=L»MN 

I] = I + 1 

I  2= I  +  i 

13=1+1 

XMUM=*  ( I ) -A(  I]  ) 

XDEN=a ( I  2  )-A(  I  3  ) 

X=XNUm/XDEN 

XX  =  AL^G( X  ) 

PF  = ( T-XX) / ( 2.*DS*C ) 

VARX=( VARA+ ( X*X*VARA ) ) / ( XDEN*XDEN) 

V  a  RR=  t  {  VARC/  (  C*C  )  )  +  (  VART/  (  T*T  )  )  +  (  VARX/  (  XX*XX*X*X  )  )  )*  (RE 

D  F  V  R  =  c  Q  R  T  (  V  A  R  R  ) 

X<=2.^*PI*VC*c/T 
XT  =  I 

D  « ■  =  -  (  (  X  I  /  ? .  ~ )  -  .  ?  c  ) 

P  P  =  X  #  *  r    -' 

dC=X**(PA-0.5  ) 
p  0  =  X  *  *  (  P  A  - 1 .  C  ) 
PF=X**( PA-l.5 ) 
D|g  =  XK*  (P3  +  PC  ) 

VARK=(  (4.*PS*C*C*V  \RV) /(T*T ) )+(  (4.*PS*V<  *V<  *VARC)  / ( T*T  ) 
2+  (  (  4.*PS*VC  *  V  *C*C*VART  )  /  (  T**4.  )  ) 


R  F 


136 


V 

2  ■•■' 


i  o 


r 


•  ■  ; 


■ 

(  ( 
M  T 

DH 

FV 

FV 
r<  I 
PI 
RI 
PI 
PI 

<~  M 

TO 
N1 


D  a  ..  i    i 
'  '  /X 

D  =  c  0  R  T 
f/=M  v 
P  =  ^f  V-' 
T  E  (  3  >  9 
T  c:  f  3  »  1 
TE ( 3  » 1 
T  E  f  3  » 1 
T  F  f  ^  » 1 

T  I  mm: 
P 


+  or  )  -*  (  P  p  +  PC)  )  +  (  V  A  P  >   x  K  *  X  K  *  {  (  ?  '  '   '  )  +  {(PA-^.5)#P'  ' 

-  r  .  «=  )  ■  ;  '  )  )  ; 


.'  <•  (  1 1 
+  ( ( p 

M  ' 

F 

(  V  A  P 

/XNU 

*F 

9)  J 

)  I 
LI)  I 
•'  ?  )  I 
3)1,/  PH»r"  VF 


0  ) 


,  r;r  ,DEVR>VARR 
,ON»DEV[  »\"  '■  i 
9  rs v  j  r  ,  o  r  V ' ' 


FXEQ  L I v"l OAD 

CALL  CURCAL 

EXEC  CURCAL»MJB 

JOB   ACTSSCLIVER  BLOCK 


MUC  ENGG   C2  37C4CM2 


137 


Figure  38.   Logic  diagram  for  extraneous  circuit  resistance 
and  current  calibration  computer  program. 


C       Start:     j 


Calculate: 
RE,DEVR,VARR 
DN,DEVD,VARD 
AMIC,DEVM 
APH,DEVP 


Output : 
R£,DEVR,VARR 
DN,DEVD,VARD 
AMIC.DEVM, 
APH,DEVP 


Read: 
T , C , VO , VART 
VARC,VARV 
VARA,F 
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APPENDIX  D 

Description  and  Explanation  of  IBM-1410  Computer  Program  for 
Calibration  of  the  Current  Derivative  Measuring  Circuit 

This  computer  program  was  written  to  calibrate  the  current 
derivative  measuring  circuit.   The  experimental  procedure  used  to  ob- 
tain the  data  is  discussed  in  Section  4.2,  and  the  equations  used  in 
the  calculations  are  discussed  in  Appendix  A.   Equations  used  in  the 
calculations  include;  (A-24) ,  (A-25) ,  and  the  propagation  of  variance 
formula.   These  equations  were  programmed  in  FOR.TRAN  IV  language  for 
the  IBM-1410  computer.   The  source  program  is  listed  and  the  logic  dia- 
gram is  shown  in  this  appendix.   Table  D-l  defines  important  symbols 
used  in  the  program. 

The  program  described  in  this  appendix  and  the  program 
described  in  Appendix  C  have  nearly  the  same  form;  the  main  difference 
is  the  interchange  of  statements  so  that  equations  for  the  current 
derivative  calibration  are  used  in  this  program  and  equations  for  the 
current  calibration  and  extraneous  circuit  resistance  are  used  in  the 
program  described  in  Appendix  C.   (Note  also  that  the  program  described 
in  this  appendix  included  the  calculation  of  the  weighted  average, 
equation  (A-26) ,  for  each  oscilloscope  trace.   The  weighted  average  was 
calculated  using  a  hand  calculator  for  the  current  calibration.) 

The  input  to  the  program  was  arranged  as  follows: 
FIRST  CARD:   FORMAT (8F10.0) 

1)  T:   Period  of  oscillation,  2)  C:    Capacitance  of 
capacitor  bank,  3)   VO:   Initial  charging  voltage  on  capac- 
itor bank,  4)  VART :   Estimated  variance  of  measured  period, 
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Table  D-l:   Symbols  used  in  the  IBM-1410  computer  program  for  calculat- 
ing the  calibration  of  the  current  derivative  measuring 
circuit . 

Symbol  Explanation 

A(I)      Vertical  position  of  I    current  derivative  peak  as  measured 
by  the  traveling  microscope 

AMIC      Current  derivative  per  microscope  division 

APH       Current  derivative  per  photo  division 

B         Weighted  average  of  current  derivative  per  photo  division  as 
calculated  from  Equation  (A-26)  for  each  trace 

BB        Estimated  standard  deviation  of  B  which  was  calculated  using 
the  propagation  of  variance  formula 

C  Capacitance  of  capacitor  bank 

DEVD  Estimated  standard  deviation  of  DN 

DEVM  Estimated  standard  deviation  of  AMIC 

DEVP  Estimated  standard  deviation  of  APH 

DN        Difference  in  value  of  current  derivative  between  peaks  I 
and  1+1 

F  Conversion  factor  (.363  microscope  divisions  per  photo  division) 

I  Number  to  identify  peaks  of  current  derivative  waveform 

J  Number  to  identify  photograph  (data  set  for  sub-run) 

L  Number  of  first  current  derivative  peak  entered  as  data 

N  Number  of  last  current  derivative  peak  entered  as  data 

T  Period  of  discharge  circuit  with  shorting  bar 

VARA  Reading  variance  of  photograph  measurements 

VARC  Estimated  variance  of  capacitance 

VARD  Estimated  variance  of  DN 

VART  Estimated  variance  of  period  of  discharge  circuit 

VARV  Estimated  variance  of  static  voltmeter  measurements 

VO  Initial  charging  voltage  on  capacitor  bank 
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5)  VARC:   Estimated  variance  of  measured  capacitance,  6) 
VARV:   Estimated  variance  of  static  voltmeter  measurements, 
7)  VARA:   Estimated  reading  variance  of  photograph  measure- 
ments, and  8)  F:   Conversion  factor  from  measuring  units  to 
photograph  units. 

CARDS  USED  FOR  SUB-RUNS 

INITIAL  CARD  FOR  EACH  SUB-RUN:   FORMAT (212) 

1)  L:   Number  of  first  current  derivative  peak  entered 

as  data,  and  2)  N:   Number  of  last  current  derivative  peak 

entered  as  data. 

REMAINING  CARDS  FOR  EACH  SUB-RUN:   FORMAT (8F10 . 0) 

A(I) :   Measured  vertical  position  of  I    current  derivative 
peak  as  measured  by  the  traveling  microscope.   (Data  for 
multiple  sub-runs  were  fed  in  order  with  one  set  of  cards 
for  each  of  the  sub-runs.) 

The  output  of  the  program  is  set  up   for  printed  form. 
Values  printed  were:   DN,  DEVD,  VARD,  AMIC,  DEVM,  APH,  DEVP,  B,  AND  BB. 
Some  of  these  values  were  used  for  further  calculations  as  described  in 
Appendix  A. 
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Figure  39.   Logic  diagram  for  current  derivative  calibration  computer 
program. 


Calculate: 

DN,DEVD,VARD, 

AMIC,DEVM, 

APH,DEVP, 

B,BB 


Output: 

DN,DEVD,VARD, 
AMIC,DEVM, 
APH,DEVP, 
B,BB 


f        Start: 


I 


Read: 

T , C , VO , VART , 

VARC,VARV, 

VARA,F 
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APPENDIX  E 

Description  and  Explanation  of  IBM  360/50  Computer  Program  for 
Calculating  Physical  Quantities  Associated  with  an  Exploding  Wire 

This  program  was  written  to  calculate  the  resistance  of  an  explod- 
ing wire,  the  power  input  to  an  exploding  wire,  the  energy  deposited  in 
an  exploding  wire,  the  action  integral  of  an  exploding  wire,  the  voltage 
across  an  exploding  wire,  and  the  current  through  an  exploding  wire. 
Measured  distances  on  photographic  records  of  voltage,  current,  and 
current  derivative  curves  versus  time  were  used  as  data  for  the  calcu- 
lations.  The  theoretical  equations  used  in  the  calculations  were  dis- 
cussed in  Section  3.0  with  the  exception  of  the  action  integral.   The 
action  integral,  defined  as  the  integral  of  the  current  squared  over 

the  time  interval  of  interest,  has  been  used  in  correlating  results  in 

19 
exploding  wire  studies.     However,  in   this  study  the  current  was  not 

uniquely  determined  during  the  initial  portion  of  the  wire  explosion; 
thus  the  action  integral  did  not  give  meaningful  results.   The  equations 
used  for  the  calculations  were  algebraic  except  for  the  calculation  of 
the  energy  deposited  and  the  action  integral.   Simpson's  rule  was  used 
for  the  numerical  integrations  and  the  estimated  integration  error  was 
estimated  as  discussed  in  Section  5.2. 

The  propagation  of  variance  formula  (equation  A-7)  was  applied  (to 
equations  used  to  calculate  the  resistance  and  power)  to  obtain  an  esti- 
mate of  the  variance  of  the  resistance  and  power  due  to  propagation  of 
variance  in  the  experimentally  determined  quantities  (voltage,  current, 
and  current  derivative) .   The  propagation  of  variance  formula  was  applied 
to  the  summation  of  terms  in  the  numerical  integration  by  Simpson's  rule 
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to  determine  the  variance  of  energy  and  action  due  to  the  propagation 
of  variance  of  the  experimentally  measured  quantities. 

These  equations  were  programmed  in  FORTRAN  IV  language  for  the 
IBM  360/50  computer.   The  source  program  is  listed  and  the  logic 
diagram  is  shown  in  this  appendix.   Table  E-l  defines  important  sym- 
bols used  in  the  program. 

The  input  to  the  computer  program  was  arranged  as  follows: 
FIRST  CARD:   FORMAT (5 12) 

1)  N:   Number  of  horizontal  positions  entered  (must  be  odd 
integer),  2)  1A:   Month  of  year,  3)  IB:   Date  of  month,  4)   IC: 
Photograph  number,  and  5)  ID:   Trace  number. 

(IA,  IB,  IC,  and  ID  were  used  to  identify  data  and  relate  data 
to  photographs  and  results.) 
VOLTAGE  DATA:   FORMAT (8F10.0) 

V(J) :   Voltage  parameter  data  which  were  entered  as  the  vertical 
position  of  the  voltage  curve  as  measured  by  the  traveling  micro- 
scope.  Data  were  read  in  eight  values  per  card  for  values  of 
J  =  1,  2,  ...,  N. 
CURRENT  DATA:   FORMAT (8F10.0) 

A(J) :   Current  parameter  data  which  were  entered  as  the  vertical 
position  of  the  current  curve  as  measured  by  the  traveling  micro- 
scope.  Data  were  read  in  eight  values  per  card  for  values  of 
J  =  1,  2,  ...,  N. 
CURRENT  DERIVATIVE  DATA:   FORMAT (SF 10. 0) 

Y(J):   Current  derivative  parameter  data  which  were  entered  as  the 
vertical  position  of  the  current  derivative  curve  as  measured  by 
the  traveling  microscope.   Data  were  read  in  eight  values  per  card 
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Table  E-l:   Symbols  used  in  the  IBM  360/50  computer  program  for 
calculating  physical  quantities  associated  with  an 
exploding  wire. 


Symbol 


Explanation 


AA  Vertical  position  of  current  zeroline  reference 

ACON        Conversion  factor  on  oscilloscope  record  of  current 
curve  (kA/photo  division) 

ACZERO  Action  integral  from  initiation  of  discharge  to  TZERO 

A(J)  Value  associated  with  current  at  horizontal  position  J 

AMAX  Estimated  maximum  action  (Used  in  scaling  for  PLOT) 

CMAX  Estimated  maximum  current  (Used  in  scaling  for  PLOT) 

CONF  Conversion  factor  (.363  microscope  divisions/photo  division) 

DELT        Horizontal  distance  between  successive  points  measured  on 
photograph  (inches) 

DEVA(J)     Estimated  standard  deviation  of  current  at  J   horizontal 
position 

DEVV(J)     Estimated  standard  deviation  of  voltage  at  J   horizontal 
position 

EMAX        Estimated  maximum  energy  (Used  in  scaling  for  PLOT) 

ERROR       Estimated  integration  error  (See  Section  5.2) 

EZERO       Energy  deposited  from  initiation  of  discharge  to  TZERO 

IA,  IB, 

}   Numbers  to  identify  data  and  relate  to  photograph  and  results 

J  Number  used  to  designate  horizontal  position  at  which  curves 

were  measured 

N  Number  of  horizontal  positions  used  to  measure  vertical 

positions  of  voltage,  current,  and  current  derivative  curves 

PMAX        Estimated  maximum  power  (Used  in  scaling  for  PLOT) 

RMAX        Estimated  maximum  resistance  (Used  in  scaling  for  PLOT) 

TCON        Conversion  factor  of  horizontal  scale  of  oscilloscope  record 
(usee/  photo  division) 
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Table  E-l: 

Symbol 

TMAX 

TZERO 

VCON 

V(J) 

VMAX 
W 
WL 
YCON 

Y(J) 


Symbols  used  in  IBM  360/50  computer  program  for  calculat- 
ing physical  quantities  associated  with  an  exploding 
wire  (continued) 

Explanation 

Estimated  maximum  time  (Used  in  scaling  for  PLOT) 

Time  after  initiation  of  discharge  for  initial  data 
point  T(l) 

Conversion  factor  on  oscilloscope  record  of  voltage  curve 
(kV/photo  division) 

Value  associated  with  voltage  at  horizontal  position  J 

Estimated  maximum  voltage  (Used  in  scaling  for  PLOT) 

Vertical  position  of  voltage  zeroline  reference 

Inductance  of  wire  (microhenries) 

Conversion  factor  on  oscilloscope  record  of  current  deriv- 
ative curve  (kA/ysec^photo  division) 

Value  associated  with  current  derivative  at  horizontal 
position  J 


YY 


Vertical  position  of  current  derivative  zeroline  reference 
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for  values  of  J  =  1,  2,  ...,  N. 

CONVERSION  FACTORS  DATA:   FORMAT (8F10. 0) 

1)  DELT:   Horizontal  distance  between  successive  data  points 
measured  with  the  traveling  microscope,  2)  TCON:   Conversion 
factor  of  horizontal  scale  of  oscilloscope  record,  3)  VCON: 
Conversion  factor  of  vertical  scale  of  voltage  curve,  4)  ACON: 
Conversion  factor  of  vertical  scale  of  current  curve,  5)  CONF: 
Conversion  factor  from  microscope  divisions  to  photograph 
divisions,  6)  AA:   Vertical  position  of  current  zeroline  ref- 
erence, 7)  VV:   Vertical  position  of  voltage  zeroline  reference, 
8)  YY:   Vertical  position  of  current  derivative  zeroline  refer- 
ence, 9)  WL:   Inductance  of  the  wire,  and  10)  YCON:   Conversion 
factor  of  vertical  scale  of  current  derivative  curve.   Data  were 
read  in  with  eight  values  on  the  first  card  after  the  current 
derivative  data  and  two  values  on  the  second  card  after  the  cur- 
rent derivative  data.) 

REFERENCE  POSITION  OF  INITIAL  DATA  POINT:   FORMAT (8F10 . 0) 

1)  TZERO:   Time  after  initiation  of  discharge  for  initial  data 
point,  2)  EZERO:   Energy  deposited  from  initiation  of  discharge  to 
TZERO,  and  3)  ACZERO:   Action  integral  from  initiation  of  discharge 
to  TZERO. 

SCALE  FACTORS  FOR  SUBROUTINE  PLOT:   FORMAT (8F10.0) 

1)  TMAX:   Estimated  maximum  value  of  time  -  used  to  scale  hori- 
zontal axis  of  plots,  2)  RMAX:   Estimated  maximum  value  of  resis- 
tance calculated  -  used  to  scale  vertical  axis  of  resistance  ver- 
sus time  plot,  3)  PMAX:   Estimated  maximum  value  of  power,  4)  EMAX: 
Estimated  maximum  value  of  energy,  5)  AMAX:   Estimated  maximum 


149 


value  of  action,  6)  VIIAX:   Estimated  maximum  value  of  voltage, 
and  7)  CMAX:   Estimated  maximum  value  of  current. 
TITLE  AND  AXIS  LABELS  FOR  SUBROUTINE  PLOT: 

CARD  1:   Title  card  for  graphs,  and  CARDS  2-7:   Axis  label  cards 
for  vertical  and  horizontal  axis.   The  horizontal  axis  for  each 
plot  was  time.   The  vertical  axis  were  in  the  following  order: 
resistance,  power,  energy,  action,  voltage,  and  current.   More 

detailed  information  of  the  content  of  title  and  axis  label  cards 

44 
can  be  found  with  write  up  for  subroutine  PLOT. 

The  output  of  the  program  was  obtained  both  in  printed  form  and 
graphs.   The  subroutine  PLOT,  available  for  use  with  the  360/50  computer, 
was  used  to  plot  the  results  in  graphical  form.   The  printed  form  for 
resistance,  power,  and  energy  was  arranged  in  columns  as  follows:   J, 
T(J),  physical  quantity  (neglecting  the  induced  voltage  across  wire), 
estimated  standard  deviation,  estimated  deviation  due  to  measurement 
of  the  photographs,  and  corrected  physical  quantity  (corrected  to  in- 
clude the  induced  voltage  drop  across  the  wire) .   The  printed  form  of 
the  action  integral  was  as  follows:   J,  T(J),  action,  estimated  stand- 
ard deviation,  and  estimated  standard  deviation  due  to  measurement  of 
the  photographs.   The  printed  form  of  the  voltage  and  current  was 
arranged  as  follows:   T,  T(J),  voltage  or  current,  and  estimated  stand- 
ard deviation.   The  printed  values  of  resistance,  power,  energy,  vol- 
tage, and  current  as  functions  of  time  were  used  to  plot  the  results 
shown  in  Figs.  8-33. 
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IF    S  FUNCTION  Or  T  I  ME  ,5Xj  12,  lh 
•:  FUNCTION  Cr  TIMEt5X»I2»lH-»I2 


]   |  J      p  D  /    •        r  ,  p  HS   r;  p 

DLOTTf  ;  •  ■■         \i        THI      FOLLOWING.        (  1  )  ESI        M         *     ( 2 )     PC 

C  (3)     EN^R     ',     (4)     ACTION j     (5)     VOLTAGE,     (6)  r  NT. 

//i         '  iO!       )     3  540412'         ■   .  ]  :  •  ':      •  -    *  »     I  ■  LEVEL=1 

//  ^XrC    FTGCLGKS 

DIMFNcION     V(  ]     n)»A(]     0),T(]     0) »DEVA( 100 ) ,DEVV( IOC) 

p  }  •  •  r  -  v  r  j  /-  • ,     y  ,  t     r, )  ,  \/  a  R  {  ]     P),0'rV(10n)»Yn0O|),RDF 7  (  1  P  0  ) 

o  i     r  o  Q  y  A  t  (  5  I  ?  ) 
101     FORMAT ( 1H1 »15X,4<    ICURREN1     Tl  H    WIKF     AS     FUNCTION    OF     TIME»5X»I2,1 

2H-,  I?,1H-,  I2,1H-»I2) 
10?    FORMAT ( 1H1 ,15X»39HVCLTAi 

2-,I2,iH, I2»1H-, 12) 
1    3    FORMAT  (  1H1  *  1  5X  ,  35i       I  RE      RES  I  ST/  ' 

2>lh-,T2»-lH-,I2) 

I  /,  FORmAt  (  1  HI  »15X,44H ENERGY  DEPOSIT:::  IN  WIRE  AS  FUNCTION  OF  TIME,5X» 

2I2,1H-, 12, 1H-, I2»1H-,I2) 
ire,    FORMAt{  1H1  ,15X»35HACTI0N  INTEGRAL  AS  FUNCTION  OF  T  [ME  ,  SX  ,  I  2  ,  1  H-  ,  I  2 

?,1H-,T?,1H-, I?) 
Pf  FORMAT  (  ?Hl  J^X.i^TI"r("I^"rCr)  *  -  X  ,  1  f  '■'  5  I  STANCE  (OHMS)  »5X»9Hi  E\  " 

IT  I  ON »  ox  »  8  IF  '■ r>  ^F  V  *  S  X  ,  ?  ]  HR  F  c  I  S  T  <= f  !Ci  (  CORRECT  ED  )  ) 
ir7  F0RMAt(  2HLJ»5X  ,14HTIME(MICR0SEC)  ,5X,16i-  NERGY  (  JOULES  )   ,5X,9H0  /I/ 

ITICS-X^H^.^  Dc  V  »  5  X  »  1 7HF  NFRG  Y  {  COI   I  C  T  r  D  )  ) 
i  •  P,  f  r,J,,Ai  (  2riLJ,^X.  14.L  II  I  -!E(  MICROS  EC)  » 5X , 16HACT I  ON  (  COUL-AMP)  »5X»9HDEVI  A 

1 T  I  0 N ,oX»8HREAD  D r  V ) 
ICO  FORN'At  (  ]■■     ,  1 4 , 4X  » F 10 . 5 , 8X , F 14 . 6 ,4X , F 12 . 8 , 4X  » F 12 . 8 ) 
11^  F  OR  WAt  (  1  H  ,  I  4,  4X  »F10 . 5 , 8X , Fl 4 .6  »4X  ,  F]  ;;  •  8  ,  4X  »  F  12  .  8  »4X  ,F  14.  6  ) 
111  FORMAT  (  i  HI  ,  15X»25HPCWER  AS  FUNCTION  OF  T  I  ■"'.   ,  5X  ,  I  2  •  1  I I-  ,  I  2  ,  1  H-  ,  I  2  •  1 

2-,  12) 
11?  -a-''at  (  ?HL  J  ,  5X,  1  4HT  I"r  (MICROSEO)  ,  5X,2]HPOWEF  (  MFGAJOULFc/SEC  )  »3X,oi 

2nF\/IATlON»5X»8HREAn  ^  r  \i ,  5  X  •  ]  6  H I  0  v.'  E  R  (  C  OR  R  E  C  T  E  D )  ) 

I I  a  F Op  M A  t  (  ?  HL  J  ,  5 X  *  1 4H  T  I  v E  (  ,-"  I  0  °  0  c  E  0  )  >  6  X  » ]  5H  VO  L  T  A  G  E (  <  \ '  ^ .  T  S  )  •  5  X  ,  9 H  D  c  V  I  '  T 

?IO?  ! 
1U  F  0°  m  A  t  (  1  H  ,  I  4  ,  F 1  (J :  .  5  ,  8  X  ,  F 1 4  .  6  , 4  X  , F  1 2  .  8  ) 
]]c  FORN'A.t{2HLJ»5X»  14HTIME  (MICRO  EC)»7X»14HC      '.Ti-  »MPS  )  »  5X  »9H  EVIATI 

2  o  N  ) 
H£  FORMAt(  /,i  8HINTEGRA.TI0N  -  RRCR  =  •  -  I '-  .  8  ) 
D  0  3  0  K  =  1  » ]  . 


z  c  ■  -  ( 

R  F  5  D  ( 

9  F  A  D  ( 

p  p  a  r>  ( 

p  c  t  0  ( 
p  p  A  n  ( 
p  r  *  n  ( 

u'p  T  T  F 


,  9  1  )  ,\  ,  I  A  ,  I  R  ,  I  C  »  I  D 
,9  )  (  V  (  J  )  •  J  =  1  »  M  ) 
,9  )  (  A  (  J  )  ,  J  =  1  ,  N  ) 
,°  i)  ( Y( J)  ?J=1»N) 

)DELT,TCON,VCON»ACOf  »COf 
,g  )TZER0*EZER0»AC7EPC 


'  •  .  VV  •>   YV  9  V L  «  VC0M 


^  ,  1  r  ^  )  I  *•  ,  I  r-  •  I  0  ,  I 
3  « 1  C  6  ) 
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P\/ARV  =      .      '/.or      6r  •/.'■'>  yCC  ■  -:    .'         N/l •    ) 

■•  :  R£  -     .       49*::  .00^  '  f  /  (     G    f  *CCNI   ) 

D  r  L  T  =  (  D c  L  T  *  T  C  0  N  )  /  O  N I 
DCK  J  =  ]  » N 
v j=J-i 

T  (  J  )  =  ^FLT*XJ-i-TZERC 
T  P  (  J  )  =  T  (  J  J 

A  (  j  }  =  a  P  <;  (  (  A  A  -  A  (  J  )  )  *  t  c  C  N  /COM!   ) 
Y ( J ) =aRc ( ( ( VV-V ( J ) ) *VCCN- ( YV-Y ( J ) ) *YCCN*VL ) / ( CCMF *A ( J ) ) ) 

Y  P  (  J  )  =  Y  (  J  ) 

\ /  (  j  )  -  a  n ,  c  (  ( \i  v _ \j  (  j  )  )  -* \j (- q ^  /r q m  r  ) 

n  c' V  V  (    |  )=■">.  P049*  VfCN /CCN  F 

prvM   n  =SORT  (  .  00&9*. '  C49*ACCN*ACCN/  (  COI  l  *CCNF)+.0332*.0332*A 
2J  )  ) 
X  (J  )  =  \/  (  J  )  /  A  (  J  ) 
RDEV( j) =  SQR1 ( ( RVARV+ ( RVARA*V ( J ) *V ( J ) ) / ( A( J) -a ( J j ) ) / ( A ( J ) *A ( J 

V  a  R  (  J  )  =  (  DEW  (  J  )  *F>E  yv  (  J  )  +  (  D  EVA  (  J  )  *D EVA  (  J  )  *V  (  J  )  *V  (  J  )  )  /  (  A  (  J  )  *  A  ( 
1  A  (  j  j  *  a  (  j  )  ) 

1 n    DFV( J^ =  SQ9T ( VAR( J] ) 

T  P "  •  A  Y  _  T  V /"■  X 

i,'P  ITPf  3»  1  1  '  )  (  J»T(  J)  •>  y  (  J)  .^r>'  (  J  )  *ROEV( J)  *v  (  J  )  ,  J=l  »N) 


(  J  )  *  /   ! 


)   )  ) 

J  )  )  )  /  ( 


CALL    ^t_CT(TP»'   . »TP^AX,0,YP,0. ,RWAX 


*  *    '       •     S 


•  »o.»n»N»] « Q « n ) 


'.vP  I  TF  f  3  »  1  1  1  )  T  A  ,  I  R  ,  I  r  ,  I  0 

VRI  TFf  3  j  11  2) 

DC] 3  J - 1  » N 

T  P  (  J  )  -  T  (  J  ) 

X  (  J  )  =  A  (  J  )  *  V  (  J  ) 

R  D  E  V  (  J  )  =  S  3  R  T  (  R  V  A  R  V  *  A.  (  J  )  *  A  (  J  )  +  R  V  A  R  A  *  V  (  J  )  *  V  (  J  )  ) 

VAR  (  J  i  =  DE  VV  (  J  )  *DEVV  (  J  )  *A  (  J  )  *a  (  J  )  +  DEVA  (  J  )  *  DEV  "  (  J  )  --  V  (  J  )  * V  (  J  ) 

DFV( Ji=SGRT( VAR ( J ) ) 

Y  (  J  )=v(  J)  *A  (  J)  *A  (  J) 

Y  P  (  J  )  =  Y  (  J  ) 
i-3    r  ^  m  j  i  m  !  i  f 

T  p '  '  A y  _  t  >  i A,  y 

u'PI  TFr  3»11      )  (J^T(J)  ,  X  (  J)  ,DEV(  J)  »ROE  V(  J)»Y(J)»J=1»N) 


TALL    PLCT(TP»  -  .  ,TD,-'AX 


n.vo 


•     •  » F    iAX>u»n»5    ■•».•>■ 1  ?   :  «  I  •  ^  j 


ypiTFf  3  »1CA.)  I  A  »  IB*  IC  j  ID 
WRI  TE(3»107) 
X(l  )=c'ZFRC 

PRVAR  =  -  (  2  .  *DELT*DF  LT/9  .  ) * ( A ( 1 ) *A (  1  )  *  RVA  RV  +  V  (  1  )  *V  (  1  )  *RV  -'  RA  ) 
VAR ( j  ) =-( 2.*DELT*0ELT/9. )* (  !  (  1 )  * A (  1  )*DEVV(  1  )*DEVV (  1  ) +V (  1 ) *V ( 
1  '  (  1  )  *  ^  E  V  \  (  1  )  ) 
FRPCR=(4.*(Y(1)-Y(N)  )  +  ( Y ( 2 ) + Y ( MN )  !  )  * ( DE  LT /9  G •  J 
f,r:Tc  =  /  N-]  )  /2 
DrilJ=3»N»2 


1  )  *  D  E 


JJJ=J-2 

JP= ( J-] ) /? 

rpRQR_cRp3^^  ( P . * 0 E I  T / 9 0 ■ 


)*  ( Y ( J)-Y (  JJ  )  ) 
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X  (  J  )  =  t      :  I  t  /  3  .  )  •;:  (  <\  (  J  J  J  )  x  v  (  J  J  J  )  +4  .  :;  A  (  J  J  )  *  V  (  J  J  )  +  A  (  J  )  *  V  <  J  )  )  +X  (  J  J  J  ) 

Y  (  J  )=v(  J)-(      .  r  "  •  '.    '   '  (  J  !  */  (  J  )  ) 
VD( jp V=Y(     I) 

TP  (  JD  ,rT[     I  ) 

r,nP\/(    ,)  =qqpt  (  (ric-LT*DFLT/Q.  )*{  (  ?  .#A  (  J  J  J  )  *A  (  J  J  J  )  *r      :  ?V  +  3  .  *  V  (  JJJ  )  *V  (  J 
?JJl  '     ft  )  +  (  i^.*/1  (  JJ  )*A  ( JJ  )*PV*     V+16.*V(  JJ)*V(  JJ  }*      ''     M  +  (£  ( J)*A(  J 

■5  )  ^oyaoV-f  V  (  J  )  >  V  (  J  )  *  RV  '  ';  '■  )  )+'  '  P  ) 

PPVAR-ROFVt J ) *RDFV( J  ) 

.  "  R  (  J^  nELT*RFLT/9.  )*  (  (3.-:  (  JJJ)*A(  JJJ)     [  :  VV  (  JJJ)*  )EVV  (  J  J  J  )  +  3.*V  < 
i  J  JJ  )  *\/(  JJJ  )   •  DFVA  (  JJJ  i     [}!  VA  (  JJJ  )  )  +  (  16. *A  (  JJ  ) ":   "  (  J  J  )  '•■   :':  VV  (  JJ  )  *  DEW  (  J  J 
?  )  +1 6 .  -"-V  (  J  J  )  *  V  (  J  J  )  *D!  v"  (  J  J  )  *  "  :  V  \  (  J  J  )  )  +  (  A  (  J  )  * A ( J ) *DE  VV  (  J  )  *  D  i  V  V  (  J  )  +  \  ( 
3J )*V(    I)  *DEVA  (  J  )  *DF  VA  (  J  1  )  )  f  VAR  (  JJJ  ) 
]  1     DEV  (  J  ^  =  .SQRT  (  VAR  (  J)  ) 

D  C  2   U  =  3  » N  >  2 

».'P  T  TP  f  3  » 11     )  (  J  »  T  (  J  )  »  X  (  J  )  » D E  V  (  J  )  » I     »EV  (  J  )  » Y  (  J  )  J 
p  r     f  f»Mf  J  mI  If 

'■°  !  Tr  f  3  »'l  1  6  )  rr?PpP 
TP'-'  A  X-T  VIA  V 

CALl      r>'  CT  (  TP  ,  ■    .  ,TP^AX»0»YP».'i.  ,pv|AX»C»^.  »0«  ,     .»0»N»1»2«01 

1»'R  ITF(3»105)IA,IP,IC»ID 

'.•-'RITF(3»1CS  ) 

X  (  1  )  =  a  C  Z  E  R  C 

\tpjS=  (N'-l  )  /2 

R  R  V  ••  ■  p  r  -  (  *  .     *  D  E  L  T*DE  L  T  /  9  .  )  *  A  (  1  )  *  A  (  1  )  *R  VA  R  A  ) 

VAR  (3  )=-(8. J*DELT*0ELT/9.  )*(  U  1  )*A(  1)*   >i  V-   (  l)*l   rvA(  1  )  ) 

0R= (DELT/9     . )*( 4.*( A ( 1 ) -A ■  (N  )  )+A(2) +A ( NN )  ) 
RC3  2J  =  39\>9? 
JJ=J-i 
JJJ=J-2 
JP= ( J-i )/? 
TP( JP^=T( J) 

ERRCR=ERRCP+(8.*DELT/9C.)*(  -'  (  J  )-A(  JJ)  ) 
X  (  J  )  =  i  HFLT/3.  )  *  (  A  (  JJJ  )  *  A  (JJJ  )  +4.  *  "-  (  J  J  )  *A  (  JJ  )  +A(  J  )  *A  (  J)  )  +X  (  JJJ  ) 

Y  P  (  J  P  )  =  X  (  J  ) 

EV(    I)  =SQRT(  (4.*RVARA*DELT*DELT/9.  )*  (3.*A(  JJJ)  *A  (  JJJ  ) +16.*  A  (  JJ  )  *  •' 
2( JJ)+a ( J) *A( j) ) +RRVAR ) 
o  q y ;  r  =  p r) E  V  (  J  )  *  R D  E  V  (  J  ) 

VAR  (  J\  =  (A..*rc|_T*-nFl_T/9.  )*(  *.*A(  JJJ  )*A  (  JJJ)*r»PVA  (  JJJ)  *r^VA  (JJJ)  +  3  6. 
1  -  M  JJ«  «  A  (  JJ)  *i>FV«  (  J  J  )^rV^  (  J  J  )+A ( J)  *A  ( J)  *OEVA  (  J)-r  rV-A  (  J)  )  +VAR  (JJJ  ) 
1  ?    pf\/(  j>  =  <;opt  {  VAP  (  J)  ) 
0C21J-3  9M  9  2 
i,.c?  t  rr  r  3  9  3   %Q  )  (  J  9  T  (  J  )  9  X  (  J  )  9 OEV  (  J  )  9  RHEV ( J )  ) 

VPlTFf 3  9ll6)ERRCR 
TP"iY_T'/AX 

C^LL     DLCT  (  TP91    .  «TDV'AX  ,  C  9YP  9'    .  *  "■'■'■'  X 
■WRI TE(3»1C2)IA»IP»IC»ID 
W  R  I  T  E  (  3  9  1 1 3  ) 
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?  ?      T  D  {  J  )  -  T  (  J  ) 

TDVAy-JV' x 

>/.'R  I  T  F  {  3  ,  1  ]  4  )  (  J  ,  T  (  J  )  »  V  (  J  )  »  0  F  V  V  (  J  )  »  J  =  1  >  N  ) 

CALLnLCT(TP»     .  »  T I       :      •     » V  *  0  .  » V?>  A  X  ,  C  ,     .  ,  '   .  , r  .  ,  ':  .  ' .  ,  1  ,  2  ,  Q  ) 
D  0  2  3  J  =  1  >  N 
2  ^    T  P  (  J  )  =  T  (  J  ) 
yo:-A  y-T"  '  X 

«/9ITF.  f  3  » 1  C 1  )  IA,  18,  IC»  ID 
'•■,£?  I  T  p  (  3  » 1 1  5  ) 
i.'r  i  j  F  f  3  , 1  ]  u  )(J,T(J)»A(J)»DP  V  -  (  J  )  »  J  =  1 »  N  ) 


CALL    n\  CT  (  T P  , '  .  ,  TP''  A  X  ,  P  ,  A  ,  n  .  ,  r>' 

■an     rA'.|T  n-mp 
c  j  q  p 

/* 


) 
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Figure  40.   Logic  diagram  for  physical  quantities  of  an  exploding 
wire  computer  program. 


Start 
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Resistance 
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Calculate: 
Power 
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Calculate: 
Resistance 


Output: 
Power 
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Conversion  factors 
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Calculate: 
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J 


Convert: 
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voltage,  current,  and 
current  derivative 


Output: 

Energy 
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APPENDIX  F 
Energy  Losses  by  Conduction  and  Convection  from  Exploding  Wires 

An  upper  limit  of  the  heat  losses  by  conduction  and  convection 

from  an  exploding  wire  can  be  obtained  as  described  below.   Since  the 

wire  length  is  much  greater  than  the  diameter,  the  heat  losses  through 

the  ends  may  be  neglected  and  the  wire  considered  as  an  infinitely 

long  cylinder.   Assuming  there  is  a  step  change  in  temperature  from 

ambient  conditions  to  the  maximum  temperature  (T   )  and  further  assuming 

r  max  ° 

the  wire  surface  temperature  (T  )  is  equal  to  the  temperature  of  the 
ambient  medium  gives  an  upper  limiting  value  for  heat  losses  due  to  con- 
duction and  convection  for  the  time  interval  from  when  the  wire  initially 
begins  to  be  heated  until  the  wire  explosion. 

The  transient  temperature  distribution  in  the  wire  is  obtained 
by  solving  the  heat  diffusion  equation, 

(  H  }   =   av2  T  C*'1) 

where: 

T  =  T(r,t)  is  the  time  dependent  temperature  distribution, 

a  is  the  thermal  diffusivity,  and 

r   is  the  initial  wire  radius; 
o  ' 

with  the  boundary  conditions 

T(r,co)  =  T  , 

00 

T(r   t)  -  T  , 

O  co' 

T(r,0)   =  T     ,  and 

max 

T(0,t)   -  finite. 
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Introducing  a  change  of  variables  such  that 

6(r,t)   =  T(r,t)   -  Tw,  (F-2) 

and  substituting  equation  (F-2)  into  equation  (F-l)  gives 

_39(r,tX  .  o[  1    »  (r  |£  j   _  (F_3) 


3t         L  r   9r   v   3r' 

The  boundary  conditions  are 

8(r,°°)   =  0  ,  (F-3a) 

6(ro,t)  =   0  ,  (F-3b) 

6(r,0)   =  T(r,0)  -  T^  =   0    ,  and  (F-3c) 

6(0, t)   -  finite.  (F-3d) 

Equation  (F-3)  can  be  solved  by  separation  of  variables  using  the 

substitution 

6(r,t)   =  R(r)x(t).  (F-4) 

Substituting  equation  (F-4)  into  equation  (F-3)  and  separating  the 
variables  gives  the  following  differential  equations: 


~^-   +   aA2T(t)   =   0   ,  and  (F-5) 


13   [  r  i|(ri  j   +  A2   }  m     Q  _ 

r   dr       3r 

Equation  (F-6)  is  Bessel's  equation  of  order  zero  with  parameter   A. 
With  the  boundary  condition  (F-3d)   the  solution  of  equation  (F-6)  is 


R(r)   =    I  b   J  (A  r)  (F-7) 

*%     n  o  n 
n=l 
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where: 

J_(^  r)  is  the  Bessel  function  of  the  first  kind  of  order  zero. 
(J  n 

The  solution  of  equation  (F-5)  is 

T(t)   =  a   exp(-aA2t)  (F-8) 

n 

where  the  boundary  condition  (F-3a)  has  been  used.   By  substituting 

equation  (F-7)  and  (F-8)  into  equation  (F-4)  gives 

2 
6(r,t)   =    )   K   exp(-aX   t)   J  (A  r)   .  (F-9) 

L-       n    c         n  '   o  n 
n=l 

Applying  boundary  condition  (F-3b)  gives  that  A  is  the  n   root  of 

J  (A  r  )  =  0.   To  evaluate  the  constants  K  the  solution  is  written  at 
o  n  o  n 

time  t  -  0  (Using  F-3c) . 


6(r,0)   =0   =    Y   K  J  (A  r).  (F-10) 

o      '•-non 
n=l 


From  the  orthogonality  relations  the  constants  K  are  given  by 

26 
o 

K   = r-r^ r  (F-ll) 

n       A  r  J-  (A  r  ) 

n   o   1   n  o 


where  J, (A  r  )  is  the  Bessel  function  of  the  first  kind  of  order  one, 
1  n  o 

The  temperature  distribution  is  given  by 

e(r,t)  =  I        [  °  i    ]      Jo(Xnr)  exp^aA^  t)         (F-12) 

n=l     n  o  1  no 

The  average  temperature  rise  of  the  wire  above  ambient  conditions  at 

time  t  is  given  by 

[      o   8(r,t)  rdr 
8  (t)   =   J  Q  .  (F-13) 

rr 
o 

rdr 
o 
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By  substitution  and  integration  we  obtain 


40n  9 

6(t)   =   I     -z-5 exp(-aXZ  t)  ,  (F-14) 

n-1  *n  rQ2         * 


which  is  the  average  temperature  rise  above  ambient  conditions  at 

time  t. 

The  ratio  of  G(t)  to  6   gives  an  estimate  of  the  maximum 

o 

energy  loss  by  conduction  and  convection  at  time  t  for  an  exploding 
wire.   Representative  values  of  0(t  )/9   are  given  in  Table  F-l, 

where  t   is  the  maximum  time  interval  of  interest  during  a  wire  ex- 

m  ° 

plosion.   For  wires  which  exploded  on  the  first  current  pulse,  t 

was  taken  to  be  the  time  of  burst.   For  wires  which  exploded  on  the 

second  current  pulse,  t  was  taken  to  be  the  time  from  when  the 

m 

voltage  crosses  the  zeroline  to  the  time  of  burst. 


Table  F-l:   Estimates  of  maximum  energy  losses  from  exploding  wires 
by  conduction  and  convection. 

Wire  Material     Wire  Radius         a  t™      0(t  )/0 

m   o 


Aluminum  .0160  .232  10  .79 

Copper  .005  .189         2.5  .71 

Copper  .010  .189         3.0  .84 

Iron  .007  .0328        2.5  .91 


Wire  Radius 

a 

Cm 

(in.) 

2 
(in. /sec) 

(Vsei 

.0160 

.232 

10 

.005 

.189 

2.5 

.010 

.189 

3.0 

.007 

.0328 

2.5 
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APPENDIX  G 

Skin  Effect  and  Magnetic  Pressure 

Since  the.  time  to  burst  of  an  exploding  wire  is  characterized 
by  times  in  the  microsecond  range  or  less  and  by  high  currents  which 
may  change  rapidly,  the  skin  effect  may  be  important  and  must  be  con- 
sidered.  The  application  of  a  step  function  increase  in  current  to 

a  cylindrical  wire  has  been  considered  to  obtain  a  measure  of  the 

27 
transient  skin  effect.    The  ratio  of  the  current  at  the  wire  center 

to  the  steady  state  current  (the  current  in  the  wire  when  a  uniform 

density  of  current  is  reached)  is  given  by 

i    (r=0)  exp    (-c*2    t/r§au)  (r    n 

$  =     _^ =     1     +  — j— '  V*~L> 

l  J    (a   ) 

00  on 


where: 


a     is    the  n        root   of  J.,  (a)    =   0, 
n  1 


a  is  the  electrical  conductivity, 

y  is  the  permeability  (4tt  x  10   H/m) ,  and 


r   is  the  wire  radius. 
o 


Equation  (G-l)  was  used  to  calculate  time  intervals  at  which  the 
skin  effect  would  be  important;  the  calculated  results  for  wires  of 
interest  in  this  investigation  are  given  in  Table  G-l.   Time  intervals 
longer  than  the  reported  values  would  have   $>.95.   Since  the  times  to 
burst  for  wire  explosions  in  this  study  were  greater  than  the  time  inter- 
vals in  Table  G-l,  consideration  of  the  skin  effect  was  not  necessary  in 
this  study. 
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Table  G-l:   Time  intervals  at  which  the  ratio  of  the  current  at  a  wire 
center  to  the  steady  state  current  equals  .95  for  wires 
used  in  this  study. 


Wire 

Wire 

Material 

Radius 
(m  x  104) 

Aluminum 

4.06 

Copper 

1.27 

Copper 

2.54 

Iron 

1.78 

(mho-m   xlO 

-7) 

t 
(ysec) 

3.8 

2.09 

5.8 

.31 

5.8 

1.25 

1.0 

.11 

At  high  current  densities  it  is  expected  that  magnetic  pressure 
will  have  a  significant  effect  on  a  wire  explosion.   If  the  self- 
magnetic  field  is  to  compress  a  wire  without  ejecting  it  from  the 
circuit,  it  is  necessary  that  the  return  current  be  symmetrical  about 

-I  £ 

the  wire.     Figure  41  shows  an  enlarged  drawing  of  the  sample  cell 

used  in  this  study.   Since  the  return  current  was  symmetrical  about 

the  wire   the  magnetic  force  was  directed  towards  the  center  of  the 

wire  and  the  magnetic  pressure  would  tend  to  confine  the  wire  material, 

When  the  current  flowing  in  the  wire  is  uniform  throughout  the 

wire,  the  forces  on  the  moving  charges  result  in  a  pressure  distri- 

bution  (dyne/cm  )  through  the  cross  section  of  the  wire  given  by. 

2  2 

p(r)   =  ~Y~2 (  1   -  ~  )  (G-2) 

uc  r  r 

o  o 

where: 

I  is  the  current  flowing  in  the  wire  (statamps) , 
c  is  the  velocity  of  light   (cm/sec) ,  and 

r  is  the  distance  from  the  center  of  the  wire  to  the  point  where 
the  pressure  is  evaluated   (cm) . 
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Figure  41.   Sketch  of  Magnetic  Pinch  Sample  Cell 
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If  the  current  is  limited  to  a  very  thin  layer  on  the  surface  of 

the  wire,  i.e.,  skin  effects  dominate,  the  pressure  is  constant  inside 

45 
the  cylinder  and  equal  to 

p   -   (  l2/2Tfr2  c2  )   .  (G-3) 

o 

Note  that  the  pressure  calculated  by  equation  (G-3)  is  one  half  the 

pressure  calculated  for  an  axial  position  of  the  wire  using  equation 

(G-2). 

Calculating  the  pressure  for  a  10  mil  diameter  wire  (r   =  5  mils) 

at  the  axial  position  using  equation  (G-2)  gives  the  following  results 

Current  Density  Pressure 

2 
(A/cm  )  (atm) 

108  50,000 

107  500 

106  5 
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Abstract 

Exploding  wire  experimental  results  are  given  in  terms  of  voltage, 
current,  resistance,  power,  and  deposited  energy  as  a  function  of  time 
after  the  initiation  of  an  explosion.   Voltage  and  current  measurements 
were  made  for  wires  at  atmospheric  conditions  as  well  as  wires  exploded 
in  a  sample  cell;  this  cell  tended  to  confine  the  wire  material  and  to 
increase  energy  deposit  by  employing  the  magnetic  pinch  principle.   Cir- 
cuit parameters  varied   included   capacitance  (4  and  8  microfarad 
capacitor  banks  were  used),  and  charging  voltage  (8,11,  and  15  kilovolt 
charging  voltages  were  studied).   Wire  parameters  varied  included  wire 
material  (aluminum,  copper,  iron,  and  nickel-chromium  wires  were  ex- 
ploded) ,  wire  diameter  (10-32  mil  diameters  were  used) ,  and  wire  length 
(1/2  -  3  inch  wires  were  studied) .   Experimental  results  show  that  rapid 
matching  of  wire  parameters  to  circuit  parameters  (to  obtain  the  most 
energetic  wire  explosions)  is  possible  knowing  the  qualitative  effects 
of  the  wire  parameters  on  the  deposited  energy.   It  has  been  shown  ex- 
perimentally that  about  six  times  the  desired  energy  deposit  must  actually 
be  stored  in  the  capacitor  bank  to  produce  the  temperatures  or  energy 
densities  predicted  by  the  theoretical  design  method  presented.   Further 
work  required  to  simulate  the  chemical  environment  of  a  nuclear  explosion 
by  exploding  wire  techniques  is  discussed. 


